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ON  THE  RESONANCE  CONCEPT  IN  SYSTBIS  OF  LINEAR  AND 
NONLINEAR  ORDINARY  DIFFHiaJTlAL  EQUATIONS 

Rahird  Ibrahim  Ibrahim  Abdel  Karim 


Theorems  on  the  resonance  cases  for  linear  and  nonlinear 
ordinary  differential  equations  of  the  first  to  the  n^*^ 
order  are  derived  and  proved  in  detail,  using  an  earlier 
report  by  the  same  author  as  partial  basis.  Minimal 
orders  of  magnitude  of  the  solutions  and  their  derivatives 
are  given  and  methods  for  the  formation  of  examples,  with 
sanple  calculations  in  matrix  notation,  are  described. 
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PART  I 


THE  RESONANCE  CONCEPT  IN  SYSTEJiS  OF  n  LINEAR  ORDINARY 
DIFFEREKHAL  EQUATIONS  OF  THE  FIRST  ORDER 

Section  1.  Problem  Fonnxilation,  Principal  Results 


R.Iglisch  investigated  the  resonance  concept  in  linear  ordinary  differen¬ 
tial  equations  of  the  second  order  (Bibl.5)«  These  considerations  will  be  ex¬ 
tended  here  to  systems  of  n  ordinary  linear  differential  equations  of  the  first 
order  which  will  be  written  in  matrix  form  with  r(t)  as  the  sought  vector: 

J[  =  «(/)r  +  f(/);  (1) 

where  the  square  (for  exan?)le,  real)  matrix  2((t)  is  to  be  continuous  in  t  and 
periodic  with  the  period  P,  i.e.,  all  n®  elements  auCt)  are  continuous  functions 
in  t,  periodic  with  P;  the  (for  exanple,  also  real)  vector  f(t)  is  assumed  as 
also  being  continuous  and  periodic  with  P: 

«(/  +  P)  =«(/);  f(/  +  P)  =  f(/).  (2) 

The  homogeneous  system  conjugate  to  eq.(l)  reads 


while  the  "adjoint"  system  is 


(3) 

(A) 


where  the  superscript  T  is  to  denote  the  transition  to  the  transposed  matrix. 

Definition  1.  In  the  inhomogeneous  differential  equation  syston  (l),  the 
resonance  case  is  present  if  the  adjoint  systen  {if.)  has  at  least  one  solution 
vector  3(t)  periodic  with  P,  for  which 

Ji^(D\(i)it  =  c=^o  (5) 
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is  valid 


In  Section  4,  we  will  prove: 

Theoran  1»  In  the  resonance  case,  any  solution  vector  ^(t)  of  eq.(l) 
assumes  arbitrarily  large  values,  increasing  without  bounds  with  increasing  t. 

Definition  2«  In  eq.(l),  we  have  the  principal  case  if  eq.(4)  has  no 
solution  vector  ^(t)  periodic  with  P. 

In  Section  5»  we  will  prove: 

Theoran  2»  In  the  principal  case,  solutions  |'(t)  of  eq.(l),  remaining 
restricted  for  all  t,  are  in  existence,  for  example,  the  uniquely  existing 
solution  periodic  with  P. 

Definition  3«  Equation  (l)  represents  the  exceptional  case  if  eq.(4)  does 
have  solution  vectors  ^(t)  periodic  with  P  but  if  the  following  relation  applies 
to  all  these  solutions  periodic  with  P: 

/Vwf(0-*<  =  o.  (6) 

In  Section  6,  we  will  prove: 

Theorem  3«  EJven  in  the  exceptional  case,  restricted  solutions  f(t)  of 
eq.(l)  exist  for  all  values  of  t,  for  exanple,  solutions  periodic  with  P  but  no 
longer  uniquely  defined. 

Sections  2  and  3  contain  auxiliary  considerations  on  the  inhomogeneous 
system  (l),  specifically  on  the  correlation  between  the  homoga:ieous  systan  (3) 
and  the  adjoint  system  (4).  In  this  case,  the  periodicity  stipulation  (2)  will 
be  introduced  only  in  Section  3» 

Many  of  the  results  are  already  known,  but  they  are  here  derived  in  a 
manner  that  requires  no  specialized  knowledge. 
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Section  2.  Auxiliary  Considerations  on  Systems  of  Linear  Differential 
Equations  of  the  First  Order 


In  this  Section,  the  periodicity  stipulation  (2)  will  not  be  used. 

Let  hi(t),  hnCt)  be  a  linear  independent  solution  syston  (fundamental 

syston)  of  eq.(3)  at  the  point  to,  which  can  be  combined  into  the  solution 
matrix 

?){/)  =  (|),W . 0.W) 

We  then  have 

?){<•)  =4=  0. 

Theorem  A.  From  eq.(8)  it  follows  that 

V(0  =  Det3(/)4=O 


(7) 

(8) 

(9) 


for  all  values  of  t. 

Proof.  If  the  conjugates  [subdeterminants  with  correct  sign  (-1)  ]  for 

the  element  yiv(t)  in  eqs.(7)  or  (9)  are  denoted  by  Yiv(t)  and  putting 


0  for  $4=*, 
1  for  i  —  k. 


(10) 


the  following  result  will  be  obtained,  taking  eq.(3)  into  consideration  and 
denoting  the  differentiation  to  t  by  a  prime: 


rfV 

di 


=  Z  y:.  v;.  ==  Z y*.  V.. = Z  1 


r,  i.A 


and  thus 


y{l)=  V'(gexp(/ (a,,  +  a„+ ••• 


Z2 


From  this,  theorem  U  follows  directly. 

Theorem  5.  In  addition  to  eq.(7),  the  expression 


(11) 


at  an  arbitrary  constant  matrix  5  with  a  determinant  C  differing  from  zero, 
represents  a  fundamental  system  of  solutions  of  eq.(3). 
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Proof*  Equations  (7)  and  (3)»  together  with 


=  . «C) 

can  be  combined  into 

^  =*WD. 

From  this*  we  can  calculate 

(f)(1)' =  ?)'«=«  DC  =  *(9tt). 

This  means  that  each  fS  is  a  solution  system  of  eq.(3)» 
That  we  also  have  det(?®)  f  0,  follows  from 


det{f)C)  =  y.C+0. 

Theorem  6.  The  vectors 


V(<) 


»i=  I.....*. 


which  can  be  combined  into  the  matrix 


form  a  fundamental  system  of  solutions  of  eq.(4). 
Proof*  Equation  (15)  indicates  directly  that 


(12) 

(13) 


(U) 


(15) 


(16) 

i*e*,  that  it  is  the  lonit  matrix*  On  introducing  the  reciprocal  matrix  ,  we 
can  then  also  write 

=  3  (17) 

From  eqs*(l6)  or  (17)  it  follows  directly  that 

Z  =  det3  +  0. 

The  statement  that  3  is  the  solution  matrix  of  eq,{l^.)  is  derived  as  follows; 
According  to  eq*(l6),  we  have 

3’‘'f)  +  3*^f)'  =  o. 
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l.e>»  with  eq«(13} 


80  that 


-3"  =  3''«?)9‘  =  3’’*. 


3'  =  -«*'3. 


which  coincides  with  eq,{u)» 


Thus,  theorems  5  and  6  yield  directly: 

Theorem  7.  A  general  fundamental  system  of  solutions  of  eq.(iv)  is  repre¬ 


sented  by  the  matrix 


3*  =  3  «  =  (»-’)’■  ff 


at  arbitrary  constant  matrix  S  with  a  determinant  C  differing  from  zero. 
Denoting 

the  following  is  valid  for  two  arbitrary  vectors  h(t)  and  3(t): 

j’‘L(t)) -(!•(*))' 9=  (: 

and 

Therefore,  the  (expanded)  Lagrange  identity 

is  valid  for  the  solution  vectors  t)  and  j|  of  eqs.(3)  and  (4),  respectively. 
Ana3.ogously,  the  following  is  valid  for  the  solution  vectors  of  eqs.(l) 


and  (4): 


bi¬ 


section  3 .  Systems  of  Linear  Differential  Equations  with 
Periodic  Coefficients 


In  this  Section,  the  periodicity  stipulation  (2)  is  essential. 

Theorem  8.  The  homogeneous  systems  (3)  and  (4)  have  the  same  number  of  p 
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linearly  independent  solutions  periodic  vdth  P  (0<  p  <  n). 

Proof*.  Let  P  be  the  number  of  linearly  independent  solutions  of  eq.(3), 
periodio  with  P.  It  is  merely  necessary  to  demonstrate  that  eq.(4)  has  exactly 
P  linearly  independent  solutions  that  are  periodic  with  P,  since  the  conclusion 
can  be  reversed  in  view  of  the  fact  that  eq.(3)  is  the  adjoint  system  to  eq.(L). 

Case  1.  Let  P  =  n.  Here,  the  entire  matrix  pKt)  is  periodic  with  P  so 
that  also  the  matrix  according  to  eq.(15)  will  be  periodic. 

Case  2.  Let  0  <  p  <  n.  Assume  that  the  hv(t),  with  v  =1,  2,  ...»  P,  /II 
are  periodic  with  P.  Then,  the  following  is  valid  for  any  arbitrary  point  ti : 

1^(01!:"’  +  /<  =  e+i  —  «•  (26) 

If  the  equal  sign  were  present  in  eq.(26),  the  periodicity  of  thishu(t)  with 
the  period  P  would  follow  from  eq.(3)  with  consideration  of  eq.(2).  In  the 
following,  let  ti  be  an  arbitrarily  selected  but  then  retained  point.  Integra¬ 
tion  of  the  first  equation  in  the  system  (24)  over  ti  and  ti  ♦  P  will  yield, 
for  each  k  =  1,  2,  ...,  n, 

hi  ~0  for  v^h.2 . p.  (27) 

These  are  n  linear  homogeneous  equations  which  have  at  least  the  p  linearly  in¬ 
dependent  solutions,  tiv(ti  )  for  v  =  1,  2,  ...,  P.  It  will  be  demonstrated  that 
no  further  solution  vector  Ui ,  linearly  independent  of  this,  can  exist.  Con¬ 
versely,  let  us  assume  that,  for  such  a  ui ,  the  following  is  also  valid 
> 

=  o  (28) 

Then,  let  us  define  a  vector  u(t)  from  eq.(3)i 

-«{/)u  *t«/i  u(/,)  =  u,.  (29) 

at 

*  By  R.Iglisch. 
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For  this  vector  u(t),  the  first  equation  of  the  system  (24),  integration 
over  ti  and  ti  P  for  k  =  1,  2,  •••»  n,  vdll  yield 

(<i  f  P)  u  (<,  f  P)  -  jl  (<«) « (4)  =  0 
and,  after  subtraction  of  eq*(28), 

iir(4  +  P)[u(<,+  P)-u(4)]  =  o. 

This  represents  a  linear  homogeneous  system  of  equations  with  a  determinant 
differing  from  zero;  consequently, 

u(4  +  P)  =  u(4) 

must  be  valid  from  which,  because  of  eqs.(29)  and  (2),  the  periodicity  of  u(t) 
with  P  is  obtained.  Consequently,  because  of  eq.(26),  Ui  =u(ti)  is  linearly 
dependent  on  hi(ti),  ...,  ho(ti). 

Thus,  we  know  that  the  system  of  equations  (27)  has  exactly  0  linearly  in¬ 
dependent  solutions.  Therefore,  the  matrix 

3(0ir  "  =  (hWir  •••.j.wr') 


has  the  reink  n  -  p.  By  suitable  numeration,  it  becomes  possible  that  exactly 
^l(t)!*^*',  ...,  are  linearly  independent.  The  remaining  of  these 

vectors  can  be  linearly  expressed  by  these  vectors; 


Consider  now  the  vectors 


Since 


(i.(0 

i:(0=  V 

liwW  ~  2  ^<*,9 


for  w  =  1 ,  2, . . . ,  —  p 

for  m  =  »«  — p  +  1 . 


(30) 

(31) 
Z12 


Det  (j?{0 *:  (0)  =  Det  (*I  (0 »-  W)  +  0 

it  follows  that  the  n  vectors  (31)  are  linearly  independent  solutions  of  eq.(4) 
for  all  t  (see  also  theorem  4)»  '^^v(t)  with  v=l,  2,  ...,  n  —  P,  because 
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of  the  linear  independence  of  the  vectors 

iTWir" 

including  their  linear  combinations,  are  not  periodic  with  P.  Conversely,  the 
^*(t)  with  u  =  n  -  P  +  1,  n  have  exactly  the  period  P,  -as-readily  seen 

from  eq8.(30),  (4),  and  (2).  This  yields  the  proof  for  our  case  2. 

Case  3«  Let  p  =  0.  If  eq.(4)  had  a  solution  periodic  with  P,  also  eq.(3) 

would  have  such  a  solution  according  to  the  above  statements. 

Section  4.  The  Resonance  Case 

Proof  for  'Eheorem  1.  Thus,  let  ;j(t)  be  a  solution  vector  of  eq.(4)  peri¬ 
odic  with  P,  for  which  eq.(5)  is  valid. 

In  contrast  to  the  argument,  we  make  the  assxmq^tion ; 

\m^B  (32) 

for  all  t  and  for  an  arbitrarily  selected  solution  ^■(t)  of  eq.(l).  An  integra¬ 
tion  of  the  first  equation  in  the  system  (25)  between  t  and  t  mP,  with  arbi¬ 

trary  positive-whole  m  and  taking  eq.(5)  into  consideration,  will  yield 

»’’(<)  (33) 

Because  of  the  periodicity  of  ^^(t),  a  restriction  of  the  following  fom  exists 
for  all  t: 

(34) 

On  the  basis  of  this  and  from  eq.(33)»  an  estimate 

D-2E^m\C\, 

is  obtained.  At  sufficiently  large  m,  this  furnishes  a  contradiction  to  eq.(32). 
The  result  can  also  be  formulated  as  follows: 

Theorqn  9.  In  the  resonance  case,  for  each  solution  vector  r(t)  of  eq.(l) 
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and  for  each  Interval 


at  least  one  point  t  exists  for  which 


vanishes  with  C  from  eq.(5)  and  D  from  eq»(34)» 


(35) 

(36) 


Section  5»  The  Principal  Case 


In  this  Section,  it  is  assumed  that  eq.(4)  has  no  solution  3(t)  periodic 
with  P.  First,  we  will  prove: 

Theoran  10«  If  ^i(t),  •••,  5|B(t’)  are  linearly  independent  solutions  of  /n 
eq«(4)»  we  have 

/irwir'v 

D(/)=Detj  14=0  f0raU  t. 

Virwrv 

Proof*  If,  for  a  special  qioantity  t  =  tg,  we  would  have  D(tc  )  =  0  in 
contrast  to  the  argument,  then  n  numbers  org,  or,  with  orf  +  ,,,  +Q'f  >0 

could  be  so  determined  that 


Putting 


.-I  »-• 


we  would  have 


j,(/, +  P)  =  *,(<»)■ 


Ifowever,  according  to  eqs.(4)  and  (2),  it  would  follow  that 

i,(<  +  P)=».W 

for  all  t,  in  contradiction  to  the  assun?5tion* 

Proof  of  Theoran  2»  If  J■(t)  were  a  solution  of  eq.(l)  periodic  with  P,  an 
integration  of  the  first  equation  in  the  system  (25)  between  a  fixed  point  to 
and  ts  +  P  would  yield 
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(38) 


if  (Oli* '’*{<.)=/ *7 for  I- =1,2 

I. 

Since  the  coefficient  determinamt  of  this  systan  of  equations,  according  to 
theorem  10,  differs  from  zero,  the  quantity  fCts )  can  be  uniquely  determined 
from  this.  Let  us  assume  that  eq.(l),  with  these  initial  values  at  the  point  ta 
is  solved,  thus  yielding  a  solution  vector  p(t).  It  only  remains  to  be  proved 
that  this  r(t)  is  periodic  with  P.  In  view  of  eq.(2),  it  will  then  merely  be 
necessary  to  determine  that,  automatically, 

f(<.  +  P)  =  iW  ^^9) 

An  integration  of  the  first  equation  of  the  system  (25)  over  ta  and  ta  +  P  will 

yield,  for  v  =  1,  2,  ...,  n, 

jf  (/,  +  P)  ■  £(<,  +  P)  -  if  W  ■  I W  =  /if  (0 f  (0 dt. 

t. 

On  deducting  eq.(38)  from  this,  the  linear  system  of  equations  will  read 

with  a  determinant  differing  from  zero,  analogous  to  eq.(8).  This  leads  to 
varnishing  of  the  bracket  and  thus  of  eq.(39). 

Section  6.  The  Exceptional  Case  /L 

If  eq. (4)  has  exactly  0  linearly  independent  solutions  ^i(t),  ^(,(t) 


periodic  with  P,  for  which 


r=1,2 . e 


is  valid,  then  the  determinant  (37)  will  exactly  have  the  rank  n  -  p.  Of  the 
system  of  eqi;iations  (38),  the  first  0  equations  (for  v  =1,  2,  ...,  o)  are 
automatically  satisfied  since  they  contain  only  zeros.  From  the  remaining 


equations  (38)  vdth  v  =  p  +  1, 


n,  a  total  of  p  linearly  independent  vectors 


rfW.rfW . j;(4)  (42) 

can  be  determined,  supplemented  by  the  corresponding  initial  values  (42)  by 
solving  eq.(l): 

?.•(<).  jfw . irw- 

.y. 

That,  for  each  of  these  ru(t)  (u  =  1,  2,  .•*,  0), 

(44) 


is  automatically  obtained  follows  from  the  system  of  equations  which  is  derived 
analogous  to  eq.(40): 

iT  (4  +  P)  [j:{4  +  P)-  tfW]  =  0  ^  ^ 

with  a  determinant  differing  from  zero.  This  means  that  all  these  r(r(t)  will 
have  the  period  P.  Thus,  we  can  make  the  following  statement: 

Theorem  11.  If  eq.(4)  has  exactly  p  linearly  independent  solutions  peri¬ 
odic  with  P,  for  each  of  which  eq.(6)  is  valid,  then  eq.(l)  has  a  P-parametric 
family  of  solutions  periodic  with  P. 

This  resxilt  agrees  with  the  trivial  fact  that  all  solution  vectors  of 
eq.(l)  periodic  with  P  are  obtained  by  adding  to  one  of  these  vectors  all  solu¬ 
tion  vectors  of  eq.(3)  periodic  with  P. 
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PART  II 


m 

THE  RESONANCE  CASE  IN  SYSTE21S  OF  n  NONLINEAR  ORDINARY 
DIFFERQJTIAL  EQUATIONS  OF  THE  FIRST  ORDER 

Section  1.  Problem  Formv>lation,  Principal  Results 

As  an  extension  of  investigations  made  liqr  R.Iglisch  (Bibl.l),  we  consider 
here  the  nonlinear  differential  eq^lation  systan 

. •*.;<)  +  *.(<)  {i  =  i.2.....n)  (1) 

which,  \mder  introduction  of  the  vectors 

(*i(<)\  /?i(“i . 

**Wj  and  9(u,<)=r*^“”  (2) 

*»(<)/  '«.(“! . •'..0/ 

can  be  written  in  the  following  form: 

4“  =  + MO  (3) 

Let  the  periodicity  assim9Jtion 

be  valid.  For  the  functions  gi  ,  in  view  of  the  variable  Uj ,  let  an  expansion 
in  Taylor  series  up  to  terms  of  the  second  order  be  possible,  and  let 
exan^jle,  be  continuous. 

Assume  that  a  solution  uo (t)  of  eq.(3)  periodic  with  P  is  known: 

4?=flK.0  +  M0.  u,(<  +  P)  =  11,(0.  (5) 

Then,  at  a  minor  modification  of  b(t),  we  will  investigate  the  solution  vectors 
u(t)  adjacent  touo(t)  of  the  system  of  equations 

47-  =  fl(«.0  +  M<)+/»fW  (6) 


17 


at  sufficiently  small  |B'  ;  here  again,  we  assume 

f{<+P)  =  f(0- 


(7) 

Z18 


Using 

UW  =  !!,(/)+ 1(0  (8) 

eq.(6)  will  be  transformed,  after  subtraction  of  eq.(5),  into 

^  =flKW  +  EW.<)-9(«.(0./)+/»f(0.  (9) 

Of  these,  the  solutions  with  small  lr(t)t  at  |3|  will  be  investigated. 

As  will  be  demonstrated  in  Section  2  where  a  transformation  of  eq.(9)  is 
to  be  made,  the  role  of  the  homogeneous  linear  system  of  equations*  is  taken 
over  by  the  system 


dt) 

dt 


(10) 


(11) 


with  the  matrix  (i,  k  =  1,  2,  n); 

I  \  du,  ‘  til,  tiu  I 

Obviously, 

«{/  +  /’)  910.  (1-2) 

The  homogeneous  system,  adjoint  to  eq.(lO),  will  then  read  as  follows  if 
the  transition  to  the  transposed  matrix  is  denoted  by  the  superscript  T: 


5*  --  -9l’'0j(/). 


(13) 


Theorem  1  (Resonance  case).  If  eq.(13)  has  a  solution  vector  r(t)  periodic 
with  P,  for  which 

(U) 


f*  P 


is  obtained,  each  solution  J’(t)  of  eq.(9)  -  independent  of  the  initial  values  - 
will  assume  values  with  increasing  t  whose  absolute  amounts  are  at  least  of  the 
order  of  magiitude  of  V'ls].  The  proof  is  given  in  Section  3. 


*  See  2ilso  another  report  by  the  same  author  (Bibl.2). 
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Then,  we  obtain  the  following  in  Section  4: 

Theorem  2  (Principal  case).  If  eq.(l3)  has  no  solution  vector  1  (t)  peri¬ 
odic  with  P,  thai  eq.(9)  will  have  solution  vectors  rCt)  for  each  sufficiently 
small  0  which,  for  all  values  of  t,  retain  the  order  of  ma^itude  |3|,  i.e., 
for  which 

IfWl^  Const.  |/?i  (15) 

is  valid  for  all  t,  for  example,  the  uniquely  existing  small  solution  periodic 
with  P: 

?(<+/')=?«• 

This  only  leaves  the  exceptional  case  to  be  defined  which  states  that 
eq.(l3)  will  have  solution  vectors  .^i(t),  ...,  ^r(t)  periodic  with  P  at  1  ^ 

^  r  ^  n  but  that,  for  all  these  periodic  solutions,  the  following  expression  is 
valid:  /19 

fiU*)Ht)dt  =  o  (e  =  i . r)  (16) 

I 

In  this  case,  it  is  iirpossible  to  iriake  general  statements  on  the  functions 
gi  (U,  t)  without  further  assxmptions. 

Section  2.  A  Transformation 

Putting 

u((,A)  =  u.(0  +  Aj(0,  (17) 

will  yield 

u(<,o)  =  u,(/),  u(/.i)  =  u,(0  +  jW  (18) 

When  applying  the  Taylor  series  with  respect  to  X  to  the  difference  on  the 

X 

right-hand  side  of  the  i^  conponent  equation  of  the  system  (9),  we  obtain 

?i(iio+  r.  <)  =  0  +  +  J  ')  .  (1  -  ;) 

or,  if  3(k(i)  are  the  conponents  of  the  vector  f(t). 
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(19) 


*(».  +  f.  0  -  &  («t.  /)  =  i;  — **  + 


In  addition  to  the  matrix  ?i(t)  from  eq.(ll),  we  next  introduce  the  ”taisor” 
T(u(t,  X),  t)  whose  n  coii?)onents  are  to  represent  the  matrices  (k,  =  1, 

n): 


»=  1, 


(X,{u(/.  !)./)> 

ItlutU)./) 

X.{u(U).0/ 


(20) 

(a) 


Then,  eq.(9)  can  be  written  in  the  following  form: 

=  «(/)!(<) +/»f(o  +  /VwT(u(/,A).OjW(i-X)ii  (22) 

Of  these,  at  sufficiently  small  e  and  3o  $  solutions  will  be  sought  with 


liWlSe  at  1^1 

It  should  be  noted  here  that  the  tensor 

T(u(/.O)./)=T(ii,(0./), 


(23) 


has  the  period  P: 


T(u,(<  +  P).<  +  P)=T(u,(0,<). 


(24) 


From  this  fact,  in  combination  with  the  first  relation  (23),  the  existence  /20 
of  a  finite  constant  M  follows  in  such  a  manner  that,  for  the  integral  in 
eq.(22),  the  following  is  valid; 

1/  f  (0  T  (u  (<,  A).  0 1 W  (1  -  A)  i a]  g  M  ■  (max  |  j (01)»  (  25  ) 

For  this,  it  merely  must  be  assumed  that  all  second  derivatives  on  the  right- 
hand  side  in  eq.(20),  within  an  interval  of  ti  ^  t  ^  ti  ♦  P,  are  restricted  for 
all  values  |u(t,  X)  -uo(t)|  ^  e  of  the  first  argument. 
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Section  3 •  The  Resonance  Case 


Let  the  adjoint  system  (13 )  have  a  solution  vector  ^(t)  periodic  with  P, 
for  which  eq.(l4)  is  vzilid.  In  contrast  to  the  argument  of  theoren  1,  we  are 
making  the  following  assumption; 

for  an  arbitrary  solution  r(t)  which  satisfies  the  first  stipulation  (23)  for 
all  t  with  a  finite  constant  N,  to  be  deteraiined  later.  Analogous  to  eq.(33) 
in  another  paper  (Bibl.2),  eqs.(22)  and  (13),  for  an  arbitrary  positive- whole  n 
with  C,  will  yield  on  the  basis  of  eq.(lA)i 


i^(0[f(<  +  »«P)  -  i(i)]  =  fitC  ■  fi  + 

t-^mP  I  /  \ 

+  /«^(T)/i'(T)T(u{T.A).T)j(T)(l-A)iArfT.  (27) 


Because  of  the  periodicity  of  ^^(t),  a  constraint  of  the  following  form  applies 
to  all  t: 

(28) 

If  eqs.(28),  (26),  and  (25)  are  used,  an  estimate  according  to  eq.(27)  will  be 

D-2S\P\*  +  DMN*\fi\mP:^m\P\\C\ 


or 

with 


2DN'^m\P\iB 

(29) 

b  =  \c\-dmn*p>o. 

(30) 

If  a  sufficiently  small  constant  is  substituted  for  N,  then  B  will  be  positive. 
For  each  0^0,  the  quantity  m  can  be  selected  so  large  that  eq.(29)  will  con¬ 
tain  a  contradiction.  This  proves  that  the  assumption  (26),  in  combination 
with  the  first  relation  (23),  is  inpossible  and  that,  therefore,  the  following 
theoron  applies; 

Theoran  3.  In  the  resonance  case,  four  constants  e.  M,  D,  N  exist,  from 
which  because  of  eq.(30)  a  positive  B  can  be  determined  in  such  a  manner  that 
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f 


each  solution  r(t)  of  eq.(9),  within  each  interval  of  the  length 
assumes  at  least  once  a  value  so  that 


is  justified. 


T.«mma.  If  3  is  restricted  by 

,  *ith  N*<  [see  eq.(30)3 

then,  in  each  interval  of  the  length  p  ,  the  following  estimate  will 


apply  for  at  least  one  value  of  t: 


Section  4.  The  Principal  Case 


In  the  following,  we  will  require  the  main  theorem  on  inplicit  functions 
as  an  auxiliary  theoran. 

Auxiliary  theoran  1»  Let  the  vector 

(»'iK . 

•'.K.  .«../»)  j 

. »,.P)/ 

in  a  certain  vicinity  of  the  quantities 

a,  =  o,  «,  =  o,  ....  «,  =  o.  fi  =  o 

possess  continuous  first  derivatives  to  av  and  let  this  same  vector  be  continu¬ 
ous  in  all  n  +  1  variables.  In  addition,  let 

»(o,  ...,0.0)  =  0;  (34) 

finally  be  the  ftmctional  determinant 


a^(o . 0. 0) 


+  0  (».*=1,2 . n). 


Thus,  for  each  sufficiently  small  a©  a  constraint  3o  will  exist  such  that,  for 
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each  0  with  13 1  ^  Bo  »  a  solution  vector  a  -  a (B )  of  the  systems  of  equations 

‘•(a.^)  =  o  (36) 

exists  uniquely,  for  which  ]  a|  ^  ao  •  For  the  proof,  see  for  exaJi?)le  another 
paper  (Bibl,3)« 

This  avixiliary  theorem  is  applied  as  follows:  Each  solution  vector  r(t) 
of  eq.(9)  is  characterized  by  its  initial  values: 


i.  •  9 


(37) 

(38) 


With  respect  to  t,  this  vector  has  the  period  P  if  and  only  if 

n(a./J)  P)  -  0)  =  o  (^^) 

This  corresponds  to  the  argument  (36)  of  the  auxiliary  theorem.  For  3=0,  /22 

the  null  vector  represents  the  solution  of  eq.(9);  consequently,  eq.(3^^)  is 
satisfied  since  then  also  a  is  the  null  vector.  The  assumptions  as  to  deriva¬ 
tion  and  continuity  of  the  auxiliary  theorem  are  ensured  by  the  following 
auxiliary  theoran,  in  accordance  with  known  theorems  as  to  the  dependence  of  the 
solution  vectors  of  eq.(9)  on  the  initial  values  a  and  on  the  parameter  3 
(Bibl.^). 

Aiuciliary  theorem  2.  Within  the  interval  0  s  t  ^  P,  for  a  given  e,  one 
positive  Sj  and  3i  each  can  be  determined  such  that  the  solutions  F(t)  of 
eq.(9)  satisfy  the  estimate  Ir(t)I  ^  e  in  the  entire  interval  if  only  |r(0)|  ^ 

^  61  and  |3|  ^3  1  are  selected. 

In  accordance  with  eq.(39),  it  is  necessary,  for  application  of  the  auxili¬ 
ary  theorem  1,  that  ei  ^  ao  as  well  as  3i  s  ,  a  differentiation  of  eq.(9)  to 
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indicatea  that  the  vectors 


-,*») 

(40) 

are  solutions  of  the  system  of  equations 

(a) 

with  the  matrix 

(42) 

Therefore,  the  vectors 

(43) 

are  solutions  of  eq.(lO)  since,  in  that  case,  eq.(iV2)  passes  over  into  eq.(ll). 
The  vectors  (43)  are  linearly  independent  since  their  determinant  for  t  =  0  is 
the  unit  determinant  [see  also  theoran  4  in  another  paper  (Bibl.2)].  Since,  in 
the  principal  case,  we  make  the  assunption  that  eq.(13)  and  thus  also  eq.(lO) 
[see  theorem  8  (Bibl.2)]  has  no  solution  vectors  periodic  with  P,  it  follows 
that  [see  theorem  10  ( Bibl.2 )] 

fltwr' 1=^0. 

A  brief  examination  of  eq.(39)  indicates  that  this  is  exactly  the  condition  (35) 
of  the  auxiliary  theorem.  Since  all  conditions  of  the  auxiliary  theorem  are 
proved  to  be  valid,  its  argument  -  in  our  case,  eq.(39)  -  can  be  considered  as 
also  proved. 

Therefore,  the  following  applies; 

Theorem  4.  If  eq.(13)  has  no  solution  vector  5<(t)  periodic  with  P,  con¬ 
straints  «i  and  01  will  exist  for  a  given  e  such  that  an  initial  vector  r(0) 
with  |j’(0)l  ^  exists  for  each  3  with  |3(  ^  3i ,  so  that  the  solution  p(t)  is 
periodic  with  this  initial  value  f(0)  and  satisfies  the  uniform  estimate 


Ir(t)l  e. 

This  represents  a  portion  of  the  theorem  2.  To  prove  this  theorem  /23 

completely,  we  will  need  still  another  theorem, 

Theoran  5»  Assume  that  a  finite  constant  E  exists  such  that,  at  suffi¬ 
ciently  small  1(^1,  the  following  will  be  valid  for  this  solution  f(t)  of  eq.(9) 
periodic  with  P: 

I?(01S£|/»I-  (A5) 

Proof,  If  :^i(t),  ,,,,  3n(t)  are  linearly  independent  solutions  of  eq,(13), 
the  following  expression  is  obtained  under  utilization  of  the  periodicity  of 
r(t)  analogous  to  eq,(38)  (Bibl,2)  for  v  =1,  n: 


UJ iJ (<)]  jr  (<)  -  P  I  il (t)  Hr)dr  + 

^  (a6) 

i- 1  (»■)  /  X  (t)  T  (u  (r.  A),  t)  j  (t)  (1  -  A)  rf A  it 

t  0 

[see  also  eq,(27)].  The  determinant  of  the  coefficient  matrix  on  the  left- 
hand  side,  according  to  the  assumption  of  the  principal  case,  differs  from  zero 
[see  eq,(AA)]«  Therefore,  the  linear  system  of  equations  (46)  can  be  solved 
for  the  conponents  Xi(t),  ,,,,  Xn(t)  of  xit)  on  the  left-hand  side,  using 
Cramer’s  solution  formula.  If,  at  fixed  t. 


Max  \x^{r)\  (v  =  1,  2, fi). 


(47) 


then  an  estimate  of  the  following  form  will  be  obtained  in  this  manner: 

f  L  (maxUWI)*  (48) 

with  two  finite  constants  K  and  L,  taking  eq,(25)  into  consideration. 

Since 

max  |j^(4|  S  V**  * 

it  follows  from  eq,(48)  even  more  so  that 

*  max  |j(/)| -g  |fi  A'l^j  f  JnL  •  (max|{(/)j)* 

max|f(/)|[i  -  V»«^  (niax|f(/)|)]gJnA|^|. 
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On  the  basis  of  theorem  4,  we  can  then  assume  that  the  expression  in  brackets 

,  1  1 

is  positive.  For  exanple,  for  expression  becomes  ^ 

From  this,  the  estimate  (45)  with  E  =  /nK  is  imnediately  obtained. 

Instead  of  using  Cramer’s  rule  for  the  solution,  it  is  naturally  also 

possible  to  solve  eq.(46)  for  ?(t)  by  left-hand  multiplication  with  the  matrix 

inverse  to  [^v(t  +  P)  -  ?»v(t)J  and  then  to  make  the  estimate  by  means  of 

Schwarz*  inequality  lr(t)|.  This  will  also  yield  eq.(45). 
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PART  III 


/27 


STUDY  OF  THE  RESONANCE  CASE  IN  SYSTEMS  OF  LINEAR 
ORDINARY  DIFFERENTIAL  EQUATIONS 


Section  1.  Introduction 

Previously  (Bibl.l)  we  investigated  the  linear  inhomogeneous  differential 
eqxiation  systan 

with  a  n-row  matrix  S((t)  which,  for  sinplicity,  was  assumed  as  being  continuous 
and  having  a  continuous  n-component  vector  f(t)  of  the  same  period  P 

/^(t  ♦  P)  -  ^(t)  ,  t  ♦  p)  -  ^(t)  (2) 

The  corresponding  homogeneous  system  then  will  be 

(3) 

and  the  adjoint  system 

g  ■  -  ,  (^) 

where  the  superscript  T  denotes  the  transition  to  the  transposed  matrix. 

In  all,  three  cases  were  differentiated;  The  principal  case  is  present  if 

eq.(4)  ’has  no  solution  periodic  with  P;  the  resonance  case  is  present  if  eq.(A) 

has  at  least  one  solution  vector  5?(t)  periodic  with  P,  for  which 

p 

J  -f  (r)dr- c  t  o  (5) 

is  valid;  the  exceptional  case  is  present  if  eq.(/4,)  does  have  periodic  solu¬ 
tions  ^i(t),  ^3(t),  ...,  Sp(t)  (1  ^  p  "S  n)  periodic  with  P  but  if  the  following 
is  valid  for  all  these  ^u(t): 

p 

/  ^  -  O  for^.  1,2 . 

Whereas,  in  the  principal  case  as  well  as  in  the  exceptional  case,  solutions 
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of  eq.(l)  periodic  vdth  P,  i.e.,  ronaining  limited  for  all  values  of  t,  are  in 
existence,  the  values  of  all  solution  vectors  f(t)  of  eq.(l)  tend  toward  in¬ 
finite  with  increasing  t  in  the  resonance  case,  independent  of  the  initial 
values.  The  individual  steps  in  this  process  will  be  further  investigated  in 
the  present  paper. 

We  use  the  following  notations  [s 

^9  ■  f 

for  denoting  a  fundamental  system  of  solutions  of  eqs.(3)  and  (4),  respectively 
[see  eq.(17)  in  a  previous  report  (Bibl.l)].  In  the  method  of  variation  of  the 
constants*,  using  the  abbreviation 

®  (  y  •i.S, . . . ,n)  (8) 

the  following  argument  is  constructed  for  the  solution  of  eq.(l): 


r(t)  ■  p  -  ^(t).  ^(t) 


In  the  present  report,  an  argument  similar  to  eq.(8),  namely, 

x.  ^  ^  “ 

[see  eqs.(77)  and  (78)]  is  discussed,  where  the  p^^t'’^(t)  are  composed  in  a 
suitable  manner  from  the  pv(t).  In  the  resonance  case,  statements  can  be  made 
on  the  behavior  of  these  vectors  p^jj,'*^(t)  for  t  increasing  without  bounds.  Here, 


the  matrix 


^  ♦  P), 


which  had  been  discussed  in  Section  2  plays  a  decisive  role.  By  a  siiitable 
transformation,  the  system  (l)  can  be  brought  to  a  "normal  form"  (Section  3) 
with  constant  coefficients,  by  means  of  which  the  study  of  the  vectors  (t) 


V.dth  respect  to  the  method  of  variation  of  the  constants,  see  footnote  on 
P*  58. 
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can  be  made  much  more  concise  (Section  4).  Making  use  of  this  method,  four 
numerical  examples  will  be  given  at  the  end  of  this  paper  (Section  5)» 

A  brief  report,  to  be  published  soon,  will  apply  the  above  considerations 
to  the  case  of  an  ordinary  differential  equation  of  the  n^^  order  with  periodic 
coefficients.  In  that  paper,  additional  examples  will  be  given. 


Section  2.  The  Matrix  SB 
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Theorem  1;  The  matrix  defined  in  eq.(lO),  is  a  constant  regular  matrix, 
i.e.,  a  matrix  independent  of  t,  vmder  the  first  assunption  [eq.(2)J. 

Proof;  In  view  of  eq.(7),  eq.(3)  can  be  written  in  the  form  of  » 


Then,  a  fundamental  solution  matrix  of  eq.(4)  will  be  as  follows  [see  (Bibl.l), 
eq.(17)]; 


which  means  that 


J2* 

-  -/t'ct)  ^  . 


(12) 


(13) 


is  valid.  According  to  eqs.(lO),  (12)^  (13)»  and  (11)  we  then  obtain 


^  *  P)J  +  p)  ^’(t  ♦  P)  - 

-  ♦  P)  ♦  P)^(t  .  P)  -  IT. 

when  taking  the  first  relation  (2)  into  consideration.  Consequently,  ^  is  a 
constant  matrix.  That  its  determinant  differs  from  zero  follows  from  the  non¬ 
vanishing  of  the  two  determinants  of  the  matrices  on  the  right-hand  side  of 
eq.(lO). 

Theorem  2;  If  we  pass  from  a  fundamental  solution  matrix  5)(t)  of  eq.(ll), 
on  multiplying  on  the  right-hand  side  by  a  regular  constant  matrix  5,  to  a  new 
fundamental  system 
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^‘(t) .  (14) 

[see  (Bibl.l),  Theorem  5],  the  following  relation  will  be  valid  for  the  matrix 
^  which  had  been  formed  in  accordance  with  eq.(lO); 

(15) 

Proof;  From  eq.(lO)  and  (14)»  we  obtain 

-j^o  \  t  )^®(t  ♦  P)  -  j:'^  2£‘~^(t)^(t  * 

This  means  that,  by  a  siiitable  selection  of  the  fundamental  system  sKt),  the 
matrix  ^  can  be  transformed  into  any  normal  form  which  can  be  produced  by  a 
similarity  transformation  of  the  type  of  eq.(15)» 

Theorem  3 ;  Let  the  constant  matrix  Sf  be  given  in  the  Jordan  normal  form  /30 
[see  (Bibl.2),  Sect.l9«l] 


where  one  elementary  conponent  of  the  order  my  ^  1  in  the  main  diagonal  has 
the  eigenvalue  Xy  and  contains  ones  in  the  next  higher  diagonal  (zeros  are  not 
entered).  It  is  obvious  that 

“  "  “i  ♦  *2  ♦  •••  ♦  *8  •  (17) 

■» 

Then,  it  is  possible  to  form  a  matrix  Hq  with  any  real  nramber  q,  so  that 

*  .  (18) 

is  valid  and  that  fi,  has  a  normal  form  of 
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r 


h' . 

4 


th 


\ 


where 


and 


i  In  I,  with  -  <  7  (  •‘y  )  ^  — 


(•i 


fi  .(-1)*'“"’  - - — p-  f"'  '*• 


Proof:  Let  us  form  [see  also,  for  exan5)le,  (Bibl.3),  pp«333/4] 

In  ^  -  ln(  +JJ;  )  , 

where  ®v  is  the  unit  matrix  of  rank  my  while  the  my-row  matrix 


N 


has  the  rank  my  -  1.  By  expansion  in  series,  we  obtain* 


In 


if  -  ln(  >^y€y)  ♦  ln(  ^  “ 

■  (In  Ay)  ify  *  1b(  "Ay""^  ^ 


*  The  formula 


In(Ayi^)  .  (InAy).^ 


used  here  is  equivalent  to 


Vy  C 


.  *  Ay  * 


(19) 


(20) 

(21) 

Z21 

(22) 


(23) 


as  is  readily  verified  by  means  of  a  power  series.  On  the  basis  of  this  fact, 
we  can  also  write  in  abbreviated  foim; 
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m 

V 

[see  also,  for  exan^jle  (Bibl.4),  p.82].  For  this,  it  is  suggested  to  use  the  /32 
principal  value  of  the  natural  logarithm  In,  i.e., 

-T<  :7(ii»Ay  )*<  T  /  (24) 

which  agrees  with  eq.(20).  Taking  eqs.(2l)  and  (19)  into  consideration,  we 
will  obtain 

from  which,  by  means  of  expanding  the  e^onential  function  in  a  series  [with 
respect  to  the  convergence,  see  for  eacanple  (Bibl.4),  p.ll9], 

*  -  ’  (26) 

is  obtained.  From  this,  eq.(18)  with  can  then  be  taken  from  eq.(19)  since, 
in  forming  the  powers  of  ,  the  elementary  components  Jin  v  do  not  exert  a  mutual 
influence. 

an  additional  collineatoiy  transfoitnation,  the  matrix  in  eq.(19)  can 
be  brought  to  the  Jordan  normal  form.  If  5  represents  the  matrix  of  this  simi¬ 
larity  transformation,  i.e., 

then  SP  is  transformed  into  5^5,  as  readily  demonstrated  by  expanding  the  ex¬ 
ponential  function  in  a  series.  We  can  now  write 
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(27) 


<k 


4 


•i  tk 


1 


This  will  yield  the  following: 

Theory  Li  The  matrix  ^  can  be  brought  to  the  following  form  by  a  sind-  /33 
larity  transformation: 


*0 


(28) 


with  from  eq.(27)« 

L«nma:  If,  specifically,  all  elanentary  components  of  the  matrix  S|l  have 
the  order  1,  also  fiq  will  be  a  diagonal  matrix. 

Definition:  Below,  we  will  r^slace  the  arbitrary  number  q  by  the  period  P 
of  the  coefficients  of  eq.(l);  in  this  case,  we  can  write  in  abbreviated  form: 


(29) 


Theoran  5:  By  means  of  the  constant  matrix  ft  [see  eq.(29)  and  (27)],  the 


fundamental  system  n(t)  of  eq.(3)  can  be  written  in  the  form 

I^(t)  .  »i»k  ^(t)  .  (30) 

where  4(t)  has  the  period  P.  AnaJLogously,  the  fundamental  system  3(t)  of  eq.(A) 
can  be  brought  to  the  form 

^(t)  .2^(t)(e  (31) 

where 

2^(t)  -  (  ^  (32) 

also  has  the  period  P. 

Proof:  It  is  to  be  demonstrated  that,  in  the  argument  (30), 

(p(t)  -  j^(t)  •  (33) 


3U 


has  the  period  P*  Making  tise  of  eq8*(10),  {29),  and  (30),  we  obtain 

I2k 

^(t  4  p)  -  l^(t  ♦  p)«-^*’.  g-rft  .yt)^ 

-At  r 

•  e  -  ^(t). 


Th®n,  the  relations  (31)  and  (32)  are  directly  obtained  from  eqs.(12)  and  (30). 

The  conventional  way  of  constructing  the  constant  matrix  described  in 
developing  the  theorem  3»  usually  is  quite  time-consuming.  That  this  procedure 
can  be  less  cumbersome  in  certain  special  cases  is  indicated  by  the  following 
theoran*. 

Theorem  6;  If  the  matrix  a(t)  with  its  integral  from  0  to  t,  periodic 
with  P,  is  transposable,  i.e.,  if  the  relation 

t  t 

y  A  (r  )dr-  / /t  (r)dr  • /l(t)  (34) 

«  0 

exists  identically  in  t,  thai  the  constant  matrix 


will  yield 


(35) 

(36) 


To  this  belongs  the  fundamental  solution  matrix 

J* ^tndr  , / 

l^(t)  -  c  •  f  C 

with  the  following  matrix,  periodic  with  P, 

**  /  ( lk<rt- 

§  (t)  -  / 


(37) 


(38) 


Proof:  Below,  the  two  fomulas  which  can  be  proved  by  expanding  the  ex¬ 
ponential  fiHiction  in  a  series  are  used  [see  also,  for  example  (Bibl.4),  /35 


*  I  wish  to  express  ny  thanks  to  Dr.H.Eltermann  for  developing  the  idea  of 
theorems  6,  7»  and  8. 
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pp.121  and  138]: 


.  H 


S-itt  •  '  i-«»  ,  / 

(e  /.  C  PWMfi-ut  e  ,  if 


First,  it  can  be  demonstrated  on  the  basis  of  eq.(4.0)  that  r^(t),  according  to 
eq.(37),  is  a  fundamental  solution  matrix  of  eq.(ll).  [But  note:  ?1(0)  =®.J 
The  constant  matrix  (see  theorem  l)  can  then  be  calculated  from  eq.(lO)  vdth 
t  =  0,  yielding  ? 

if  -  l^(P)  -  e  ,  (^1- 

from  which  eq.(36)  is  obtained  with  eq«(35)« 


It  only  remains  to  prove  that  the  matrix  $**( t  }  defined  by  the  argument 
(37)  in  accordance  with  eq.(30),  can  be  written  in  the  form  of  eq.(38).  Because 
of  eq«(35),  ^^^(t)  has  the  period  P,  since 


MT 

j  iSl(r)-  A  )4t  Ct 


The  relations  (37)  and  (38)  can  be  reduced  to  the  identity 


j  d  *)Vr 


if  it  can  be  demonstrated  [see  eq.(39)]  that  the  matrix  ft  is  commutative  with 

t 

the  matrix  f  Sl(T)dT,  i,e«,  if 


J  Cl <xMr  ■  J  =  J 


is  valid.  Since  this  eqioation  is  directly  understandable  for  t  =  0,  eq.(44) 
will  be  obtained  by  differentiation  to  t,  based  on  the  commutability  relation 


J  ■  Cl (n  «■  h 
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Taking  the  periodic! t7  of  8{(t)  into  consideration,  we  can  calculate 

^  I 

rT>Vr  jm 

-  Dt<  TjVr 

and,  taking  eq.(3i4.)  into  consideration. 


Z26 


(U6) 


/  {•tiJx  ^  <  t -tp)  -  J  T  ■  ^(t)  * 
Q[iXi^r  ~  J  z'^  ^llti  -  j  ^<xUt  ^lt>  , 


which  proves  eqs.(i4.4)  and  thus  also  eq.(43)« 

By  a  similarity  transformation,  the  matrix  (35)  can  now  be  brought  to  the 
Jordan  normal  form,  indicated  in  eq.(27),  in  such  a  manner  that  the  sequence 
of  Oj  ,  QTa ,  ...,01  is  the  same.  This  will  yield  a  matrix  f{  ,  for  which 


e  .  e 


(47) 


with  ft  from  eq.(29)  is  valid.  The  eigenvalues  of  the  matrices  (47)  thus  are 


mutually  equal.  Accordingly,  the  eigenvalues  ory  of  H  and  H  can  differ  at  most 
2^1 


by  nv  * 


with  integral  ny .  Thus,  we  have 


ith  'Co  »  It  ■ 


(48) 


[see  also  eq. (27)],  if  ®v  denotes  the  my-row  unit  matrix  [see  eq.(17)].  This 
will  lead  to  the  following  theorem; 

Theorem  7:  By  a  similarity  transformation,  fi^"‘  from  eq.(35)  can  be  /37 
transformed  into  a  matrix  ^  which  is  correlated  with  eq.(48)  over  from 
eq.(29),  so  that  Si  can  be  calculated. 

Theorem  8;  The  assun^stion  (34)  for  "the  theorans  6  and  7  can  be  replaced 
by  the  stronger  assumption 
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(X  ^  •  IX(r)  /^(t) 


(49) 


at  arbitrary  t  and  t. 

Proof;  It  is  to  be  dmonstrated  now  that  eq.(34)  (with  arbitrary  lower 
limit)  follows  from  eq.(49)«  We  calculate: 

Plat  u  IX>Jr  (50) 

4  t  4  4 

Section  3»  Transformation  of  the  Differential  Equation  System  /38 

to  a  Normal  Form 


Theorem  9;  Ey  means  of  the  transformation 

Z(t)  -  $(t)  •  v(t)  (51) 

with  ^(t)  from  eq.(33)»  the  system  of  differential  equations  (1)  is  transformed 
into  the  following  system  [Floquet’s  theoran,  see  for  exanple  (Bibl.3), 
Chapt.III,  Sect. 5,  and  specifically  p.75]s 


(52) 


with  [see  eq.(32)] 


T 

-^(t)  -  l^(t)  (t) 


(53) 


which  has  the  constant  matrix  fi  according  to  eqs.(29)  and  (27)  and  also  has  a 
6(t)  periodic  with  P.  We  will  call  eq.(52)  the  normal  form  of  our  system  of 
equations  (l). 

Proof;  In  view  of  eq.(32),  a  substitution  of  eq.(5l)  into  eq.(l)  and 
consideration  of  eq.(53)  will  yield 

19' •  4^  I  -  §')  t  ^  ■ 


It  then  ranains  to  demonstrate  that  the  matrix  to  be  applied  to  b  is  equal 

to  «,  i.e.,  , 

f  ~  .  (5/J 

In  fact,  repeatedly  taking  eqs.(33)  and  (11)  into  consideration,  we  obtain 
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-  Dlige^'i^  e^‘  e'‘*‘^)-a€. 


Theorgn  10;  The  homogoieous  system 

MU*  •<&  (5< 

conjugate  to  eq.(52)  has  no  other  solutions  periodic  vdth  P  than  solutions  of 
the  form 

-*4»  .  (5' 

Naturally,  an  analogous  statement  applies  also  [see  eq.(12)]  to  the  system 


^  ^  ^  -tL  • 


(58) 


adjoint  to  eq«(56). 

Proof;  The  system  of  equations  (56)  is  resolved  into  the  mutually  inde¬ 


pendent  systans 


Here,  we  then  have 


Such  an  individual  syston  (59)  possesses  the  following  my -dimensional  solution 


vectors 


-  e 


with  an  arbitrary  constant  vector  Cv  Because  of  [see  eqs,(27),  (23),  and 


footnote  on  p.32] 


^  9ith 

y  ^y 


eq.(6l)  will  be  transformed  into 


,  \  (62) 

/I  »  f  £-l_\ 


'‘r*  X*  It  ...  *  i 

■HffyiXi)  -  e  -e  ^  =  «  I- A'" 

1...  iii/ 
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Because  of  eq.(20)  (with  q  -  P),  the  vector  (63)  is  periodic  with  P  if  and  /hO 
only  if 

o  and  <;-  /  o  \ 


where  Cv  has  been  normed.  The  only  nortned  solution  vector  of  eq.(59),  periodic 
with  P,  will  thus  be 

-  (:)•  »/  -  o.  (65) 

y  o  j 

The  corresponding  solution  vector  tD(t)  of  eq.(56),  which  is  periodic  with  P  be¬ 
cause  it  is  constant,  will  have,  except  for  a  1  at  the  point 


(  y )  -  ♦ 


♦  ■  y-1  ♦  ^ 


with  Ofy  =  0,  only  zeros  as  con^jonents.  The  n-conponrant  solution  vector  of 
eq.(56),  conjugate  to  eq.(65),  will  be  denoted  byh)(v)  •  This  means  that  in 
eq#(60),  all  h>u  “  0  must  be  used  at  M-  /  v,  so  that  fcv  will  have  the  value 


[eq.(65)]: 


:  ,  i/  o 

r-' 


This  will  yield  the  auxiliary  proposition: 

Lemma:  The  constant  solution  vectors  of  eq.(56),  i.e.,  the  vectors  peri¬ 
odic  with  P,  are  linearly  conposed  of  the  vectors  (6?)  withttv  from  eq.(65). 
Analogously,  the  normed  constant  solution  vectors  for  eq.(58)  are 


.  i/  -  O  , 


with  the  mv -component  vector 


(69) 


where  Ujvj  ,  in  addition  to  zeros,  has  only  one  component  1  at  the  point 


*  *2  *  * 


with  •<i>  ■  o  . 


Taking  eqs.(30)  and  (31)  into  consideration,  we  obtain  directly: 

Theoran  11:  The  solution  vectors  h(t)  of  eq.(3)  and  3(t)  of  eq.(4),  peri¬ 
odic  with  P,  can  be  written  in  the  form  of 

-  i/")*  (71) 

and,  respectively, 

Ut)  =  2^(t)  (t)  -  y(t) 

rvj-*  M 

where  (v)  and  [v],  respectively,  have  the  meaning  of  eq.(66)  and  (70). 

Theorem  12:  For  the  systans  of  differential  equations  (l)  and  (52),  we 
always  have  simultaneously  the  principal  case  or  the  resonance  case  or  the  ex¬ 
ceptional  case. 

Proof:  According  to  theoran  11,  the  solutiofis  of  eqs.(A)  and  (58),  periodic 
with  P,  are  at  a  one-to-one  correspondence.  In  addition  [see  eqs.(5)»  (53\ 
and  (71)],  the  following  applies: 

p  ^ 

«  •  y  (72) 

where  the  conventional  notation  was  used.  Now,  the  correctness  of  theorem  12 
can  be  read  directly  from  the  second  paragraph  of  Section  1. 

The  matrix  ^(t),  defined  in  eq.(33),  is  broken  up  into  the  sxim  of  s  n-row 


square  matrices 


#  ■  (V)*  -  t  f*'  , 


where  the  matrix  ,  in  the  columns  with  the  numbers 

(y),  (  y)  ♦  1,  ...,£:yj  contains  the  vectors 
but  otherwise  only  zeros*.  The  system  of  differential  equations  (52)  decomposes 


->  Here,  the  notation  of  eqs.(66)  and  (70)  has  been  used  for  the  first  time 
without  the  restriction  Q^v  ^  0. 


into  the  independent  subsystems 

^  y-  1,2 . s,  (74) 

where  vie  always  have  the  dimension  mv ;  here,  bv(t)  contains  only  the  conponents 
(v),  [v]  of  b(t).  Then,  we  have 


I 

V2(t) 


<>  (t) 


(75) 


if  the  n-con^ionent  vector 


(»> 


tC 


/?  ' 
6 

KV'w 

o 

U/ 


,  y  -  1,2, ...8 


(76) 


is  defined* 

It  would  seem  logical  to  introduce  the  vectors 


c'Jt)  = 


(77) 


C<J  _ 


in  which,  however,  all  n  con?)onents  may  differ  from  zero,  for  which  reason  the 
denotation  was  supplemented  by  an  asterisk.  This  will  yield: 

Theor«n  13:  Each  solution  j(t)  of  eq.(l)  can  be  written  in  the  form 


with  f^'’^(t)  from  eq.(77). 

Lemma:  If,  specifically,  SP  has  a  diagonal  form,  eq.(78)  with  (v)  =  [v]  =  v 
specializes  to 


■C 


L 

y=i 


/‘-T 

U2 


(79) 


Resolving,  in  accordance  with  eq.(73)> 

/LL 

♦  ...  ♦  ^ 

(80) 

as  well  as 

iH  - 

(81) 

the  following  theorem  is  obtained: 

Theoran  iZi:  The  vectors  r^'’^(t),  defined  in  eq.(77)f  satisfy  the  s  systans 
of  n  differential  equations  each 

y-  (82) 

with 


^(t)  (83) 

\ander  use  of  the  above  ssrmbolism. 

Proof:  By  differentiation  of  eq.(77)  and  with  consideration  of  eqs.(7/4.) 
and  (76),  we  obtain 


It  is  readily  verified  that 

(84) 

in  which  case  the  form  of  the  matrix  fi  must  be  taken  into  consideration. 

This  will  further  yield 

from  which,  according  to  eq.(77),  the  correctness  of  eq.(82)  with  the  first  /45 
relation  (83)  is  obtained.  The  second  relation  (83)  can  be  perceived  as  fol- 

T 

lows:  The  zero  columns  of  are  transferred  into  zero  rows  of  .  Conse¬ 

quently,  it  follows  from  eq.(53)  that 


-f'-. 


(85) 
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Section  4.  Discussion  of  the  Syst^  of  Differential  Equations 
in  the  Normal  Form 


Theorem  15;  Let  the  matrix  ^  have  a  diagonal  form  and  let  be 

the  independent  solution  vectors  of  eq.(4)  with  the  period  P  and,  finally,  let 


Ceq.(72)] 


I  ir  •  I  m 


fo  for 


0  for  J 


Then,  for  1  ^  M-  ^  (partial  resonance  case)'*^  each  component  of  f^.j.^^(t)  [see 
eq.(79)]  for  which  the  corresponding  component  Hu(t)  does  not  vanish  identi¬ 
cally  assumes,  with  increasing  t,  values  of  the  order  of  magnitude  of  t;  for 
M-  =  Cl  +  1,  p  (exceptional  subcase)^,  all  these  (t)  are  periodic  with  P; 

•if 

for  M'  =  D  •*•  1,  ...,  n  (partial  main  case)  exactly  one  r^j.^'(t)  periodic  with  P 
exists. 

Proof;  Since  ^(t)  is  periodic  with  P  (theorem  5),  the  analogous  theorems 
need  be  proved  only  for  the  (t).  These  have  only  one  conponent  ’-ii(t)  dif¬ 
fering  from  zero  which,  according  to  eqs.(52)  and  (74),  satisfies  one  equation 


each  of  the  form 


Vytt  *  ^  *1,2,  ...,11) 


The  general  solution  reads 


yCt)  =  J  C  ^  C  e 


with  an  arbitrary  constant  Cu .  In  the  case  /  0  (u  =  p  +1,  ...,  n),  accord¬ 
ing  to  eq.(88)  and  from  the  condition 

(t  *  p)  .  ^  (t) 

the  constant  q^  can  be  uniquely  defined  as 


TMs  case,  however,  must  not  necessarily  occur. 


J- 

‘r-  f-.*  J‘ 


1 


- 


M 

(89) 


In  the  case  =  0,  according  to  eq.(88),  we  have 

t 

%2 . f.  (90) 

For  M-  -  ^  *  1,  ...,  P,  this  is  periodic  with  P  at  any  selection  of  the  con¬ 
stants  qi .  At  M.  =  1,  •••f  CJ,  eq.(90)  can  be  written  in  the  following  form  [see 


eqs.(86)  and  (88)]; 


v^(t)  •  J-a  t  -e 
-p  ^ 


Since  the  integral  is  periodic  with  P,  the  linear  increase  of  Mii(t)  with  t 
follows  from  eq.(91)  and  the  corresponding  statanent  fory^'^^(t)  follows  from 
eq.(79). 

This  theoran  is  a  special  case  of  the  following; 

Theorau  16;  Let  ^  not  necessarily  be  a  diagonal  matrix  but  let  it  be 

ftp 

given  in  the  form  of  ^  =  e  [see  eqs.(29)  and  (27)3,  with  the  elementary  com- 
ft^P 

ponents  ^  =  e  ,  where  the  quantities  ftv  are  the  elementary  components  of  the 
Jordan  normal  form  of  the  matrix  ft  [see  eq.(27)  with  q  =  P].  If  0°^**  f  1 
(principal  subcase),  then  eq.(7A)  will  have  a  uniquely  defined  solution  vector 
’*(v)(t)  periodic  with  P,  so  that  [see  eq.(76)]  the  vector  (t),  defined  in 
eq.(77),  has  the  period  P.  If  e*'’^  =  1  and  [see  eqs.(70),  (71),  and  (72)  as 


well  as  theorem  12] 


r  ^ 


(exceptional  subcase),  theft  a  one-parameter  family  of  solution  vectors  *\(t)  of 
e<3»(7/4),  periodic  with  P,  will  exist  from  which  a  corresponding  family  of 
vectors  will  resrrlt.  Finally,  if  the  following  relation  exists  in 


45 


addition  to  =  1 


r  r 

[  r  -  fbco 

»  A 


(resonance  subcase),  then  each  solution  vector  Dv(t)  of  eq.(74)  and  thus  also 

each  r^'’^(t),  with  unbo\indedly  increasing  t,  will  take  values  of  the  order  of 
■  V  of 

magnitude  t  ,  provided  that  my  is  the  order  of  and  fiy .  (A  supplanent  to 
this  is  contained  in  theoran  17. ) 

Proof;  In  the  case  Oy  0  (principal  subcase),  let  us  successively  solve 
the  system  of  equations  (74),  starting  with  the  last  equation. 


J  l"  ,  (90 

m  J  ^  ^  “M-b 

The  condition  Xi(t  +  P)  =  •vLi(t),  analogous  to  eq.(89),  will  successively  yield 

r  “*4* 

I  f-rr  , 

Thus,  the  vector  Dy(t)  periodic  with  P  and  thus  also  (t)  is  viniquely  defined 

w 

in  this  case  as  well  .  If  ofy  =  0,  then 


^C¥l  *  CrJ  > 


il  ]°°" 


‘to'”''’' 


*  If  gfl(ory)  /  0,  the  integrals  in  eq.(94)  with  t©  =  sign  5R(Q'y  )  *  ®  instead  of  0 
as  lower  limit,  will  automatically  have  the  period  P,  as  can  be  readily  veri¬ 
fied.  In  that  case,  the  last  summands  in  eq.(94)  must  be  omitted  [qj  =  0  for 
V-  =  l^h  [v]  -  1,  ...,  (v)]. 


[see  eq.(86)]  is  used  for  deconposing  the  vector  in  eq.(74) 


(97) 


xhere 


where,  because  of  theorem  10,  it  still  can  be  stipulated,  for  example,  that  the 
first  component  ofQv(t)  satisfies  the  condition 


0 

(100) 

to  the  two  systems 

(101) 

u>*  (t)  -  ♦  ZlfCt). 

(102) 

Solution  of  eq,(102)  in  a  manner  analogous  to  that  used  for  eq.(94),  will  yield 


t 

+  cjf.  /•<.'  W-v... 


(103) 


The  condition  D  (t  +  P)  =  D  (t),  because  of  the  fact  that  1ii(t)  =  tti(''’)  is  valid 
for  M-  =  [v]  -  1,  ...,  (v),  will  yield 


?  ? 

/  *0  and  *0- 


(104) 


Because  of  eq.(98),  the  first  condition  is  automatically  satisfied.  Conversely, 
the  remaining  conditions  (IC4)  successively  yield  vinique  values  for  the  constants 
c^cvj  ,  c*tv3-i,  .«•»  cfv;  +  i>  whereas  c*v>  remains  conpletely  arbitrary.  This  / 50 
means  that  eq.(102)  has  a  one-parameter  family  of  solutions  Ov(t)  periodic  with 


P.  If  now,  in  addition  to  ofy  ^  0»  also  eq.(92)  is  satisfied  (exceptional  sub¬ 
case),  then  Ov  of  eq.(96)  is  the  null  vector  eind  eq.(102)  has  the  zero  vector 
as  the  only  solution  vector  periodic  with  P  which  satisfies  eq.(lOO),  This 
means  that  the  above  calculated  one-parameter  family  of  solutions  Ov(t),  peri¬ 
odic  with  P  and  having  cfv)  as  parameter,  furnishes  the  only  solutions  Ov(t)  = 
=  Ov(t)  of  eq*i7U)t  periodic  with  P.  From  this,  the  quantities  r^jj.'’^(t),  peri- 
odic  with  P,  can  be  calculated  according  to  eq.(77).  Conversely,  if  eq.(93) 
coexists  with  Oy  =  0,  the  system  (101) 


(105) 


is  solved  in  the  form  of 


“  Polynomial  of  tke  degree  >1^  •••  *y'*’  ^ 

-  M  ♦  1  -e-  (  #-.  (y )  ,  (y)  +  1 ,  .  .  .  ,  M  ) 


(106) 


T  CO 


with  the  highest  coefficient  .  Consequently,  specifically  q(y)(t) 

will  be  a  polynomial  of  the  degree  my*  Since,  because  of  eq.(98),  the  terms 
VCT(t)  in  eq.(103)  increase  at  most  like  t  ,  meaning  that  v^yj(t)  at  most 

■  y“  1  . 

will  increase  like  t  ,  each  solution  t)y(t^  with  increasing  t,  will  take 
values  of  the  order  of  t  in  accordance  with  eq*(99)*  The  same  statement  then 
also  applies  to  the  vectors  (t)  to  be  calculated  from  eq.(77)  since,  accord¬ 
ing  to  eq.(30),  cp^yj(t)  cannot  be  the  null  vector. 

Without  further  proof,  these  same  considerations  show  directly: 

Lemma:  The  components  Vj(t)  (o  =  (v),  (v)  +1,  *,.,  [vj),  in  the  resonance 

C  V3  +  1-c 

subcase,  take  values  of  the  order  of  t  with  increasing  t.  The  same 

statement  applies  to  the  vectorial  conponent  svunmands  'Pa  (t)v^(t),  occurring 
in  X’^.^^(t)  in  accordance  with  eq.(77),  for  each  scalar  component  for  which  the 
corresponding  con^)onent  of  9_(t)  differs  from  zero  (see  theorem  (17).  / 53 
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Theoram  17:  In  the  resonance  subcase,  each  con?>onent  (o  =  1,  2, 

•  n)  of  from  eq.(77)  takes  arbitrarily  large  values  with  increasing  t, 

unless,  for  exaaple  for  o  =  p,  the  components  c^Q^u(M•  =  (v),  (v)  ♦  1,  [v]) 

of  the  row  vector  of  ^  [see  eq.(33)]  all  vanish;  in  this  case,  we  have 
x^'’p(t)  2  0.  More  accurately,  it  can  be  stated:  If  we  do  not  have  xl^^.'^^Ct)  =  0 
(see  above),  then  x*''p(t)  (p  =1,  2,  ...,  n),  with  increasing  t,  will  take 
values  of  the  order  of  t*'*  ^  provided  that  in  the  p^^  row  of  [see  eq.(73)], 
the  qiwntity  9o^v)'*'3^  first  conponent  differing  from  zero. 

Proof:  According  to  eq.(77),  (99),  and  (106),  we  can  write 

if  Q^'*l[t)  or  D^'^^t)  represent  the  n-conponait  vectors  supplemented  by  zeros 
from  q^(t)  and  Dy(t).  Rather  than  by  eq.(lOO),  the  resolution  of  eq.(99)  is  now 
restricted  by  the  stipvilation  that  Dv(t)  is  to  be  periodic  with  P  and  normed. 


In  addition,  we  have 


<(t).  V 


(107) 


[see  eq.(33)].  In  this  case,  •••»  ^ivj(t)  are  linearly  independent 

vectors  with  the  period  P  while  c^(t),  according  to  eq.(106),  is  a  polynomial 
in  t  of  the  degree  [note:  o  =  (v)  +  A]: 

Cia+  1  (108) 

Thus,  the  p^^  conponent  of  eq.(107),  taking  eq.(105)  into  consideration,  will 

/52 

become  W  - 


Since  the  first  stim  in  eq.(l09)  is  periodic  with  P,  it  will  remain  finite  for 
all  t.  If  9pg[v)  +  A  is  the  first  nonvanishing  coefficient  in  the  second  sum, 

/  y  \  .  *  y"  ^ 

th«i  x^  (t)  will  take  values  of  the  order  of  magnitude  t  with  increasing  t. 
We  still  note  the  following  (slightly  weaker)  alternative: 


Alternative ;  Either  x^^^(t)  assumes  arbitrarily  large  Veilues  with  increas¬ 
ing  t  or  else  (t)  =  0  is  valid. 

Note;  If,  corresponding  to  eq.(73) » 

^  (110) 

is  defined,  then  from  the  vanishing  of  all  elanents  of  the  row  of  be¬ 

cause  of  eqs.(30)  or  (33 )»  the  same  statanait  follows  for  and  vice  versa. 


Section  5.  Examples 


Example  1;  (Re  theorem  15,  resonance  case  and  exceptional  case,  and  theo- 
ran  8). 

Let  us  consider  the  system  of  differential  equations 

cost  ♦  Sint  ♦  f^(t) 

x^  -  -  x^  Sint  Xg  cost  ♦  fgCt)  (111) 

vmder  the  following  two  assun?5tions: 


and 


f^(t)  •  *co8  (l  -  cost) 

^2^^)  •  -  sin  (1  -  cost  ) 

•fyj(t)  «  ^  sin  t  sin  (1  -  cos  t  ) 

fpCt)  •  -  ^  sin  t  cos  (1  -  cos  t  )  • 

■  ^ 


(111a) 


(Ulb) 


In  eq. (ill),  we  have 

ft(t) 


vd.th 


cos  t 

sin  t 

Sint 

cos  t  i 

-  cos  t 


♦  sin  t  y 


(2.1) 


(112) 


*  The  numerals  in  front  of  the  period  refer  to  ntimbers  in  the  main  text  while 
the  ntmierals  behind  the  period  designate  the  number  of  the  exaiiple  involved. 
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On  the  basis  of  this  latter  presentation  (2.1)  the  validity  of  eq.(49)  can  be 

directly  verified.  Therefore,  at  P  =  2n,  the  quantity  from  eq.{35)  vdll  /5A 
« 

be  the  zero  matrix,  i.e., 

(36.1) 

(30.1) 


and 


as  well  as 


if  •  r 

^ (t)  •ig(t) 

2/.(t)  (31.1) 

In  addition,  eq.(37)  must  be  used  as  the  solution  matrix  of  eq.(3),  because  of 
eq.(il4.9)  (see  theoran  6).  Taking  eq.(2.1)  into  consideration,  this  will  yield 

i  .  .e  e 

e  .  ^  if  t 


(t)  -  a 


(37.1) 


8l.n(1  -  cos  t  )1 

I^Bin  t 

008(1  -  cost  )l 

—^slnTj 

Therefore,  it  follows  that 


/cos(^  -  cos  t ) 

sia(l  -  cos  t ) 

I-  sia(l  -  cos  t  ) 

cos(1  -  cost  ) 

l^-sia  t 

0 

L 

•6int 

r 

e 

=  2n. 

(12.1) 


In  the  case  (111a),  the  resonance  case  is  present  since  it  is  calcxilated 


that 


(resonance  subcase) 

T 

f  i  =  ^  (exceptional  subcase) 


(86.1a) 


*  Ey  substituting  eqs.(37.l)  and  (12.1)  into  eq.(lO),  it  is  possible  to  confirm 
eq.(36.1)  a*  control. 
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Taking  eqs.(53)»  (75),  (76),  aaid  (31»1)  into  consideration,  we  obtain 


€  r 

»^(t)  .  -1 


(74.1) 


and  thus,  under  consideration  of  eq.(30.1). 


v^u>  »  fcj- 


(79.1) 


Calculating  out,  eq.(37.l)  in  view  of  eq. (86.1a)  will  yield 


/COSO  -cost)) 


sinCl  -  cos ' 


^BXJU  I*  /  .  V 

•  (t  ♦  c^). 


/ Sin(1  -  cost  )\  , 

“(cosO  -  cost)/* 


(113) 


In  accordance  with  theorem  15,  each  ccmqponent  of  r^^^(t),  with  increasing  t, 
will  again  and  again  take  values  of  the  order  of  t,  while  T  (t),  as  a  vector 


periodic  with  2n,  will  remain  restricted.  The  general  solution  of  eq.(lll) 


then  becomes 


r(t)  . 


(78.1) 


In  the  case  (lUb),  the  exceptional  case  exists  because  of 

[?  -  r  ** 

If  “  2  j  sinT  sin(2-2cosr)dr.-co8(2-2co8r)  ^  -o, 

7  *»  * 

2y  8inTcos(2-2cosT)dT.-sin(2-2co8T)|  -o. 

Analogous  to  eq.(113),  calcvilation  yields 


(86.1b) 


(>  /  cos(1  -  cost  )\ 

2  “I  ^  p-co8(2-2  costO-t^^l 

slnCl  -  cost  )J  ^ 

I  /  slnO  -  cost  )\  _ 

-  I  .  ^  I-  sin(2-2  cos  tJ.KS2]. 

^  008(1  -  cost  )/  ^  -* 


(114) 


sin(2-2  cos  • 


The  general  solution  of  eq.(lll),  con?5osed  in  accordance  with  eq.(78.l),  thus 


has  the  period  2n 


Im  (Re  theorem  15,  principal  case,  and  theorem  8), 


The  system 


possesses  the  matrix 


■  (1  ♦  COB  t  )x^  ♦  (2  ♦  Bint  )X2 
}x^  .  -(2+Bint)x^  ♦  (1  ♦  coBt)x2  - 


/l  ♦  COB  t  2  ♦  Bln 

\-2  -  Bln  t  1  ♦  COB  y 


■°^coa  ^(1*2t)-co8^ 
8ln^(1+2t)  -  cos  ^ 


|(1+C08  t)  ^ '4'(2«Bin 


(115) 


(2.2) 


with  3  from  eq,(ll2)  and  satisfies  the  condition  iu9)»  As  in  the  above  example, 


we  find 


so  that 


(35.2) 


(36.2) 


since  derivation  is  made  by  the  method  of  power  series  which  consists  in  expand¬ 
ing  the  exponential  function 

-V 

Consequently,  in  the  diagonal  matrix?,  we  have  [see  eq.(l6)3 

Thus,  the  principal  case  is  involved  here.  In  addition,  we  calculate 


l^(t)  -  B^^  ^  sint)^.  ^(1  ♦  2t  -  COS  t  ) Cf 


co8(l-f2t-cos  t)  1  sin(1+2t-cost)  I /e 


-  sin(1+2t-co8 1)  I  cos(1-f2t-co8t)/  I  o 


4  sin  ■ 


(37.2) 


from  which,  by  means  of  eq.(lO),  we  again  obtain  eq.(36.2).  Analogous  to  / 57 

eq.(113),  we  find 


w 

■£  (t) 


r  (t)  - 


/  C03(1+2t-COS  t)  1  ^ 

V  -  8in(1+2t-coa  t)  J  ' 

/  sin(Wt-co8  t3\  ^  ^ 

I  cos('l+2t-co8t)y  2 


(116) 
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With  Cl  =  -1,  the  quantity  (t)  becomes  periodic  with  2n  while  at  ca  =0  the 
vector  thus  also  the  total  vector  (78)  becomes  trivially  periodic. 

Example  3.  (Re  theorem  16,  resonance  and  exceptional  case,  and  theoron  8), 
Consider  the  system  of  differential  eq\iations 


with 


f  xjj  ■  -  x^sint  ♦  X2  ♦(sin  t  ♦  C08t)e*^“* 
I  x^  ■  -  X2Sint  ♦  ^2^*) 

r  f,(t)  -  ^ 


(117) 


and,  respectively 


« 


(t)  - 
(t)  -  o. 


cos  t  -  1 


The  coefficient  matrix  becomes 


sin  t  I  1  ♦  (sin  t  ♦  cos  t )  ^ 


~  sin  t 


(117a) 


(117b) 


(118) 


if  we  define 


■  —  Sint  ^♦^♦(sin  t  ♦  cos  t  • 

’  •(: } 


(119) 


Here,  ?t(t)  satisfies  the  assunption  (/+9).  We  can  then  calculate  directly 


■  (l  'J' 


so  that,  because  of§  =  0, 


(; ;? 


(35.3) 


(36.3) 


Consequently,  ^  is  an  elementary  component  of  the  rank  m  =  2,  with  the  eigen¬ 
value  A  =  1,  Further,  we  obtain 


with 


where  eqs.(30)t  (35»3)»  and  (36.3)  have  been  taken  into  consideration,  means 
of  X 

the  following  system  of  equations  is  obtained: 


-tc*  5(t)  (52.3) 

with  [see  eq.(32)] 

■^(t)  «  ^  *  (t)  ^  <tf  •  (53.3) 

Jn  the  case  of  eq.(117a)  and  in  accordance  with  eqs.(120)  and  (117a),  we 


have 


.  L-cos  t  ' _ ^sin  t  -  2  cost  +1  -cos  tj/_sin  t  -  1  | 

- /  5 - L 

\  o  '  ''  -  / 

•(-:)■ 


(53.3a) 


Because  of 


(92.3a) 


the  resonance  case  is  present. 

From  eq.(52.3)  [see  also  eq.(35.3)]»  we  then  obtain 

.  -  ^ 

-  V2  e 

«  1  ,  V2  •  t:  +  Cp  , 

-  1 

V  a  fir  ♦  C-jt  -^6  t  ♦  C  y.  f 

*  1  t  ' 
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(121a) 


so  that,  according  to  eq.(5l»3)» 


'(j  ♦  C2  t 


♦  e  t  ♦  c 


cos  t  —  1 


(t  +  Cp) 

In  accordance  with  theoran  16,  r(t)  takes  values  of  the  order  of  magnitude  t^ ; 
in  accordance  with  the  lanma  of  theorem  16,  Vi(t)  and  thus  also  xi(t)  again  /60 

and  again  are  of  the  order  of  magnitude  t^,  while  va(t)  and  xzCt)  are  of  the 
order  of  magnitude  t. 

Analogously,  in  the  case  of  eq*(117b),  we  obtain 

■^(t)  .  0-  (53.3b) 

7l.<TUr.o  .  (92.3b) 

;  ‘ 

so  that  the  exceptional  case  is  involved  here.  From  eq.(52.3)  we  obtain 

VjJ  •  Vg  •<  1 

▼2  ■  ®  »  ^2  ■  ®2  » 

-  (1  ♦  C2)t  ♦  , 

and  thus  also,  in  accordance  with  eqs.(5l»3)  and  (120), 

/ 


r(t)' 


^^coet  —  1_^ 


C2  • 


COS  t  -  1 


-2]c2  + 


(121b) 


Obviously, 

f(t  +  2n)  =  p(t)  is  valid  for  cg  =  -1 
Example  A.  [Re  theorem  15,  resonance  and  exceptional  case,  at  general 
«(t)]. 

Consider  the  system  /6l 


xjj  "  (1  ♦  sin  tj  +  Sint  X2  + 

x^  «  Sint  X2  + 


(122) 


iidth 


1  -  CO«t 


and,  respectively. 


f  f-,(t)  -  -  y  « 
[  f2(t)  . 


f^(t)  - 


(sin  t  -f  cos  t} 


(122a) 


rpCt)  *  o  • 


(122b) 


The  coefficient  matrix 


1  Sint  I  Sint 


(123) 


is  rot  commutative  either  with  St(T)  or  with  f  8(('>‘)dT,  as  can  be  readily  checked 

b 

by  calculation*  From  the  second  equation  of  the  system  (3)  we  find,  together 


with  eq.(l23), 


IzCt)  - 


.1  -  cost 


and  thus  from  the  first  equation  of  the  system  (3), 


t  *  ^  -  cost 


*  ^2  f  djrj  ■ 


t  1  -  cost  t  1  -  cost/  -t  cost  ♦Sint  1] 

-  c^s  e  ♦  c^s  e  ^-e  5  ♦ 


SO  that  we  can  select 


^  -  cos  j  (sin  t  ♦  cos  ^ 


Hence 


•4?  -i^'''(o)  -l^(2T)  - 


’  i\h‘' 


o  1/  \  o 


?»■  > 

e  o 


1/  \o  1 


(7.4) 


(lo.z.) 


The  matrix  Sp  is  partitioned  into  the  two  elementary  conponents  with  the 

3TT 

eigenvalue  =  e  and  ^  with  X3  =  1, 


Further,  we  calculate 
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-“t 


.^1 . *  - '« 


j  (sint'*'  cost  ‘  1, 


(12.4) 


In  the  case  of  eq.(122a),  because  of 

IT 

J  }^cc>  (86.4a) 

0 

the  resonance  subcase  is  present  for  ^  and  the  principal  case  for  5|Si .  We  then 
use  the  method  of  variation  of  the  constants*  for  further  calculation.  In  the  ' 
argument 


we  then  have 


;*(t)  -  - ««« *  (**)*" 

Z*2Ct)  -  s'*  -  /-  ?  (Sint*  cost)) 


o1 


(9.4) 

(124) 


(t  *  c^2). 


The  quantity  |ja(t)  again  and  again  assumes  values  of  the  order  of  magnitude  t 
while  |‘i(t),  on  selection  of  cqi  =  0,  is  periodic  with  2n. 

In  the  case  of  eq.(122b),  because  of 


J  €-r> 


I  so 


(86.4b) 


*  On  substituting  eq.(9)»  under  consideration  of  eq.(8),  into  eq.(l),  we  will 
have 

r'.  igv,  , 

i.e.,  according  to  eq.(12) 


t"  (t)  s  J  ^  (Ti  Affi  • 


this  will  yield 

For  the  sximmands  of  eq.(9)  in  view  of  eq.(8),  this  furnishes 
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the  exceptional  subcase  exists  for^a  end  the  principal  case  for  As  above. 


eq.(9*4)  is  obtained  vdth 


1  -  cost 


»  e 


r*(t)  -  s'*  - 


' -  2  (sint  ♦  cos  t 


(125) 


'o2  • 


Here,  r*(t),  for  any  Coa ,  has  the  period  2rT  while  rf(t’)  has  this  only  for  -  o. 
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PART  IV 

THE  RESONANCE  CASE  IN  LINEAR  ORDINARY  DIFFHIEJJTIAL 
EQUATIONS  OF  THE  ORDER 


Section  1.  General  Considerations 


In  an  earlier  paper  (Bibl.l),  we  investigated  the  resonance  case  in  systems 
of  n  linear  ordinary  differential  equations  of  the  first  order  with  periodic 
coefficients.  Here,  this  theory  will  be  applied  to  linear  ordinary  differential 
equations  of  the  n^^  order  with  periodic  coefficients.  Consequently,  let  the 
differential  equation 

L[30»  ajj(t)x  -  f(t),  (l) 

be  given  in  which,  for  all  coefficient  fvinctions  and  for  the  function  f(t), 
reality,  continuity,  and  periodicity  with  the  period  P  are  assumed: 


a^(t  ♦  P)  •  a^(t)  f(t  +  P)  «  f(t). 


Using  the  notations 


^  2^  "  ^*****»  _x  ■  X 


(n  -  1) 


the  differential  equation  (l)  is  transformed  into  the  system 


with 


r  '  ■ 


0  1 


•  0 


«  >•-»  t  '  1 


To  the  homogeneous  differential  equation  conjugate  to  eq.(l) 

*  a^(t)y^””"'\...+  an(t)y  -  o 

the  following  system  of  differential  equations  will  then  correspond 
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exists,  having  a  matrix  4(t)  periodic  with  P  and  a  constant  matrix  ft  which 
latter  is  written  in  the  Jordan  normal  form  [see  (Bibl.l),  eqs.(27),  (28),  (29)], 
from  which,  for  the  fundamental  solution  matrix  (13),  the  following  presentation 
is  obtained  [see  (Bibl.l)  eqs.(3l)  and  (32)]: 

Bits  Lf^ ‘iff,.- 

We  will  treat  the  general  case  [see  (Bibl.l)  eq.(27)]  that  the  matrix  ft  is 
partitioned  into  s  elementary  coirponents 


>  y  •  •  • • >  A 


(16) 


with  the  orders  my,  in  which  case  the  following  is  assumed  to  be  valid 


«  o  for  *  -  1,...,  f 

f  o  far  ...  ,8 


(17) 


(p  =  0  and  P  =  s  is  admitted). 

For  the  general  solution  x(t)  of  eq.(l),  using  the  method  of  variation  of 
the  constants  [see  (Bibl.l),  eq.(9)]  and  the  footnote  on  p.5S,  we  obtain 

/20 

»{t)  -  •  I  */.  ct)  fCt)  4t  (^) 

[see  eqs.(5),  (10),  and  (13)]  which  can  be  resolved,  in  the  form  of^ 


into  the  components 


*(t)  ■  21  *  *(*) 


(19) 

(20) 


*  We  changed  the  symbol  into  '’x  from  that  used  in  our  first  paper  (Bibl.l) 

so  as  to  prevent  confusion  with  the  derivatives. 


w 


I 


(21) 


with  the  sumnatlon  indices  [see  (Bibl.l),  eqs.(66)  and  (70)]: 

(y)  -  ^  *  •••  * 

For  the  components  ''x(t)  with  v  =  p  +  l,  , . , ,  s,  we  can  assvmie  that  they 
are  functions  periodic  with  P.  These  are  imiquely  determined  [see  (Bibl.l) 
theorem  16], 

In  the  case  v  =1,  2,  P,  the  following  is  obtained  in  accordance  with 

our* first  paper  [(Bibl.l)  eqs.(109)  and  (106)]: 


where  ''(^(t)  has  the  form 


•4  ^ 

(22) 

for/^.o,! . .^-1, 

(23) 

*e  f  tr*  rt)j 

’  fy»fr  / 

(24) 

The  v^^j4.Y(t)  defined  elsewhere  [(Bibl.l)  eq.(99)]  are  of  no  inportance  here,  /71 
It  is  further  found  that  the  constants  aside  from  depending  on  depend 

only  on  the  difference  m-  -  y;  thus,  using 


y  y 


(25) 


eq,(23)  can  be  replaced  by 


'  tmO 


0,1, 1; 

with  the  procedure  being  the  same  for  eq.(24).  The  quantities  ^  ~ 

-  Y  >  0  are  arbitrary  integration  constants,  whereas 


(26) 


V/*  “  P  “cy^  0,1,...,*^  -  1)  (27) 


with  [see  (Bibl.l),  eqs.(86),  (92),  (93)  and  note  ntvj  (t)  =  Ztvj  (t)  according 
to  eq.(lO)] 


W 


frt 


ftTtUr  . 


(28) 
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In  the  resonance  subcase  (v  =  1,  a)  the  quantity  '^xCt)  will  always  be  a 

polynoniial  of  the  degree  my  with  the  coefficients  ''S^(t)  which  is  periodic 
with  P,  while  the  coefficient  '*'3g(t)  of  the  highest  power  is  obtained  from 
eqs.(l4),  (26),  and  (2?)  as 


-  p  <  *'  -  'I . •')•  (29) 

which  means  that  it  is  an  eigenfvmction  of  the  homogeneous  differential  equa¬ 
tion  (6)  periodic  with  P.  If,  in  the  exceptional  subcase,  '’d0v  with  the  small¬ 
est  3v  >  0  is  the  first  nonvanishing  coefficient  (25),  we  will  obtain  in  an 
analogous  manner,  provided  that  0v 

/72 


for 


(30) 


a) 


A  <'~y 


In  what  follows,  eq.(22),  omitting  the  highest  t-powers  with  vanishing  coeffi¬ 
cients,  will  be  written  in  the  following  form: 


%  .Ml  ^  »*^«"iyforV»  1,,.,,^ 

for  /■  ^4  . 


(31) 


A  solution  (19)  of  eq.(l)  which  is  represented  in  the  form  (3I)  for  u  =  1, 
P  and  which  is  periodic  with  P  for  v=p+l,  s,  will  be  denoted  here 

as  the  "normal  solution”.  Such  a  normal  solution  can  be  written  in  the  form  of 

— / 


x(t) 


.  y  ^  w  i _ 


(32) 


where 


W  .  i-*,  *0) 


(33) 


is  used  auid  where  the  ^  are  con^josed  of  the  '"^(t)  in  eqs.(26)  and  (2aV  to 
yield  functions  periodic  with  P. 
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It  should  be  eophasized  again  that,  in  eq.(32),  only  one  special  solution 
of  the  principal  subcases  is  used,  namely,  the  solution  periodic  with  P,  whereas 
in  the  resonance  and  exceptional  subcases  the  conplete  solution  is  used  with  my 
arbitrary  integration  constants*  /73 

In  view  of  eqs*(29)  and  (30),  the  following  theorem  can  be  confirmed: 

Theorem  It  If,  in  eq.(33)*  («  >  0,  then  the  factor  Yo(t)  of  the  highest 
power  t  occiirring  in  eq#(32;  will  be  a  nonidentically  vanishing  solution  of  the 
homogeneous  differential  equation  (6),  periodic  with  P. 

It  is  directly  obvious: 

Theorem  2:  If  the  resonance  subcase  is  valid  for  v  =  l,  ,.,,a>0,  the 
order  of  the  powers  of  x(t)  in  eq.(32)  will  be  at  least 


and  at  least  one  solution  x(t)  will  exist  for  which 


This  latter  statement  follows  directly  from  a  consideration  of  eqs.(26),  (27), 
and  (31 )• 


A  successive  differentiation  of  eq.(32)  to  t  yields 


(35) 


where  the  formally  written  functions  Tg®^(t)  with  0  <  0  must  be  replaced  by 
zero.  From  this,  the  following  is  directly  obtained: 

Theoran  3:  If  To(t)  in  eq.(32)  is  not  constant,  all  derivatives  x^‘‘^(t)  of 
eq»(32),  with  k=l,  ...,n-l,  have  the  same  sequence  of  power  incranent  t**^ 
as  x(t)  itself*  If,  conversely,  To,  Ti,  **.,  are  constant  whereas  Y^(t)  is 
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not  c<mstant»  the  power  orders  of  the  derivatives  of  x(t)  decrease  successively 
by  1  down  to  (t)  and,  from  thai  on,  ranain  constant  equal  to  oo  - 

Section  2.  The  Case  a,(t)  #  0  Z?/,, 

Considering  that,  in  the  case  of 

“  (36) 

the  hoBiogeneous  differential  equation  (6)  in  eq«(l)  cannot  have  constant  solu¬ 
tions,  it  follows  directly  from  theorems  1  and  3^ 

Theorem  L:  Under  the  condition  (36),  all  derivatives  x^''^(t)  (k  =  0,  1, 

...,  n  -  1)  of  a  normal  solution  (32)  of  eq.(l)  always  have  the  same  power  order. 
Note;  The  only  exceptional  case  of  theoran  4  can  possibDy  be  the  case  of 

x(t)  •  Const.  -  C  ^  o  (37) 

which  occurs  for 

f(t)  -  (38) 

since,  in  that  case,  x^*^^  (t)  =  0  for  ■6^1.  It  is  worthwhile  to  groMp  this 
trivial  but  interesting  special  case  with  the  general  considerations. 

Primarily,  the  resonance  subcase  must  not  occur  here  since  then  terms  with 
t-powers  would  necessarily  enter  in  eq.(32).  Consequently,  the  assunption  of  a 
solution  ^a(t)  of  eq.(ll),  periodic  with  P,  for  which  [see  eqs.(28)  and  (38)] 

T 

j  fO  (39) 

o 

is  valid,  must  be  continued  vqp  to  contradiction.  This  contradiction  is  obtained 
directly  from  the  first  and  last  equations  of  the  syston  (10),  in  view  of 
eq.(39),  since  iZo(t),  and  thus  also  3a (t),  does  not  have  the  period  P  [  see 
eq.(l2)]. 

Consequently,  for  v=l,  ...»  P  the  exceptional  case  exists  while  for  v  = 


=p+l,  ...,  s  the  principal  case  is  pres«it;  here  P  -  0  and  p  -  s  is  admis-  /75 


sible 


First,  let  us  calculate  the  vector 


'« 


for  the  special  case  involved  here*  According  to  another  paper  [(Bibl.l), 
eqs.(85)  and  (31)]  and  to  the  notations  in  that  same  paper  and  in  eq.(76)  (see 
footnote  on  p.32),  we  obtain 


# 


'l  ,  T  . 


Consequently,,  taking  eq.(38)  as  well  as  the  first  and  last  equations  of  the 

system  (lO)  into  consideration,  we  have 

r  •  1  r  . 

#  '  •  ! 

L  J  [  ^  j 

in  accordance  with  eq.(l6). 

We  then  define  the  row  vectors 

X  "  (Xy,  ’  ’  I'fr.a  )  ,  j/  «  >■  .*crj)- 


According  to  eq.(15). 


c  =  I*'  ,7-'V  - 


SO  that  the  folloviing  is  va^-id  for  the  correspordirg  column  vectors; 

-  (if.  )  -  e  iJy  '*»  • 

Equation  (7U)  (Bibl.l),  namely, 

ic/  «  XV  t 
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as  readily  verified  by  eqa^iku)  and  (4l)»  bas  the  solution 


^  t  y-  1,2 . 8.  (^+6) 

periodic  with  P.  Thus,  eq.(32)  according  to  [(Bibl.l),  eqs.(79)  and  (77)J  can 
be  written  in  the  following  fom: 


4  ^ 

x(t)  -  c 


Mtif 


where,  in  the  last  transformation, 


*■1 

=  6  was  used  in  accordance  with  eq.(l5)» 


Section  3»  The  Differential  Equation  Reduced  in  Order 


/77 


Since  the  case  (36)  in  the  previous  Section  was  taken  care  of,  we  can 
assume  below  that^'*^ 

a  a  •••  •=  a  o,  |  o, 

1  ^  n  -  1 

5(t)  -  x(J)(t)  ,  ?(t)  -  y<J>(t)  , 

eqs.(l)  and  (6),  respectively,  are  transformed  into 

LfyJ-  a^(t)  ...  ♦  -  r(t), 

L[yJ  a  a.,(t)  ...  ♦  *  °  • 

The  pertaining  adjoint  homogeneous  equation  will  then  be 


•s  ..(n-3»  A-a-V 

z  tf-’) 


•i-i 


(t)z  •o. 


(48) 


(49) 


(50) 

(51) 

(52) 


Each  solution  of  eq.(52)  simultaneously  is  a  solution  of  eq.(9). 

Conversion  of  the  differential  eqiiations  (50),  (5l),  and  (52)  into  the 
corresponding  differential  equation  systems  of  the  first  order,  together  with 


*  The  trivial  case  j  =  n  leads  to  the  differential  equation  x^°^(t)  =  f(t)  with 
the  reduced  equation  [see  eq.(49)  and  (50)]  x  =  x^"'Ct)  "  f(t).  The  resonance 
case  or  the  exceptional  case  is  present  depending  on  whether  the  mean  value 

1  f 

-  !  f(T)dT  is  not  equal  or  is  equal  to  zero. 

^  0 


will  yield 


and 


^  m  ft  f  ^  , 

^ >  &  n>  9 


with  the  matrix 


i  » 


«> 


•  1 
,  .0^(0,  _*,«> 


which  is  correlated  with  the  matrix  3t(t)  as  follows  [see  (Bibl»l),  eq»(23)] 


•  1 

0 

• 

(T 

X 

(T 

a 

cl“> 

0 

■» 

Now  let 


be  a  fundamental  solution  matrix  of  the  reduced  homogeneous  differential  equa¬ 
tion  (51)  so  that,  for  B(t),  the  following  representation  [see  eqs.(l4)  and 
(19)]  applies:  ^  ^  ^  dr 

i .. ±  'i-  ■«  -  2  -« -  2  /  (6( 

Vm1 

with  the  fundamental  subsystems: 


1  *  *  * 

^C<J 

»<v>  ... 

i>) 

ti-y 

-  - 

(*-»> 

I 

IT  - 


and  the  periodic  submatrices: 


a 

•• 

0  U)  S 

cr 

• 

. 

IT 

m 

Again,  let  us  assume  [see  eq.(17)]  that 

o  fox 
f  o  fox 


o  r or  f  * 

o  for  ^^•f  +  1«***«^t 


where,  as  always  in  special  cases,  P  can  be  equal  to  zero  or  s.  In  that  case, 
the  matrix  ^t)  can  be  expanded  in  the  following  manner  to  a  fxmdamental  solu¬ 
tion  matrix  ?)(t)  of  eq.(6)  [see  eq.(58)]j 
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Here,  the  sqxiare  matrix  tn  (t;  is  successively  constructed  by  integrating  the 
first  row  of  the  matrix  ^t)  in  which  case  the  integration  constants  can  be 


arbitrarily  selected. 

Next,  we  partition  the  matrix  (t),  analogous  to  eq.(6l),  in  the  form  of 

M 


into  the  square  submatrices  S)  • 


Then,  we  can  formulate  the  following  theorem. 

Theorem  5:  If  the  reduced  differential  equation  (51)  has  exactly  P  inde¬ 
pendent  solutions  y(i),  7(2)  >  •••»  periodic  with  P,  then  the  differential 

equation  (6)  either  will  have  also  exactly  p  independent  solutions,  periodic 
with  P,  or  else  eq.(6)  will  have  exactly  p  +  1  independent  solutions  periodic 
with  P.  In  this  case,  the  following  is  valid  [see  the  notations  of  eq.(2l)]: 
l)  If,  for  all  y(  v)  (^)  •••»  P)» 


Si 

J  .  o, 

applies,  then  eq.(6)  will  have  exactly  P  +  1  independent  solutions  periodic 


with  P. 

2)  If,  conversely,  for  at  least  one  y(^)(t)  (v  =  1,  p),  the  following  mean 


value 
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(67) 


l)  Let  •  o  (for  y-  [see  eq.(66)]. 

By  j  integrations  of  the  functions  ^  which  case  the 

integration  constants  must  be  determined  each  time  in  accordance  with  the  auxil¬ 
iary  theorem,  exactly  one  function  y(y)(L)  having  a  mean  value  of  zero  and  being 
periodic  with  P  will  be  obtained  for  each  y^ (t).  Accordingly,  based  on  the 
trdvial  solution  yo(t)  =  0,  the  fionction  yo(t)  =  1  will  be  obtained  as  a  further 


solution  of  eq.(6)  periodic  with  P.  That  this  solution  is  independent  of  the 
above-defined  solutions  y(^)(t)  (''  ~  2,  6)  follows  from  the  fact  that 

all  these  y(y)(t)  have  the  mean  value  0,  while  yo(t)  has  the  mean  value  1.  That, 
in  addition,  the  y(v)('*')  ('^  "  •**»  also  mutually  and  linearly  inde- 
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pendoit  follows  from  the  linear  independence  of  the  functions  y(v)(t-)  after  j 
differentiations • 


2)  Conversely,  if  not  all  y(v)(''  =  1»  •••,  P)  have  a  mean  value  of  zero,  new 
solutions  can  be  obtained  by  linear  combination  which,  again,  will  be  denoted 
by  ^1)  »  •••»  7{o)  i  that  can  be  so  defined  that,  for  exairple,  y(i)  has  the  mean 
value  1  while  the  remaining  y(v)(''  =  2,  ...,  p)  have  the  mean  values  0,  i.e., 


•  r  *  , 


(ym  /,  •  -  A)  • 


(70) 


By  J  integrations  of  the  functions  y(v)(v  2,  ...,  p)  and  in  accordance  with  783 
the  auxiliary  theorem,  again  P  -  1  solutions  y(v)(t)  of  eq.(6),  periodic  with  P 
and  having  a  mean  value  of  zero,  will  be  obtained  which  are  mutually  and  linear¬ 
ly  independent.  An  integration  of  y^jj  (t)  would  yield  a  solution  ycD  (t)  not 
periodic  with  P.  However,  the  trivial  solution  yo(t)  =  0  again  leads  to  the 
solution  yo(t)  =  1  periodic  with  P  which,  together  with  the  functions  y(2)(t), 
•••»  y(P)(b),  fomis  a  system  of  P  linearly  independent  solutions  of  eq.(6), 
perdodic  with  P. 

Finally,  it  is  easy  to  demonstrate  that  eq.(6)  can  have  no  further  solu¬ 
tion  periodic  with  P.  Let  y(t)  be  any  nonconstant  solution  of  eq.(6)  periodic 
with  P;  then,  y(t)  =  y^^^(t)  will  be  a  solution  of  the  reduced  differential 
equation  (51)  periodic  with  P  which,  in  addition,  has  a  mean  value  of  0  because 
of  the  differentiation  process;  consequently,  y^'*^(t)  must  be  linearly  com- 
posable  of  the  already  known  solutions  y(V)(b)  ('^  ^  1^  ••*,  P). 

It  should  also  be  mentioned  that  the  solution  vectors  Mt)  of  eq.(56)  must 
be  conqjleted  into  the  solution  vectors  i(t)  of  eq.(ll)  by  a  successive  differ- 


entlation  process  in  accordance  with  eq.(lO),  toward  components  with  smaller 


indices. 


Section  4.  Definition  of  a  Special  Normal  Form  of  the  FundamCTital 
System  of  the  Reduced  Homogeneous  Differential  Equation 


Unfortunately,  for  solving  the  reduced  syston  of  differential  equations 
(54),  it  is  not  sufficient  to  merely  bring  the  matrix  ft,  introduced  in  accord¬ 
ance  with  eq.(l4),  to  the  Jordan  normal  form,  in  view  of  the  fact  that  the  main 
purpose  of  the  theoretical  consideration  is  a  discussion  of  the  system  (4)»  For 
this,  it  is  necessary  to  obtain  the  square  matrix  ?)*(t)  occurring  in  eq.(64)  in 
as  concise  and  sinple  a  form  as  possible.  For  this  purpose,  the  fiindamental 
systan  (59)  of  solutions  of  the  reduced  partial  equation  (5I)  is  brought  to  a 
special  normal  form,  which  we  will  characterize  by  the  following  properties  of 
the  matrices  (62)  occurring  in  eq.(60): 

Definition; 

1)  At  ^  0»  either  the  following  mean  values  are  valid  for  =  (v),  ...,  [v] 

T 

p  J  <5^(t)dt  -  o  for  all  u  (71) 

or  an  index  0  ^  iv  ^  niy  -  1  (mv  =  [v]  -  (v)  +  l)  exists,  so  that  the  following 


2)  The  elementary  conponents  are  so  arranged  that,  for  v  =  1,  ...,  sub¬ 
scripts  iv  exist  whereas,  for  v  =  A.  +  1,  .. .,  p ,  all  mean  values  (71)  are  equal 
to  zero.  For  v  =  p  +  1,  ...,  s,  we  have  ofy  =  0;  here,  the  elementary  components 
are  not  restricted. 
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3)  The  orders  my  increase  monotonically  for  v  =  1,  X; 

m  >  •  I'  “  "  t  •  «  A  « 

SO  that  the  ineqiJalities 

A  ^ 

ly  >  iy_-i  1  ®y-1  ~  -1  2  , .  •  .  ,  A  . 


(73) 

Z8^ 

(7/.) 


are  valid.  The  elonaitary  conponents  for  which  all  mean  values  vanish  are  also 
arranged  in  accordance  with  increasing  values  of  niy  except  that,  in  this  case, 
only 


A  A  *  O  »» 


(75) 


applies. 

For  constructing  this  special  normal  form,  we  will  need  several  new  con¬ 
cepts  and  theorems. 

Let  [see  eq.(l4)  and  the  notation  (42)j,  under  the  assumption  of  Oy  =  0, 


f 


if'* 


i  i 


(76) 


be  a  fundamental  subsystem  belonging  to  the  v^b  elementary  conponent. 
all  mean  values  are  j 

'V  -  f  J  (t)dt  =  o  forA=  (v),.. 

0 

or  one  iy  exists,  so  that 


Either 

(77) 


M 


r  < 


o  o  for 
f  o  for 
arbitrarily  for 


(y<  <  /»<<■'»  -ti-y  ) 


Then,  the  following  applies: 

Theorem  6:  A  regular  matrix 


(^’S) 
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^5,  y 

*  C  ^  ^  C  '‘v  ^  *"*  ^  ^ 

rw  *'  '  I  tr J 


is  in  existence  [see  the  denotation  (58)J»  so  that  the  transformed  fundamental 
system  ihv^  has  the  property  1. 

V 

Proof;  Since  the  matrix  5  is  commutative  with  the  matrix 


/  .  . 
-f- 1  ^  T  ^  ^  T 


-i  -1  a 


^  '  ~~ 


y 

it  is  immaterial,  because  of  eq.(l/:(.),  whether  the  transformation  ®  is  applied 
to  the  y^^  or  directly  to  the  The  row  vector 

h  ^  (8] 


is  then  transformed  into  the  column  vector 


Ji 


If,  for  each  subscript  u,  the  conponents 


are  separated  into  the  constant  mean  value  [see  eq.(78)J  and  into  the  func- 
tion  $|j_(t)  of  a  mean  value  zero,  the  following  applies  in  a  readily  understand¬ 


able  notation; 


'•*  *  ~  V-  ytr'^ 


The  row  vector  Ev  is  then  to  be  transformed  by  the  last  summand  in  eq.(82)  into 
the  row  vector  !Rv^,  which  is  defined  by 
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This [see  eq.(79)]  leads  to  the  systan  of  equations: 


®  ■  *(y) 


O  ■  C(y)  *  •**  *  ®(»')+i  -'I  (»') 

■  *(.»)  “(»')*iy  *  ®(V)+I^“(y) 

®  “  ®(y)  '‘(y)  +  iy-»''  *  *“  * 


0  -  e 


(y)  “(y)  ♦  •**  ♦  ®(y)“(y) 


The  first  iv  eqmtions  are  automatically  satisfied  because  of  eq.(78).  The  re¬ 
maining  equations  successively  lead  to  solutions:  C(v)»  C(v)  +  i»  •••»  C[v]-iy» 
in  which  case  we  definitely  will  have  c^v)  f  0»  This  means  that  5  is  regular. 
[In  the  case  that,  instead  of  eq.(62),  all  ~  0  (see  eq.(77))  one  can  simply 

•  y 

pose  5  =  ®v  •] 

From  now  on  we  can  assume  that  the  fundamental  subsystems  (y(  v)  ,  •••, 
ycvj  ),  Tor  all  v,  have  at  least  the  first  property  of  the  special  normal  form. 

The  further  properties  can  be  established  by  means  of  the  following  theoron:  /88 
Theorem  7:  Let,  in  addition  to  the  fundamental  subsystan  (76),  another 


fundamental  subsystem  be  given: 


so  that,  in  contrast  to  eq.(74)> 


•k  "  i*' 


the  following  simultaneously  applies: 
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(89) 


The  right-hand  inequality  in  eq.(89)  is  trivially  satisfied.  Then,  by  super- 
position  of  these  two  systans,  a  fundamental  subsystem  ihv  analogous  to 
eq.(76)  can  be  formed  for  which  the  functions  •••#  ®cvf,  which  are  peri¬ 

odic  with  P  all  will  have  the  mean  value  zero.  In  this  case,  the  fundamental 

A  T 

systan  ihk  ranains  unchanged. 

Proof;  If  the  new  solutions  y^.*  of  eq.(5l)  are  formed  in  accordance  with 


Here,  eq.(89)  guarantees  the  correct  values  of  u  for  eqs.(90)  and  (91)  while 
eq.(88)  ensures  that 

(k)  <  (kH/x.  -  (fi-z-t  /  C*J 

is  applicable.  It  can  be  confirmed  that  the  functions  from  eq.(91)  have  /89 
a  zero  mean  value  since,  for  u  =  (v)  +  iy,  we  simultaneously  have  (k)  +  u  - 
-((v)  +  iy  -  j^j)  =  (k)  +  ik.  Thus,  theorem  7  is  proved. 

Now,  we  can  establish  a  certain  normal  form  by  also  satisfying  the  condi¬ 
tions  2  and  3  of  its  definition.  Each  elementary  component  will  be  denoted  by 
a  pair  of  numbers  (%,  iy)  or  (my,  my)  if  no  iy  exists.  In  addition,  we  calcu¬ 
late  the  differences  my  -  iy  and  define  a  sequence  of  number  triples  (nly,  iy , 
my  -  iv)  or  (iiiy,  niy ,  O).  First,  the  subscripts  v  with  the  smallest  iy  will  be 
defined.  Among  these  subscripts,  one  is  selected  for  which  my  -  iy  is  as  large 
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as  possible.  After  this,  the  elenoitary  components  or  the  number  triples  (m*, 
iv ,  fiv  -  iv )  are  so  rearranged  that  the  above  characterized  number  triple  will 
be  in  the  first  position,  after  which  we  set  the  corresponding  v  =  1,  Accoi*ding 
to  theorem  7,  all  other  elementary  coirponents  for  which  eqs.(88)  and  (89)  with 
k  =  1  are  satisfied,  can  be  so  transformed  by  superposition  with  the  first  ele¬ 
mentary  component  that  no  more  iy  will  exist  for  them.  The  already  determined 
first  elementary  conponent  will  ranain  tinchanged.  Since,  in  any  case  iy  ^  ii 
and  since,  for  iv  “  ii »  we  automatically  have  niy  -  iy  s  nii  -  ii ,  this  means 
that  the  possibly  remaining  number  triples  for  which  an  iy  still  exists,  will 
satisfy  the  conditions  iy  >  ii  and  my  -  iy  >  mi  -  ii .  Of  these  remaining 
triples  we  again  select  one  for  which  we  will  pose  v  =  2  so  that,  at  minimal  iy, 
the  my  -  iy  becomes  maximal.  Here  again,  a  superposition  in  accordance  with 
theorem  7  will  make  it  possible  to  transform  the  elementary  components  for 
which  eqs.(88)  and  (89)  with  k  =  2  are  valid,  together  with  the  above-defined 
second  elementary  component,  in  such  a  manner  that  no  iy  will  exist  for  them 
whereas  the  number  triples  for  which  an  iy  might  still  exist  will  satisfy  the 
conditions 

^2^  i-]  and  -  1^*2  “  ^2^  *1  ” 

We  continue  in  this  manner  and  finally  obtain  a  sequence  of  elementary  compon¬ 
ents  with  number  triples  (niy ,  iy ,  niy  -  iy)  for  v  =  1,  2,  ...,  X,  where  the  /90 

following  is  valid  for  v  =  2,  ...,  X  (naturally,  X  =  0  or  X  =  1  is  also  pos¬ 
sible):  ^  ^ 

from  which  also  ifiy  >  iSy-i  follows  directly.  Then,  for  ay  =  0  the  elementary 
conponents  with  the  characterizing  triples  (my,  my,  0)  for  v=X+l,  ...,p 
might  possibly  be  left  over,  which  can  be  so  arranged  that  niy  ^  myi  (for  v  = 

=  X  +  2,  ...,  p)  is  valid.  The  elementary  components  with  ofy  /  0  can  remain 
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unchanged.  This  finally  establishes  the  special  normal  form 


Section  5.  Completion  of  the  Special  Normal  Form  of  the  Fundamental 
System  of  the  Reduced  Homogeneous  Differential  Equation 
into  a  FunriajiiflntAl  System  of  the  Original  Homogeneous 
Differential  Eaiiation 


The  purpose  of  this  particular  Section  is  to  obtain  proof  of  the  following: 

Theorem  8t  Based  on  the  special  form  (59)  resp.  (60)  of  the  fundamental 
system  of  solutions  of  the  differential  equation  (51),  a  fundamental  systorn 

t^Ct)  -  ^(t)  •  ^(t)  (92) 

of  solutions  of  the  differential  equation  (6)  can  be  obtained,  in  which  case 

(t)  -  .  ^(t  ♦  i>)  -  ^  (t)  (93) 

but  where  the  constant  matrix  fi  does  not  necessarily  have  the  Jordan  normal  form 
[see  also  eqs.(10e)  and  (11/*.)]. 

The  proof  must  proceed  in  seyeral  steps.  Primarily,  we  establish: 

Theorem  9:  If,  in  the  case  Qfy  =  0,  eq.(7l)  is  valid  for  all  values  of  u  = 

=  (v),  ,,,,  [v]  as  well  as  in  the  case  ary  f  0*,  the  solution  submatrix  [see 
eq.(6l)3  of  eq.(5l)  can  be  thus  completed  by  a  matrix  m  into  a  solution  sub¬ 
matrix  SI  of  eq.(6)  such  that,  once  more, 

’'10  *1  f  ’'#*1  I  t 

‘  ^  '  (94) 


is  valid. 


Proof.  First,  let  us  consider  the  case  j  =  1.  So  that  the  direct  argument 


*  And  thus  also  Oy  f 


2kni 


[see  (Bibl.l),  eq.(20)  with  q  =  P]. 


[see  eqs.(8l)  and  (80)],  after  one  differentiation,  will  lead  to 

xr  It 

tit  •  c  e 


131 

(96) 


the  existence  of  the  following  differential  equations  for  <P(V)(t),  'P(v)  +  i(t), 
tptvj  (t)  will  be  necessary  and  sufficient: 


1 

A 

* 

1  - 

A 

*  1 

A 

from  which,  because  of  the  two  possibilities  of  theorem  9,  functions  'P(v)(6), 

•  ••*  (t)  periodic  with  P  can  be  successively  determined  for  Ov  ;  these  will 

be  defined  uniquely  in  the  case  ory  ^  0  whereas,  in  the  case  ofy  =  0  if  eq.(71) 
applies,  they  will  be  defined  uniquely  only  if  it  is  additionally  stipulated 
that  all  mean  values  are  zero.  In  the  case  of  j  =  2,  the  same  method  of  reason¬ 
ing  is  to  be  applied  first  to 

(98) 

instead  of  to  eq.(95),  after  which  the  iCpy^  is  determined  fix>m  eq.(95)  according 
to  the  same  syllogism,  by  substituting  eq.(96)  with  eq.(98).  The  procedure  is 
wholly  similar  for  the  remaining  values  of  j.  From  this,  eq»{9U)  of  theorem  9 
will  follow  if  ^  A 

t  ^  (99) 

is  also  taken  into  consideration  [see  the  equation  corresponding  to  the  second 
relation  (ZJ.)J. 

This  leaves  the  case  Oy  =  0,  at  validity  of  eq.(72),  to  be  considered.  Un- 

......  -  . .  .  T  . 


der  introduction  of  the  vectors  ny  with  the  con5)onents 

1  for  ^ 

e  •  ' 

„  for 


/93 

(100) 
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we  resolve 


XT  X  ' 

X  “  ♦  *- 


(101) 


As  in  the  proof  for  theoran  9j  we  can  first  determine  a  viniquely  defined  row 
vector  with  the  mean  value  0  so  that  the  derivative  (note^v  “^v) 

A  ^  ^ 

«  .  T  T 

of  icpv  •  e  is  exactly  i^Pv  *  e  •  By  j  integrations  of  nv  *  e  we  ob¬ 


tain  the  row  vector 


,  fjrryT  i 


K-  K^J  * 


(102) 


From  this,  we  have: 


Theoran  10:  In  the  case  Ov  =  0  and  eq.(72),  eq.(94)  will  be  replaced 


with 


Yc).  ; 


o,...,o,  jr  ,  -^jnTr  v" 


0,...,0,  t,  jy 


Instead  of  eq.(l03),  we  can  also  write 


y  a;  * 


t*y“  ^ 


(103) 


(104) 


(105) 


Finally,  j  integrations  of  the  trivial  solution  system 


will  yield  the  row  vector 


.  (o,...,o)  (j  zeros  ) 


^“*0  “  . fWr) 


from  which  the  solution  subnatrdx 


■»  ♦•••  t_ _  r  3  r 

<ir«’ 

/r  ‘  =  ' 


(106) 


(107) 
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is  obtained.  By  adding  the  matrices  of  of  theorems  9  and  10,  we  form 


the  matrix  ^*(t)  and  write 


1 

— 

IT 

(106) 


Similarly,  we  fom  (see  theorem  10)*  the  quantity  6'''(t)  by  addition  of  the  ma¬ 
trices  '*fe*(t)  according  to  eq.(104)  and  then,  taking  eq.(107)  into  considera¬ 
tion,  put  •<» 

■to  in  . _ 

A  At 

This  will  yield  a  fundamental  solution  matrix  m(t)  of  eq.(6)  in  the  form  of  /95 
eq.(92),  where  already  the  second  equation  of  the  systen  (93)  possesses  validity. 

To  furnish  a  complete  proof  for  theorem  8,  it  merely  remains  to  be  demon¬ 
strated  that  eq.(109)  can  also  be  written  in  the  form  of  the  first  equation  of 
the  system  (93).  This  is  automatically  the  case  if  no  ay  "0  with  the  validity 
of  eq.(72)  are  present,  i.e.,  if  the  theorem  10  need  not  be  used.  Then,  the 
matrix  §(t)  of  eq.(109)  will  already  have  the  conventional  form: 


(no) 


with 


ir 

(111) 


In  the  case  of  theorem  9,  we  must  use  '’^5"(t)  -  0, 


and 


e  -f 


(112) 


of  the  order  nb  =  j.  Here,  has  the  Jordan  normal  form. 

Generally,  the  following  theorem  applies:  /96 

Theorem  11:  The  matrix  ^(t)  according  to  eq.(109)  has  the  form  of 


with 


For  the  meaning  of  X,  see  the  definition  given  in  Section  k*  In  the  case  X  = 
=  0,  eq.(ll6)  is  transfoimed  into  eq.(llO)  or  eq.(ll4)  into  eq.(lll). 
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with  a  1  in  the  row  occvipy  the  position  1  ♦  iy.  In  the  case  ^  >  0,  eq.(ll4) 
will  not  have  the  Jordan  normal  fom. 


with  (note:  Av  ~  for  v  -  1,  X) 


» 4-' 

»#  •  --^Oj  $  V  .  1 


6-»>!  ^ '  TTJTcFv^* 

•  >■•'<•«  1  i  *  j  ■■•■«  ^  .  ** 


•j-  •  -fO  rt—l  - -  f  ■■  •»£ - - - 


>1  ^  J 


rf-i 

(?«! 


1  J>  o>.  .  •  J  s 


i.e.,  according  to  eqs.(104),  (10?),  and  (115): 

V«’  *  X 

Hence,  eq.(117)  will  yield 


U)J 


A 

e  .  4 


.  Ji. 


Z22 

(118) 


(119) 


From  this,  on  the  basis  of  eqs.(ll6)  and  (II4),  the  following  equation  can  be 
read  off: 


■^(t)  .  ^(t)  ‘A 

Further,  in  accordance  with  eqs«(ll6)  and  (IOZ4,),  we  obviously  have 


(120) 


8' 


(121) 

However,  eq.(ll3)  follows  from  eqs.(l20)  and  (121)  [see,  for  exan?»le  (Bibl,3), 
Secb*3*A-3  • 

As  a  corollary,  we  should  note:  For  the  matrix  jp,  defined  in  another  /lOO 
report  [(Bibl.l),  eq.(lO)]  and  which  is  constant  in  accordance  with  the  same 
paper  (theorem  l),  the  following  is  valid: 

t  -  ^P)  -  9^^  .  (122) 

For  this,  besides  eqs.(92)  emd  (121 ),  only  the  second  equation  in  the  system 
(93)  as  well  as  eq.(ll3)  are  needed. 

Section  6.  Construction  of  the  Jordan  KlnmiAl  Form  for  /lOl 

Now,  we  will  have  to  bring  the  matrix  ft,  defined  in  eq.(llA),  to  the  Jordan 
normal  form  ff  by  a.  similarity  transformation: 

(123) 

The  pertinent  fundamental  solution  matrix  of  eq.(6),  according  to  eq.(lA)  and 
according  to  our  first  paper  [(Bibl.l),  eq.(l4)]  will  then  read 

.  (124) 

This  consideration  is  necessary  only  for  X  >  0  in  eq.(ll4);  at  X  =  0,  we  can 
put  5  =  ®.  The  collineatory  transformation  (123)  is  performed  in  j  individual 
steps,  in  which  case,  for  U  =  1,  2,  ...,  j  the  similarity  transformation 

-"A  ^  ■ -A  (125) 

will  produce  a  chain  of  matrices  ^ft  whose  first  link  (|Ji  =  0)  is  formed  by  the 
matrix  °ft  =  ft  and  whose  last  link  is  formed  ty  the  matrix  ^ft  =  ft° .  In  this 
case,  each  matrix  ^ft  will  have  the  following  form: 
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r 


(126) 


where  the  matrix  ^  order  (j  -  M-),  while  the  matrices  = 

=  V  =  0,  1,  X)  have  the  orders* 


(127) 


For  the  qiaantities  which  occur  in  eq.(l27)  and  are  not  yet  defined,  the  follow¬ 


ing  must  be  set  4.  ^  t  m.  <  \  m. 

®o  •  ^0  •  ®  *  “A  (128) 

Obviously,  °niv  =  my.  The  matrices  ^fiy,  also  occurring  in  eq.(l26)  and  having 
^my  columns  and  j  -  IJ-  rows,  contain  a  1  for  ^^  <  iy  in  the  last  rov  and  in  the 
(1  +  iy  -  u )^^  column  and  also  contain  a  1  in  the  last  row  and  first  column  for 
iy  ^  M-  <  iy+i,  but  otherwise  only  zeros.  In  the  case  ki-  ^  iy+i,  the  quantity  fly 
is  the  zero  matrix.  Here,  the  undefined  ix+i  constitutes  no  restriction. 

Obviously,  has  the  Jordan  normal  form  while,  for  pi  =  0,  the  quanti¬ 

ty  =  Pi  according  to  eq.(llZ4.)  has  the  form  of  eq.(126).  It  should  be  noted 
here  that,  because  of  °nt)  =  0,  no  matrix  resp.  occurs.  In  what  follows, 

we  will  make  an  induction  from  -  1  to  p..  For  this,  we  define  a  subscript 

't(n)  for  each  p.  =  1,  2,  ...,  j  in  the  following  manner:  If  u  falls  into  one  of 
the  intervals  with  the  end  points  ii ,  is,  ...,  ix  (the  left-hand  end  point  is 

*  See  the  definitions  at  the  beginning  of  Section  4,. 
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Included  in  the  interval  while  the  right-hand  point  is  excluded),  we  will  dwiote 
the  left  interval  end  point  by  itCu).  Thus, 

For  the  not  yet  covered  values  of  we  assume 

-  o  lor  (129b) 

l  ^  for  ^  ^  *X  ’  n  OQr  ) 


The  following  relations,  derived  from  eq.(l27),  should  be  noted: 

where  eq.(2l)  had  been  taken  into  consideration,  as  well  as 


(129b) 

(129c) 

/m 


(130) 


'"’'v  •  -*  +  ^  'ii 


far  ^  >  L(^ 


(131) 


(132) 


It  is  useful  to  repeat  the  matrix  (126)  in  a  more  detailed  form  (for  u  > 


>  0): 


<! 

a 

N 

3 _ El _ EL 


(133) 
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The  matrices  5^  in  eq.(125),  to  be  constructed  next,  have  a  differing  /lOU 
structure  dep«r\ding  on  whether  the  indices  u  -  1  and  belong  to  the  same  or  to 
different  intervals  (129a,  b,  c). 

First  case :  We  have 


i.e.,  the  following  is  valid; 
or,  in  the  two  limiting  cases. 


In  this  case,  we  put 


/or 

for 


(134) 

(135a) 

(135b) 

(135c) 

(136) 


(137) 


orders  of  the  matrices  ^”^o»  •••»  in  the  matrix  [see 

eqs.(l26)  or  (133 )]•  The  matrix  ^(Sv*  as  already  shown  in  the  matrix  (137)  is 
a  null  matrix  if  the  matrix  is  also  a  null  matrix.  The  ranaining  matrices 

a 

5v  for  V  ^  A.  are  calculated  on  the  basis  of  the  transformation 
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(138) 


y, 

from  which  it  follows  specificedly  that  also  the  matrix  -  1)  is  a  zero 

matrix. 

Relative  to  with  v  >  tCii  -  1),  the  1  in  the  last  row  in  ^Gv  is 

shifted  one  place  toward  the  left.  In  the  case  of  ^"^f<G|j,,  the  ones  in 
with  V  >  t{}i  -  l)  are  supplemented  by  a  1  which  is  shifted  by  one  unit  toward 
the  left  top.  Since  ^Gv"^  is  formed  from  ^Gv  by  replacing  the  ones  in  the  ^Gv 
by  -1,  all  ones  inG^^”^  in  the  (j  -  row  of  vanish  except  for 

the  1  in  the  first  column  of  -  1)  in  the  case  of  u  >  1. 


This  means  that  the  matrix  G, 


fi  '  has  the  property  that  the 


ones  standing  in  the  matrices  for  u  >  -  l)  are  shifted  by  one  place 

toward  the  left  top  in  a  diagonal  parallel  to  the  main  diagonal,  while  all  other 
elementary,  components  of  the  matrix  are  retained.  Specifically,  the  one  in 
the  lower  left  comer  of  the  matrix  -  l)  remains  fixed  for  >  1. 

Next,  we  perform  another  similarity  transformation  with  the  "permutation 
matrix"  SBjj,  which  is  produced  from  the  n-row  unit  matrix,  by  supplementing  /106 
the  rows  and  columns  having  the  numbers 


»  4  1  t  ■  ■  ■,  >  ♦It  ■  >> 


(the  next  following  column  exactly  has  the  1  in  question)  with  the  cyclic  trains- 


formation  matrix 


(139) 


In  the  case  <  ii  ri*e«»  '^(M')  ~  0],  we  must  put  25^  =  G  because  of  the  fact  that 


no  excess  1  need  be  eliminated  here.  It  should  be  noted  that 


(uo) 
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In  ^  which,  until  now,  had  been  in  the  (j  -  pt  ■*■  row  at  the 

+  1)^^  place,  has  been  shifted  to  the  place  directly  diagonally  above  the 
still  to  be  eliminated  1  which  remains  at  its  place;  all  matrices  in  the 
located  within  the  range  of  the  matrix  (139),  are  shifted  by  one  unit  toward 
the  left.  In  the  case  of  same  are  shifted  by  one 

unit  toward  the  top  while  the  1  in  the  (j  -  row  remains  xmchanged;  the  1 

standing  diagonally  toward  the  lower  left,  which  is  to  be  eliminated,  has  been 
shifted  in  the  same  colvunn  into  the  last  row  of  differing  from  the  corre¬ 
sponding  row  of  the  unit  matrix,  and  thus  changes  ,  by  one  first  row 

and  column  to  a  with  an  order  greater  by  one  [note:  -tCu)  =  ‘t'(u  -  1], 

Consequently,  all  in  all  the  new  similarity  transformation  with  shifts  the 
one,  still  standing  at  an  unwanted  place,  by  one  place  toward  the  left  top,  in¬ 
creases  the  order  of  toward  the  left  top  by  one,  shifts  the  matrices 

f  •••,  li-i  j -1  by  orie  place  toward  the  left  top,  and  finally 

decreases  the  order  of  by  one.  This  means  that  receives  the  order 

j  -  M.  while  receives  the  order  “'tc  (i-i ) +i  •  '^be  other  orders 

are  retained  [see  eq.(l27)J.  The  ones  in  the  (j  -  M- row  are  now  standing  /107 
as  stipiilated  by  eq.(l33)*  Consequently,  in  this  first  case  we  have  transformed 
the  matrix  with  5,,  =  S,,T15,,  in  a  similar  manner  into  the  matrix 

Second  case:  In  contrast  to  eq.(l34)»  we  assume 

■<  K^?  t  (iZtl) 

X  •  0  • 


and,  respectively. 


i  ^ S-i  {or-i<, 

"  i  \  > 

tor  • 


(U2a) 

(I42b) 

(I42c) 


Obviously,  more  accurately  than  eq.(l4l),  we  then  have 


s  and  ^  ~^mO-  (143) 

As  in  the  first  case,  we  first  use  the  similarity  transformation  with  the  ma- 
trix  **4*  As  a  result,  we  obtain  the  matrix 

^  ,  (,y., 

in  which  the  ones  are  standing  in  the  (j  -  row  but  [see  also  the  second 

equation  of  the  systan  (143)3  are  now,  for  v  =  -(.(M')  -  1  as  well  as  for  v  =  t(M.), 

shifted  to  the  first  colxann  of  Consequently,  the  matrix  differs  from 

the  matrix  only  by  the  one  at  the  place  -^(11)  -  1  in  the  (j  -  row  which 

u 

is  not  present  in  ft. 


For  eliminating  this  one,  the  following  transformation  is  introduced: 


(145) 


with  the  ’’superposition  matrix”  which  is  obtained  from  the  n- rowed  unit  /lOe 
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In  this  case,  the  unit  matrices  have  the  same  orders  as  the  corresponding  ma¬ 
trices  ,  ...»  in  the  matrix  (133)*  The  quantity  S3 is  obtained 

by  substituting  in  the  inserted  matrix  ^  *  Already  ft 


has  the  prescribed  new  orders  of  the  d-matrices  in  ft.  Specifically,  the  order 
of  has  increased  by  one  while  ^”^9°  has  received  an  order  lower  by  one. 
The  application  of  S3 


to  ft  eliminates  the  superfluous  1,  but  brings  the  /109 


matrix  -  to  the  same  place  at  which  is  located  in  the  matrix 

(146);  application  to  the  left-hand  side  will  again  eliminate  this 


•C.  •  -^1  ^ 


Consequently,  in  this  second 

(147) 


auxiliary  matrix  so  that  we  finally  obtain  ft 
case  we  must  put 

As  the  overall  result  of  our  considerations,  we  then  obtain: 

Theorem  12:  The  matrix  ft,  by  means  of  the  similarity  transformation  (123) 
with  J"  =  4r  //I 

is  transformed  into  its  Jordan  normal  form 

U: 


(148) 


'•<^1 


(149) 


Here,  the  matrices  (5,,  are  defined  by  eqs.(136)  resp.  (147),  with  the  auxiliary 

M' 

matrices  (137),  (I46),  and  the  described  The  elementary  components  for 
0  ^  V  5  X  have  the  following  orders  mv  [see  eq.(l27)j: 

-  my  +  Uin  -  Min  (j,v),  (150) 

where,  accordingly,  we  must  put  v  =  X  [see  eq.(l28)] 


Min  (d.ij^^-,)  -  0,  (151) 

for  V  =  0: 

m.i.o  (152) 


*  Ey  means  of  such  a  superposition  matrix  S,  the  proof  of  theorem  7  can  be  con¬ 
ducted  e^qjlicitly  and  sinply,  in  a  concise  manner. 
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In  the  case  of  =0,  the  elanentary  component  vdll  vanish.  At  v  >  X,  /IIP 
i.e.,  if  either  /  0  or  if  no  iv  at  all  exists  at  ofy  =  0,  the  following  is 
trivially  valid; 

m,  -  m,  ,  tpecifically.  -  j  .  (153) 

Another  remaric  should  be  made  here  as  to  the  correlation  of  the  rows  s 
and  s.  In  accordance  with  eq.(17),  s  denotes  the  number  of  elanentary  com¬ 
ponents  in  the  original  Jordan  normal  form  of  ft  so  that  s  also  will  be  the  num¬ 
ber  of  elanentary  components  of  eq.(l49),  which  means  that 


J  .  ‘/  »o  >  ® 

1  O,  if  -  o  . 


(154) 


is  valid.  It  is  obvious  (see  theorem  5)  that  ny)  is  zero  if  and  only  if  the 
mean  value  (67)  differs  from  zero  for  at  least  one  y^v)(^)  (^  "  •••»  P)  since 

then  we  have,  correspondingly,  iv  =  0.  In  the  special  normal  form  (see  the 
definition  in  Section  4),  we  will  have  ii  =  0  according  to  eq.(74)  and  thus, 
accordLig  to  eqs.(150)  and  (152),  nb  =  0*  ■  In  the  other  cases  (ix  >  0,  or  ab¬ 
sence  of  iv),  we  have  nfe  >  0  [see  eq. (153)3 ♦ 

It  is  useful  to  give  an  e3q)licit  computation  of  the  transformation  of  ft 
into  the  Jordan  normal  form,  using  a  siJiple  example  in  which  all  possibilities 
occur.  Let 

n=2o,  i2“5,  (155) 

The  general  idea  is  demonstrated  in  the  following  matrices: 
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Section  7.  Solutions  of  the  Ad.joint  Homogeneous  Differential  /II 5 

F.r|uat1on  of  the  n^^  Order 

The  solution  matrix  f^Ct)  of  eq.(6)  conjugate  to  in  eq.(lZ4.9),  can  be 
calculated  in  accordance  with  eqs.(124)  and  (IZ+S).  Since,  however,  the  condi¬ 
tions  for  the  resonance  subcase  and  for  the  exceptional  subcase  [see  eq.(28)J 
contain  the  solutions  ^(t)  of  the  adjoint  homogeneous  system  (ll)  periodic 
with  P  [respectively  of  the  adjoint  homogeneous  equation  (9)],  we  will  discuss 
this  first.  Analogous  to  eq.(l2/4)  and  because  of  eq.(l3),  we  have 

g't)  -g(t)  ^156) 

We  then  again  decoirpose  the  transformation  C  into  the  j  subtransformations  5^ , 

according  to  eq.(l48):  . 

^  (157) 

after  which,  we  will  investigate  the  individual  matrices  ^'^(t)  for  u  0,  1, 

j,  whose  first  and  last  matrices  are 

^  resp.  ^  ’  o  ( 1  58  ) 

We  then  partition  the  matrix  b(t;,  in  accordance  with  the  structure  of  ft  in 
eq,(126),  into  s  2  submatrices 

(159) 

LX  D 

where  3  ,  in  addition  to  zero  columns,  contains  only  the  first  j  -  u  columns 
of  ^3»  while  ^70  contains  the  next  ^1%  columns,  ^3i  next  following  ^mi 
colxmins,  and  so  on. 

It  is  useful  to  introduce  the  following  notation: 

t = ‘'"t- ,/•  (160) 

In  the  last  step  (m-  =  j )  we  will  omit  the  superscript  [see  eq.(l50)j: 

^ •••  ^  °  ^  J.  ' *-( f)»  -  (161) 
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Then,  the  following  theoran  applies:  /116 

Theorem  1?  The  solutions  £tvi  (t)  (v  =  1,  2,  s)  of  eq.(52)  periodic 

with  P  are,  in  the  same  sequence,  identical  with  the  solutions  Zjtvj)  (t)  of 
eq.(9),  again  for  v  =1,  2,  s.  Consequently, 

for  y-  1,2 . S  .  (162) 

Proof:  In  the  matrix  °3(^)  according  to  eq.(l57),  the  solution  vectors 
periodic  with  P  and  resulting  from  the  periodic  ztvi  (t)  for  v  =  1,  s,  are 

exactly  the  jjtv]  (t)  which  means  that  they  are  standing  exactly  in  the  last 
colvmm  of  the  °3v(t)  for  v  =  1,  2,  ...,  s.  A  successive  application,  on  the 
right-hand  side,  of  the  transformations  described  in  the  preceding  Section 

)^ ,  =2^  and  )  (if  they  occur  at  all)  to  the  matrix 

will  transform  the  last  column  of  each  *^~^v(t)  unchanged  into  the  last  column 

of  ^3v(t)(v  =  1,  2,  s).  Consequently,  also  the  last  component  of  these 

colvunns  will  remain  unaltered.  5y  this  inductive  syllogism,  the  proof  for 

eq.(l62)  is  obtained. 

In  the  case  nfe  >0,  i.e.,  ii  >  0  or  nonexistent,  a  further  solution  of 
eq.(9)  periodic  with  P  will  occur  in  accordance  with  eq.(l54)  and  theorem  5. 

For  this,  the  following  is  valid: 

Theorem  IL:  In  the  case  ii\j  >  0,  the  auxiliary  solution  vector  periodic 
with  P  is  present  in  the  last  column  of  ^3o(t);  its  last  conponent  which,  con¬ 
sequently,  must  be  denoted  by  Z[toa]  (t),  has  the  following  form: 


i  CO 

2  --2  ^ 
Ml  1 


fti  ■  y 


^ 

'-St„ 

y*  i 


(163) 


If  no  iy  is  present,  i.e.,  for  v  >  X,  the  last  sum  naturally  is  omitted.  /117 

Proof:  First,  it  should  be  remembered  that  the  quantities  occurring 

in  eq.(l63),  originate  from  the  last  (j^")  rows  of  the  matrices  i  in  eq.(105). 


102 


whereas  the  quantities  t^^Ct)  originate  in  the  last  [(n  -  j)^^]  row  of  the  ma¬ 
trix  [see  eq.(60)] 

i-i. -llT  > 

where  the  row  index  n  -  j  has  been  omitted  here.  Since  all  these  functions  are 
periodic  with  P,  this  statement  is  valid  also  for  the  expression  (I63). 

A  coltunn  of  '*30  (t)  =3o(t)  occurs  for  the  first  time  as  the  only  column 
of  ^3o(t).  Since  this  column,  at  increasing  ij.,  always  remains  the  last  column 
of  ^3o(t)  with  M.  ^  1,  it  follows  -  as  in  the  proof  of  theoron  13  -  that  tids 
col\imn  cannot  change  anymore.  Consequently,  we  have  only  a  single  column 
vector  ^3o(t)  of  ^3o(t)  which  is  located  in  the  column  of  ^3(t);  its  last 
component  is  the  sought  Ztton  (t). 

In  the  case  ii  >  0,  formation  of  ^3(t)  =3(t)  *  /  will  directly 


yield 


,1  i  ~ 


^2.(n  - 


5? 


z  (t>  ■ 


(165) 


It  shovild  be  recalled  that  the  last  row  index  n  or  n  -  j  had  been  omitted  in  the 
elements  of  3  or  3.  In  the  case  that  no  ij  exists  (ja  >  X),  the  second  sum  is 
eliminated. 

According  to  eq.(l3)  and  in  analogy  with  eq.(64),  we  have 


T9'V 


<t)~ 


r  j»t 

e 

T"-  ’ 

-t 

e  (j  n? 

(T  7 

% 

> 

.  i  J 

with 


(166) 

/lie 

(167) 


as  is  readily  verified  from  the  relation 
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Taking  the  definition  of  ?)*(t)  in  eq.(6u)  into  consideration  and  also  consider¬ 
ing  that  nZjCt)  is  the  element  in  the  right-hand  lower  comer  of  3*^(t),  the 
following  is  obtained  from  eq.(l67): 


.1  •“  K "‘•-ZL  V 

7»i  »*-■> 

where  [see  eq.(64)] 

is  valid.  According  to  eq.(l65)  this  yields  the  intermediate  formula 


(168) 

(169) 


«> 


(170) 


Denoting,  by  C(v)+i^_i»  a  mv-con^onent  vector  which,  in  addition  to  zeros, 
contains  only  an  i  in  the  ((v)  +  iy  -  1)^^  conponent,  the  subsum  of  eq.(170) 
can,  be  written  as  follows  for  v  <  X ; 


i;  •  - >  !  «>» 


f  A  1 

*  - 

L  > 

Taking  into  consideration  the  formula  for  the  reduced  differential  equation, 
which  is  analogous  to  eq.(A3)  and,  in  its  subdivision,  analogous  to  eq.(65)  as 
well  as  considering  eqs.(103)  and  (lOZ,),  we  will  obtain 


*• 

^  - 
a 

j'^<w 

w 

e 

t 

1 

f 

;  M 

Z112 


(171) 
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where  the  1  as  the  ((v)  +  iv  -  1)^^  conponent  of  the  last  vector  originates 


from  C(  V)  *  iv-x •  This  last  vector,  however,  is  exactly  e 
that  we  obtain  further 


^ (  v) l 


t*j 


f  S  i 

jTti» 


J 


+ 


•''>t  ly-1 


i  M  “  for  t’T  •*,•••,4 


(172) 


In  the  case  v=X+l,  ...,  s,  taking  eq.(94)  into  consideration,  the  last 
vector  in  eq.(17l)  is  omitted.  This  will  yield 

^  for  »'-A.  .  (173) 

From  eqs.(172)  and  (173)»  we  directly  obtain  eq.(l63)» 

In  passing,  we  note  the  following: 

Theorem  15:  The  function  zcto]j(t),  investigated  in  theorem  1/+,  satisfies 
the  inhomogeneous  reduced  adjoint  differential  equation 

u  t-tv^.2.1.  (174) 

✓  ^  “■  ■ 

Proof In  view  of  eqs.(52),  (56),  and  (57),  it  merel/  must  be  demonstrated 
that  a  solution  vector  u 

u  .  -  -  •^(t) 


with 


r  N 


(t) 


(175) 

Zi2g 


(176) 


exists  whose  last  [(n  -  j)^*^]  component  is  Z[ioi3.(t)» 

The  variational  method  for  the  constants  [see  (Bibl.l),  footnote  9]  fur¬ 
nishes  the  following  relation  as  the  general  solution  of  eq.(175): 


Su  (t)  -  <J(t)  ih'  tT>  -^(riUT  i'C-j 


(177) 


vn.th  an  arbitrary  constant  vector  c  and,  as  the  last  component. 


(178) 


_ s  lyJ  ^  f  t  ^  V 

.(t)  -z:  w  + w  • 

^-3  ^,1  />■-*'>* 

According  to  the  general  theorons  [see,  for  example  (Bibl.l),  eqs.(l6)  and  (17) 
and  the  context  there],  eq.(175)  has  a  solution  vector  periodic  with  P  if  the 
following  applies  for  all  solutions  of  eq.(55),  periodic  with  P: 

—  /  1  r  =  o  ioT  y  t  • 

This  is  exactly  the  condition  for  >  0  [see  for  example  the  remiarks  after 
eq.(l54)]»  Consequently,  eq«(178)  can  then  be  determiined  as  a  function  periodic 
with  P.  The  lower  limit  in  the  integrals  can  be  so  defined  that  Cp_  =  1  applies 

for  u  “  (v)  +  iv  -1  ('^  =  1,  ••*,  ^)  while  it  has  a  value  of  zero  everywhere 
else.  Now,  eq.(178)  obviously  coincides  with  eq.(170)  which  means  that  theorem 
15  is  proved. 

By  j  differentiations  of  eq.(174)»  we  again  verify  that  ztco33(t),  obtained 
from  eq.(165)  or  eq.(170),  is  the  solution  of  eq.(9). 

Section  8.  The  Minimal  Order  of  Magnitudes  of  the  Solutions  and  /121 

their  Derivatives  for  the  Resonance  Case 

With  respect  to  the  resonance  case  for  the  differential  equation  (l),  the 
following  statements  can  be  made:  The  adjoint  homogeneous  differential  equa¬ 
tion  (9),  as  defined  in  the  preceding  Section,  has  the  solutions  zttojj  »  ^tiin  * 
•••>  periodic  with  P,  where  Ztnjj  ,  ...,  Ztjpjj  according  to  theorem  13 

is  identical  with  the  solutions  znj  ,  ...,  StoJ  periodic  with  P  of  the  adjoint 
homogeneous  reduced  differential  equation  (52)  and  zt[033  occurs  only  in  the 
case  of  mfe  >0.  Consequently,  if  the  resonance  subcase  exists  for  the  inhomo¬ 
geneous  reduced  differential  equation  (50),  for  an  index  v  >  0,  i.e.,  if  the 
following  is  valid  [see  eq.(28)]; 
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(179) 


/  'f-j  »  •  •  r.'j  ^  ®  ' 

then  the  resonance  subcase  also  will  exist  for  the  inhomogeneous  differential 
equation  (l)  for  the  same  index  v;  consequently,  we  then  have 

/  ^rc^4  ^ ^  “  Vj  ^  ®  •  ( 180 ) 

Speaking  generally,  we  will  denote  any  index  v,  i.e.,  also  v  =  0,  as  a  resonance 
index  provided  that  the  resonance  subcase  exists  for  this  v.  Thus,  we  have  the 
following:  If  v  is  a  resonance  index  of  the  reduced  differential  equation  (50), 
then  V  will  also  be  the  resonance  index  for  the  differential  equation  (1). 

We  must  now  differentiate  between  the  following  cases: 

I.  No  resonance  index  exists. 

II.  V  =  0  is  no  resonance  index,  but  at  least  one  resonance  index 
V  >  0  exists. 

III.  V  =  0  is  the  only  resonance  index. 

IV.  V  =  0  is  a  resonance  index,  and  at  least  one  index  v  >  0  exists 
which  also  is  a  resonance  index. 

The  points  I  and  II  also  contain  cases  in  which  the  solution  zcco33  is  not 
present,  i.e.,  in  which  ii  =  0  applies. 

In  the  Case  I.  either  the  principal  case  or  the  exceptional  case  is  /122 
involved  for  the  reduced  differential  equation  (50).  For  the  differential  equa¬ 
tion  (l),  at  ii  =0,  the  same  case  as  fo^*  the  reduced  differential  equation  (50) 
is  involved.  At  ii  >  0,  the  exceptional  case  is  present  for  eq.(l)  since  at 
least  one  solution  zttoia  »  periodic  with  P,  of  the  adjoint  homogeneous  differ^ 
ential  equation  (9)  exists. 

In  the  Case  II.  the  minimal  order  of  the  power  increment  m  of  x(t),  accord¬ 
ing  to  theorem  2,  is  determined  by 


(181) 


m  •  Max  (b^)  , 

ax.}) 

while  the  analogous  order  of  increment  m  of  the  solution  x(t)  of  eq.(l)  is  de¬ 
termined  by 


m 


-  Uax 


(182) 


with  mv  according  to  eqs.(l50)ff ,  with  the  resonance  indices  being  the  same  in 
both  cases. 

In  the  Case  III,  the  resonance  case  is  not  present  for  the  reduced  differ¬ 
ential  equation  (50)  but  is  present  for  the  differential  equation  (1).  The 
resonance  order  for  x(t)  is  defined  as  [see  eqs.(l50)  and  (152)] 

®  -  Min  (j,  i^)  ,  (183) 

where,  if  nc  ii  occurs,  eq.(l‘^l)  must  be  taken  into  consideration. 

In  the  Case  IV.  the  resonance  case  exists  for  the  reduced  differential 
equation  (50)  as  well  as  for  the  differential  equation  (l).  The  resonance  order 
for  x(t)  is  defined  by 

B  .  Max  (m^)  •  J 

In  determining  the  maximum,  the  quantity  nb  can  be  disregarded  if,  for  at  /123 

least  one  resonance  index  from  the  interval  1  <  v  <  S  one  iv  exists;  in  that 
case,  the  corresponding  mv  >  nb  according  to  eq.(l50)  because  of  mv  irii  > 

>  nn(j,  ii );  see  also  the  definition  at  the  beginning  of  Section 

Similar  statements  apply  also  to  the  derivatives  x' (t),  x"(t),  ..., 
x^'*~^^(t)  if,  in  the  cases  I  to  IV,  the  index  j  in  the  derivative  (t) 
in  eq.(l50)  for  the  mv  is  replaced  by  (.]  -  k).  The  given  increment  orders,  in 
the  case  of  resonance,  are  always  the  minimal  orders.  It  is  entirely  possible 
that,  for  exairple,  x(t)  represents  a  solution  with  minimal  order  (with  the  in¬ 
dex  j),  while  the  correspoi;ding  derivative  x*  (t)  has  a  higher  order  than  the 
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minimal  order  valid  for  j  -  1.  Conversely,  it  could  be  that  x'  (t)  is  a  solu¬ 
tion  with  the  minimal  order  valid  for  j  -  1,  while  the  once  integrated  fur.ctior 
x(t)  =  r  x’  (t)dt  has  an  order  which  is  higher  than  the  minimal  order  valid 
for  j.  This  is  due  to  the  fact  that,  for  different  indices  j  -  k,  no  differirig 
determination  of  the  parameter  constants  is  necessary  if  a  solution  with  minimal 
order  is  to  be  obtained  in  each  particular  case. 

For  explaining  our  discussions,  the  following  example  is  used: 

Let  mi  =  3>  fts  =  5,  =  8,  =  L,  nfe  =10 

ii  =  2,i3  =  3t  X-i  -  5,  and  is  do  not  exist. 


"*  If  the  definition  (195 )  is  used,  ea.(185)  can  be  made  more  rigorous  by  re¬ 
placing  the  index  X  by  Y . 
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This  results  in  the  Table: 


0-0 

1 

2 

5 

4 

5 

— 

6 

7 

8 

m 

0 

1 

2 

2 

2 

2 

2 

2 

2 

2 

0 

“1 

3 

3 

3 

4 

4 

4 

4 

4 

4 

■2 

5 

5 

5 

3 

6 

7 

7 

7 

7 

7 

*3 
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Of  these  numbers,  according  to  eq«(184),  the  maximum  must  be  formed  at  fixed  j, 
for  the  indices  for  which  the  resonance  subcase  exists.  /125 

For  exanple,  if  the  resonance  subcase  exists  for  u  =  3,  then  m  will  be  at 
least  equal  to  n^.  In  any  case,  however,  1%  is  the  maximum  for  j  ^  6  so  that 
m  =  nfe  =  j  +  3  applies  for  j  ^  6.  This  means  that  the  minimal  order  m  finally 
increases  linearly  with  increasing  j. 

Conversely,  if  the  exceptional  subcase  exists  for  v  =  3,  the  quantity  m 
will  never  be  larger  than  nfe  =  10,  If  the  exceptional  subcase  also  exists  for 
V  =  5,  the  quantity  m  will  never  exceed  nfe  =  7;  and  so  on.  It  is  easy  to  demon¬ 
strate  that  m  retains  a  constant  value  if  the  exceptional  case  exists  for  v  =  3, 
at  least  beginning  with  a  certain  index  j. 

In  general,  it  is  easy  to  confirm  the  following  theorem  on  the  basis  of 
eq.(185): 

Theoron  16:  If  the  resonance  subcase  is  present  for  the  index  v  with  the 
greatest  existing  iy,  i.e.,  for  the  minimal  order  m,  beginning  with  a 

certain  index  j,  will  increase  linearly  with  j  (see  also  the  definition  in 
Sect, 4).  Naturally,  this  index  X  is  already  defined  before  the  special  normal 
form  is  established.  Conversely,  if  the  exceptional  case  is  present  for  the 


110 


index  v  =  X,  the  minimal  order  m  vdll  ranain  constant  in  any  case,  beginning 
with  a  certain  index  j.  In  all  other  respects,  the  behavior  of  the  minimal 
order  m,  as  a  function  of  j,  follows  from  the  formula  (I84);  m  =  m(j)  is  parti¬ 
ally  piecewise  constant  and  partially  increases  linearly  with  the  slope  one. 

For  each  index  v  for  whichaniy  exists,  the  pattern  of  a  broken  curve  will 
be  obtained  according  to  eq.(185)  for  the  minimal  order  my  as  a  function  of  j, 
which  has  a  horizontal  slope  for  0  ^  j  <  iy  and  for  iy+i  ^  j,  whereas  it  has  the 
slope  1  for  iy  ^  j  ^  iy+ 1 . 

Par  o  ^  J  /  iy  pe  have  * 

<  for  ^  3  re  have  m,  m  By  *  > 

for  d  ^  pe  have  m,  -  Sy  +  ( j  -  iy)  , 

where,  if  iv  ?  no  longer  exists  (v  =  ^),  the  last  law  is  valid  for  all  j  ^  iy : 

/I26 


Since  the  following  relation  always  is  in  question  for  two  indices  vg  >  Vj  with 

existing  iy_  and  iy,  [see  eq.(74)] 

<  1 

“y2  ”  X  ^^2  " 

the  cuiTfe  for  the  index  Vs  will  begin  with  a  constant  which  is  greater  than  the 
end  constant  for  the  index  Vj ,  i.e.,  the  curves  for  the  various  indices  and 
V2  with  existing  iy  do  not  intersect.  For  each  index  v  for  which  no  iy  exists, 
my  =  my  will  be  constant  for  0  ^  j. 


Ill 


Thus,  the  following  statement  is  obtained  for  the  minimal  order  according 
to  eq.(184):  Let  v  be  the  highest  index  v  with  existing  iv  for  which  the  reso¬ 
nance  case  is  present.  In  addition,  let  v  be  the  highest  resonance  index  at 
nonexisting  iv.  Then,  the  following  obviously  applies  (see  definition  3  in 


Sect.L)! 


t, --J  fy 


(my)  m  Uax  (m  ,  m.), 

y  * 


(186) 


where  the  pattern  of  my  as  a  fiinction  of  j,  jrields  a  broken  curve  (see  Fig.l), 


while  m^  is  a  constant. 

We  wotild  like  to  mentiorTa  few  interesting  relations:  In  the  case  that  /127 
only  V  =  0  is  the  resonance  index  of  eq.(l),  i.e.,  in  the  case  III,  at  least 
one  solution  x(t)  periodic  with  P  of  eq.(50)  exists  as  we  already  know,  whose 
mean  value  can  be  either  equal  to  zero  or  different  from  zero.  Then,  the  fol¬ 
lowing  theorem  applies: 

Theorem  17:  If  ii  =  0,  i.e«,  if  no  solution  Ztcojj  periodic  with  P  of  the 
adjoint  homogeneous  differential  equation  (9)  exists,  then  -  if  at  all  -  solu¬ 
tions  x(t)  of  eq.(50)  periodic  with  P  will  exist  whose  mean  value  differs  from 
zero,  as  well  as  solutions  whose  mean  value  is  equal  to  zero. 

If  ii  >  0  or  if  no  ii  exists  and  if  a  solution  x(t)  periodic  with  P  whose 
mean  value  differs  from  zero  is  present,  then  also  every  other  solution  x(t) 
periodic  with  P  will  have  the  same  mean  value  differing  from  zero  and  v  =  0  will 
be  the  resonance  index.  Conversely,  if  a  solution  x(t)  of  zero  mean  value  and 
periodic  with  P  exists,  then  also  all  other  solutions  x(t)  periodic  with  P  will 
have  the  mean  value  zero  and  v  =  0  will  be  the  exceptional  index. 

Proof:  We  can  write  the  general  solution  x(t)  of  eq.(50),  periodic  with  P, 
in  the  following  form.:  ^ 

X(t)  .  x*(t)  ?(/)(<:)  .  (187) 
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xhere  x  (t)  is  a  special  particvilar  solution  of  eq.(50),  periodic  vdth  P,  while 
the  sum  next  to  it  represents  the  general  solution,  periodic  with  P,  for  the 
homogeneous  reduced  differential  eqioation  (51)* 

In  the  case  ii  =  0,  the  quantity  y^^^  will  be  a  solution,  periodic  with  P, 
of  the  homogeneous  differential  equation  (5O)  with  the  mean  value  1.  The  other 
solutions  ~  2,  p)  of  the  homogeneous  differential  equation  (5O) 

have  the  mean  value  0.  This  shows  directly  that,  by  a  suitable  selection  of 
the  constamts  Ci,  the  solution  5c(t)  can  be  made  to  have  a  mean  value  of  zero. 
Conversely,  it  is  also  possible  to  make  certain  that  x(t)  has  a  mean  value  /I 28 
differing  from  zero. 

For  ii  >  0  or  for  nonexistent  ii ,  all  solutions  y^ ,  periodic  with  P,  of 
the  reduced  homogeneous  differential  equation  (51)  will  have  the  mean  value 
zero.  From  this  it  follows  that  all  solutions  x(t),  periodic  with  P,  must  have 
the  same  mean  value  as  x^(t)  (the  case  $  *  0  is  included  here). 

Then  it  merely  remains  to  be  demonstrated  that 

p  ^  ®  1  resonance 

from  r  X(t)dt  it  follows  ^  r  ^  ~  0  ^exceptional*  (188) 

i  [resp.  =0J 

We  will  demonstrate  this  indirectly:  If,  for  the  differential  equation  (l), 
the  exceptional  case  were  present,  i.e.,  if  v  =  0  would  be  the  exceptional  index 
[the  principal  subcase  cannot  be  present  for  v  =  0  since  the  homogeneous  differ¬ 
ential  equation  (6)  has  at  least  the  solution  y(t)  =  1,  periodic  with  P],  a 
solution  x(t)  of  eq.(l)  periodic  with  P  would  exist  whose  derivative 
5c^'*^(t)  =  x(t)  is  a  not  identically  vanishing  solution,  periodic  with  P,  of  the 
reduced  differential  equation  (5O)  with  the  mean  value  zero.  However,  this 
woiild  mean  that  all  solutions  x(t),  periodic  with  P,  necessarily  must  have  the 

p 

mean  value  zero.  If,  consequently,  [*5c(t)dt  f  0,  it  follows  necessarily  that 

'6 
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V 


=  0  is  the  resonance  index.  If,  conversely,  f  x(t)dt  =  0,  a  solution  x(t)  of 

o 

the  differential  equation  (l)  can  be  obtained  (see  the  auxiliary  theorem  in 
Sect. 3)  by  j  integrations  of  the  fxmction  x(t),  which  is  periodic  with  P  and  has 
a  mean  value  zero.  This  mecins  that  the  resonance  case  cannot  exist  for  the 
differential  equation  (l)  so  that  v  =  0  must  be  the  exceptional  index  since  the 
main  subcase  had  already  been  excluded  above  for  v  =  0.  This  proves  theorem  1?. 
From  this  theorem,  the  following  statements  are  obtained:  Whenever  the  /I 29 

p 

qxoantity  J  znon  (t)  f(t)dt  differs  from  0,  i.e.,  whenever  v  =  0  is  the  reso- 

0  P 

nance  index,  f  x(t)dt  will  differ  from  zero  for  a  possibly  existing  periodic 

0  P 

solution  x(t)  of  eq.(50).  Whenever  J  zttoij  (t)  f(t)dt  is  equal  to  zero,  i.e., 

0  P 

whenever  v  =  0  is  the  exceptional  index,  the  quantity  T  x(t)dt  will  be  equal  to 

0 

zero  for  a  possibly  existing  periodic  solution  x(t)  of  eq.(50)»  Consequently, 
a  relation  must  exist  between  the  mean  value  of  such  an  x(t)  and  the  mean  value 
of  the  function  2((on  (t)f(t).  Thus,  the  following  theorem  applies: 

Theorem  18:  If,  in  the  case  ii  >  0  or  nonexistent  ii ,  the  reduced  differ¬ 
ential  equation  (50)  possesses  a  solution  x(t)  periodic  with  P,  the  following 
equality  will  apply: 


r  7 


(189) 


Proof:  This  can  be  proved  by  means  of  Lagrange’s  identity  [see,  for  ex¬ 
ample  (Bibl.4),  Sect. 5.3] 


•  (>  Vqj) 


(ao  -  1)  in  the  following  manner: 
Using  [see  eqs.(50)  and  (174)3 
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Lt?)  -  f(t)  and 


we  obtain 


f(t)  -  x(t)  -  at  L 


(190) 


where  L*[x,  Ztton  3»  because  of  the  periodicity  x(t)  and  Zctoij  (t)*  is  a  func¬ 
tion  periodic  with  P  whose  derivative,  obviously,  has  the  mean  value  zero.  From 
this,  the  argument  is  directly  obtained. 


Section  9*  The  Order  of  Magnitudes  of  the  Derivatives  of  the 
Solutions  in  the  Resonance  Case 


/13P 


In  statements  on  the  order  of  magnitude  of  the  derivatives  of  the  solu¬ 
tions  x(t)  of  eq.(l),  for  the  resonance  case,  it  is  of  iirportance  whether  the 
solutions  y((V))(i^  “  ^((v))(t)  [see  eq.(l4)]  of  eq.(6),  periodic  with  P,  are 
constant  for  v  =  0,  1,  ...,  $  or  whether  they  actually  depend  on  t.  Here, 
eq.(2l)  and  thus  also  eq«(l6l)  are  analogously  defined: 


A  ^ 


((>'))  -  ■o+*i'*’*2****''‘*,»l‘*’  1  ■  ■i'*^“2'*’*****y-.T  ♦ 
tUin  (d,iy)  +  1  -  (y)  ♦  Min  (d,iy). 


(191) 


For  this  reason,  let  us  start  with  an  auxiliary  consideration  which  describes 
the  transition  of  the  matrix  $(t)  from  eq.(108)  to  the  matrix  $°(t)  from 
eq.(l2/^).  Primarily,  according  to  eqs.(12/j,),  (14)»  (92),  (113),  and  (123),  we 
have 

e  ^  a  Ui  JL  s  ^  ft)  ^  -X" a 

=  inTr'’ 

1*  f 

^'(t)  s  ffl)  .jr  (192) 

with  5  from  eq.(lij,8).  In  the  elonents  of  the  column  blocks  '^^(t)  with  v  =  1, 
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I 


X  of  ^(t)  from  eq.(108),  the  partition  is  made  analogously  to  eq.(lOl), 

Then,  we  can  write  ^ 

^  t  »  ( X93 ) 

where  the  matrix  "^(t)  contains  only  zeros  in  the  columns  1,  j,  and  contains 

ic*^'’i(t)  of  zero  mean  value  in  the  columns  j  1,  ....  n;  the  matrix  £  contains 
only  ones,  etc.,  in  the  first  j  places  of  the  main  diagonal  and,  on  partitioning 
(in  a  readily  understandable  symbolism)  /I3I 

1  -  /  “  t  ,  (m) 

1 

into  £ ,  it  will  also  contain  ones  in  the  diagonal  below  -45°  beginning  at  the 
element  1,  (  v)  +  i  v  *  this  concept,  it  is  necessary  to  consider  not  only 

eq.(l05)  but  also  the  fact  that  the  rows  of  5i(t)  result  from  each  other  by  suc¬ 
cessive  differentiation  [see  eq.(8)].  An  application  of  the  transformation  ma- 
trix  ®  to  $(t)  furnishes,  in  all  columns,  always  only  elements  of  zero  mean 
value.  For  this  reason,  the  effect  of  the  individual  partial  transformations 
(148)  and  (147)  on  £  will  be  investigated  first. 

For  greater  clarity,  we  will  do  this  on  the  exanple  (155)  given  at  the  end 
of  Sect. 6,  with  the  exception  that  now  j  =  7  (instead  of  j  =  4)  is  used  so  that 

all  imaginable  cases  [see  eq.(185),  specifically  the  third  relation]  will  be 

* 

covered  by  this  exanple.  The  matrices  5^,  and  58^,  constructed  earlier, 

need  now  be  only  supplemented  by  7  -  4  “  3  ones,  placed  in  front  of  the  main 
diagonal.  The  following  19  matrices  perform  the  stepwise  transformation  of  £ 
into  £°  =£7  . 


116 


It  is  also  useful  to  perform  a  stepwise  transformation  of  the  matrix  $(t) 
into  »  (t)  on  hand  of  the  same  examples.  This  will  be  done  in  the  following  /I36 
11  matrices  where  (except  for  a  few  obvious  exceptions)  only  the  column  indices 
of  the  functions  i$v(t)  =  'Pv(t)  occurring  in  eq.(lOl)  are  given,  which  at  first 
appear  in  the  (j  +  l)^^  row  forl^Ct)  =  #o(t). 
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22 


L25 


Let  the  first  row  of  the  matrix  Qj  again  be  written  explicitly 

for  the  cases  j  =  1,  2,  U*  and  7,  i.e.,  y  “  0»  1»  2,  and  3?  ^  “  3;  ix  “  2;  /1L3 

ia  =  3;  ia  =  5: 

j-’|M  V'^WWK.*  ' 

•«„>,  '  %,jn '  *♦ 

tofi*  ?o*»l  tv**  T»+«*  ^fvis'^vt^"  tj)*/ 

'~  %tJt¥  1^'  ‘  ' 

'tv-tr  • 

Since,  in  the  general  case,  the  considerations  are  completely  analogous  to 
this  example,  we  can  formulate  the  following  theoron: 


Theorem  19;  If  j  is  located  in  the  interval 

/  IS  r,o  limit), 


(195) 


the  following  is  valid  [see  eq.(97)]5 


f  a>  = 


for  p's 


for  >/s  V+l,—  /f  , 

for  ye  f  ••'•*  ** 


(196) 


(where  the  last  row  is  of  no  interest);  in  addition,  we  have  the  following  /ILA 


in  the  case  j  /  iy 


if  muf  «r  •'  =¥.  .  .  ,^'Of  «hile  f  ■  fifth 

»a  *»»♦’  ***** 


(197) 
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where  the  function  cp 


occurs  only  in  the  case  of  j  -  iv  >1.  Here, 


( Y) ) +J-1 y-i 

we  have  ^  ^ 

^  const,  for  y-  .  (198) 

[This  latter  follows  from  eq.(95)  or  (97) •] 

Let  us  now  consider  an  arbitrary  normal  solution  x(t)  of  eq.(l),  which  we 
will  retain  in  nhat  follows,  and  let  us  make  statements  <5n  the  increment  of  its 
derivatives  x^'‘^(t)(k  =  0,  1,  n  -  l).  For  this  solution  x(t),  eqs.(19), 

(31),  (26),  and  (32)  are  valid.  Here,  it  must  be  considered  that,  in  ea.(3l)» 
the  indices  v  are  arranged  in  a  different  manner,  namely,  first  the  resonance 
indices,  theft  the  exceptional  indices,  and  finally  the  principal  case  indices. 

In  order  to  retain  our  above  notations,  we  will  replace  eq.(3l)  ly 

u)v  =  mv ,  if  V  is  the  resonance  index  • 

(199) 

(Uy  <  niy,  if  V  is  the  exceptional  index 


with 


Instead  of  eq.(29)»  we  must  then  write 

with  X  t  0. 


(200) 


Similarly,  eq.(19)  is  rewritten  as 

x(t)  »  J  ’'x(t). 


(201) 


The  order  of  increment  of  the  investigated  normal  solution  x(t)  in  eq.(32)  will 
be  denoted  by  ou  so  that,  according  to  eqs.(33)  and  (34)  but  using  our  new  nota¬ 
tion,  the  following  is  valid; 


w  «  Max  i  Max  Max  «  id  • 

(y=o,l,..  .,f )  (yRes)  (/Fes) 


(202) 


Independent  of  the  investigations,  made  at  the  beginning  of  this  Section, 


on  the  matrix  i°(t),  \.re  following  statement  can  be  made  directly  on  the  basis 
of  theorem  3»  /1^5 

Theorem  20;  For  a  definite  normal  solution  x(t)  of  eq.(l),  relative  to 
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1 


the  power  order  of  x(t),  x  (t),  x‘“"^^(t),  we  will  obtain  a 

broken  curve  of  the  following  type: 


In  the  case  =  0,  only  the  horizontal  at  the  height  u)  ronains.  [Corresponding 
patterns  are  obtained  for  all  '’x(t)  at  v  =  0,  1,  3  from  eq.(201)  and  the 

p  ertainin  g  de  ri vat i ves . ] 

On  the  basis  of  the  int\iitive  stat«nents  on  the  construction  of  the  matrix 
(t),  stricter  statements  can  be  made*  Let  j  be  located  in  the  interval  (195)* 
According  to  eqs.(23)  and  (25)»  we  have  the  following  for  a  resonance  index  or 
for  an  exceptional  index  v; 

’'0,(0  •  ^*  t  Const.  (  >' »  0,1,...  ,?  (203) 

[see  eqs.(30)  and  (31)] •  Consequently,  in  view  of  eq.(196),  we  have: 

Theorem  21:  For  an  index  v  s  3,  differing  from  y,  all  (t)  for  u  = 

=  0,  1,  *.*,  n  -  1,  have  the  same  power  order  (JUy  (see  Fig.3).  The  same  state- 


Powtr  order* 


1 


f  f  1  -  ■  • 

Fig.3 


v't  ^  ^ 
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ment  applies  also  to  the  index  v  =  y  in  the  case  of  j  =  iy •  This  theorem  /1l6 
has  its  necessary  con5)lement  in  the  following  theorem; 


Theorgn  22;  If,  in  the  case  j  /  iy  [see  eq.(195)J,  we  have  y  ^  o,  the 
power  orders  of  the  derivatives  x(t)  successively  decrease  by  1,  beginning 
with  (By ,  down  to  the  order 

ly  -  Min  (3-1^,  -V)  (20/,) 

after  which  they  remain  constant  and  equal  to  u)y  -  ( j  -  iy  )  in  the  case  of  uy  ^ 

^  j  -  iy  while  they  will  be  equal  to  0  in  the  case  of  uy  s  j  -  iy  (see  Fig./,). 


Proof;  In  eq.(199),  according  to  eqs.(23),  (24),  (25),  and  (30)  we  have  /l47 


r  m 


for  "  ■*  'y”  *  ' 


1:  i-. 


r 

•X  fX 


A-0 


■f 

r  tti'f  u>  • 


(205) 


Because  of  eq.(197). 

In  the  case  of  («y 
are  constant  for  M-  =  0, 


C “t )  .  Cons^v  for  ■  n  1  i » •  •  •  1 •  (20b) 

a' 

>  j  -  iy,  the  coefficients  9^,(1)  of  t  in  eq.(199) 

1,  •••,  j  -  (iy  1)  according  to  eq.(197),  with  the 
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(U 

coefficient  of  t  ^  being  different  from  zero  because  of  eq.(200).  Conversely, 

«»Y-(  J-lv) 

the  coefficient  of  t  is  not  constant.  An  application  of  theorem  3  to 

theorem  20  will  demonstrate  the  correctness  of  theoran  22  for  the  case  in  ques¬ 
tion  here. 

uu^-U 

In  the  case  uJy  ^  j  -  iy*  only  the  powers  t  with  U  =0,  1,  ...,  uoy  s 
^  j  -  iy  with  constant  coefficients  will  occur  in  eq.(199),  where  again  the  co- 

(U 

efficient  of  t  ^  differs  from  zero.  From  this,  the  statement  of  theorem  22 
follows  eilso  for  this  case. 

^  means  of  theorans  21  and  22,  the  theoran  20  can  be  made  somewhat  more 
rigorous.  Let  xis  investigate  a  certain  normal  solution  x(t)  of  eq.(l).  Since, 
according  to  eq.(A8)»  we  have  aB-j(t)  f  0,  all  derivatives  x^'‘^  (t)  with  k  =  j, 
j  +  l,  ...,  n-1  have  the  same  order  of  magnitude  in  accordance  with  theoran  4 
(see  also  the  ranark  made  there).  Consequently,  only  the  power  orders  for  k  * 

=  0,  1,  ...,  j  ranain  to  be  discussed.  Here,  several  cases  must  be  differenti¬ 


ated. 


If  we  have  the  following  in  eq.(202)  [see  also  eq.( 195)3 

w  > 


(207) 


or 


tj  M  simultmtously  0  ~  ly  t  (208) 

then  the  highest  coefficient  Yo(t)  in  eq.(32)  is  not  constant  emd  all  deriva¬ 
tives  x^'‘^(t)(k  =  0,  1,  ...,  n  -  1)  have  the  same  power  order  (see  theoran  21). 
Consequently,  in  Fig. 2  we  have  t  =  0.  Figure  5  shows  the  power  orders  of  the  /lZi.8 
derivatives  of  x(t)  as  a  dotted  line  and  those  of  ^x(t)  as  a  dot-dash  line*. 


If 


3  f  iy. 


(209) 


*  All  quantities  ''x(t),  not  entered  in  the  illustrations  given  below,  yield 
horizontals  located  below  the  dotted  cui*ve. 
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then  three  additional  cases  must  be  differentiated 


a)  An  index  v  =  r  with  r  ^  0  and  r  f  y  exists,  for  which 


(ao) 


Also  in  this  case,  the  quantity  'l'o(t)  in  eq.(32)  is  not  constaint  so  that  we 
again  have  t  =  0  in  Fig.2.  Now,  Fig. 6,  in  which  ^ xit)  with  its  derivatives  is 


Fig.  5 


plotted  as  a  dashed  line,  illustrates  the  power  orders  of  the  derivatives. 

b)  If,  in  addition  to  eq.(209),  the  following  applies  for  all  indices  v  <  S 


differing  from  y 


USX  ^  w.  <  w  -  ^  <  w 

(y,: . >  r  r  r 


/1U9 

(211) 


Fig.  6 


then  all  ^5(1)  for  6  =  0,  1,  ...,  -ty  “  1  eq.(32)  are  constant,  whereas  'l'^^(t) 
is  not  constant.  Consequently,  we  have  in  Fig.2  t  =  iy  =  Min  (j  -  iy ,  ^y);  see 
also  eq.(20/4.). 
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c)  If  at  least  one  index  q  exists  so  that,  in  addition  to  eq.(209),  we  also 


Fig.8 

This  yields  the  pattern  shown  in  Fig.8". 


Section  10.  Construction  of  Solutions  x(t)  of  Eq.d’ 
Given  Derivative  has  the  Minimal  Order 


in  which  a 


iiiO 


We  will  now  investigate  how  high  the  power  order  of  a  singular  solution  of 
eq.(l)  must  be  in  order  that,  at  given  k  (k  =  1,  2,  . . . ,  n  -  l),  the  quantity 
x^''^(t)  has  a  minimal  power  order. 

For  this,  a  coii5)lement  to  theorem  19  is  required,  whose  correctness  can  be 
read  from  the  matrix  examples  of  theorem  19  in  the  same  manner  as  the  proof  of 

^  With  respect  to  b)  and  c)  see  theorem  2U  in  Sect. 10, 
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theoren  19  itself.  This  can  be  formulated  as  follows: 


Theorem  23:  In  the  matrix  *°(t),  we  have  [see  eqs.(l60)  and  (191)] 


i/ 

J 

■  '  •'»,<)  - 

Ot  > 

-(w 

'^■V* V *1*  '* l[v. V t 1- 

,  12) 

i/  ojii  ^  ^ 


i/  Vs  it* 


(213) 


(210 


(215) 


‘/  y*»  ^ A 


■> 


i/ 


»'>A 


(v, . Va)'  • 


fw 


ZL5i 

(a6) 


(a?) 


In  the  cases  (a5)»  (a6),  and  (217)»  we  have  [[v]]  =  [v].  For  producing 
the  genera  proof,  it  sho^d  be  stated  that  the  zeros  in  eqs.(a3)  and  (a^)  are 
due  to  the  permutation  matrices  while  the  negative  terms  in  eq.(2l3)  are  due 
to  the  superposition  matrices  18 .  The  matrices  have  no  effect  at  all  on 


In  the  case  v  =  y  =  0,  only  (cPdo))  ,  ...»  ^((o))^j-i  )  =  (li ,  Qb  , 

•  O4  )  remains  of  eq.(a4)  which  means  that,  in  the  case  j  =  1,  we  only  have 
(^Ko)))  =  (li  ).  The  situation  is  similar  for  eq.(a3). 
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the  matrices  ’'♦ui. 

First,  let  us  repeat  Fig.l  with  respect  to  the  minimal  order  of  the  deriva¬ 
tives  (t )  [see  eq.(19)],  with  k  as  the  abscissa. 


With  a  minor  modification  of  eq.(186),  we  can  then  define: 


9 


^*0 

y  4 


n%,  Md 


(: 


A, 

sjff*,.. 

>>  ■ 


') 


Z112 

(218) 


where  Y  is  defined  by  eq.(l95)«  Then,  the  relations  (185)  apply,  provided 
that  ^  is  substituted  by  Y  (see  footnote  on  p,109);  this  had  been  taken  into 
consideration  already  in  Fig. 9.  If  we  have 

.  -  Max  (219) 

then,  according  to  the  relations  (185)  modified  in  this  manner,  the  minimal 
orders  of  all  derivatives  x^'‘^(t)  (k  =  0,  1,  n  -  l)  of  a  normal  solution 

x(t)  of  eq.(l)  are  constant  and  equal  to  my  =  my.  This  constitutes  no  problem. 

Then,  only  the  case 

B  -  Max  (b^  ,  B^  )  -  B^  (220) 

remains  to  be  considered.  If  Y  is  the  resonance  index,  i.e.,  if  y  ~  the 
power  order  of  the  derivative  x^'‘^(t)  of  a  normal  solution  is  equal  to  the  power 
ordei*  my  of  ’^x^'‘^(t)  according  to  Fig. 9b.  Accordingly,  the  only  case  of  interest 
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is  that  in  which  y  represents  the  exceptioncd  index,  i.e., 

*  <*'•  (221) 

If  k  «  j  -  i=  +  1,  then  (t)  will  yield  the  minimal  power  order  of  a  normal 

solution  x(t)  in  accordance  with  Fig. 9a.  Information  on  the  case 

^  ~  ^  (222) 

is  obtained  ty  the  following  theorem: 

Theorem  2L:  If  a  given  normal  solution  x(t)  of  eq.(l)  is  resolved,  under 
the  assumptions  of  eqs.(195),  (220),  (221),  (222)  and  in  view  of  eq.(20l),  in 
the  following  form 


x(t)  -  x^(t)  ♦ 


with 


U) 


(223) 

i22U) 


[see  eq.(32)],  in  which  case  /I 53 

■"  '  S  *(»- '})  ■  ”‘5  ♦(>-  V,;  (225) 

[see  eq.(150)]  can  be  selected,  then  the  constants  occurring  in  the  ''x(t) 

of  eq.(22i+)  according  to  eq.(26)  can  be  so  modified  that  H'o(t),  Yi(t),  ..., 

Y^i(t)  are  constant  while  ’l'^(t)  is  not  constant,  with 

1  •  J  ■  i..  (226) 

On  substituting  this  modified  fimction  X5j.(t)  in  eq.(223),  the  quantity  x^v.(t) 
can  be  so  selected  that  the  derivatives  x^’‘^(t)  of  u)  =  uuy  from  eq.(225),  for 
k  =  0,  will  decrease  with  increasing  k  by  1  each  time  until  the  power  order 

M  >  Uaxfn  «  »  ®  y  )  •  (  227  ) 

is  reached.  In  the  case  i-i  =  m_,  the  decrease  of  the  power  order  proceeds  to 

V 

the  derivative  with  t  according  to  eq.(226)  while,  in  the  case  M  =  mv  =  mv 
it  proceeds  to  the  t  derivative  with 
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dices,  which  means  that,  for  these  indices,  the  constants  ^  normal 

solution  x(t)  can  be  arbitrarily  selected.  Consequently,  because  of  [see 


eq8»i7k)  and  (150)] 


t  M  9 

r 


•^5 


we  can  select  eq,{22u)  as  the  upper  limit  of  the  sum.  If,  analogously  to 
eq.(225),  v»e  assume  for  v  =  v  +  1,  v  +  2,  ...,  Y  -  1, 

■  *.  )  =  t(*  -  * 

and  take  into  consideration  [see  eq.(150)] 

«*/  *  •*»»•  _  •  ^  r* »  .  *  * )  * 

a  -  -  \  «'/ 

then  eq.(22/i)  can  be  written  in  the  form  of 


^  It) a  y  >  i - 

*  7^  V 


**  ^  1 
- -  (t)i 


(229) 


(230) 


(231) 


[see  eq.(199)  or  (31)3  with  [see  eqs.(26)  or  (2^)3 

J  Axtf' 

■'  ^  «•'»)+*  ' 


Zi^ 


(232) 


where  the  v*^v))  periodic  with  P  and  having  a  mean  value  zero  can  be  selected 
[see  (Bibl.l)  eq. (103)3.  A  comparison  of  the  t-powers  in  eqs.(22A)  and  (231), 


using  the  notations 


will  yield 


yr  y,  T  +  l,.  tr-1 


^  yxy 


(233) 


(23A) 


which  can  be  used  only  up  to  6  =  t,  in  which  case  any  summands  with  UvS  not  de¬ 
fined  in  eq.(233)  must  be  omitted.  It  must  be  noted  here  that  the  quantities 
(225)  and  (229)  are  the  highest  occurring  exponents  of  t  at  ^0o(t)  or  '’9o(t). 
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In  addition  to  eq«(206),  we  will  also  need  the  following  relation  which 
has  been  obtained  by  means  of  eq«(2l4)  from  eq.(232): 


i  f  k  f  m  i  »  •  •  •  > 


t,iS  \  >i 


(235) 


In  the  case  m=  =  0,  the  second  sum  of  the  second  formula  in  the  system  (232) 
will  also  occur  in  eq.(235)  for  6=1,  which  must  be  specifically  taken  into 
consideration  later  for  eqs«(237),  (24l)»  (243)»  and  (2UU)»  For  6  =  0,  1, 
j  -  iy  -  1,  eq#(206)  is  contained  in  eq*(235)  since  in  that  case,  according  to 
eq.(235),  the  quantity  ^  becomes  negative,  meaning  that  the  sum  does 

not  occur  in  eq.(235)«  This  fonniala  is  applicable  no  matter  whether  j  =  iy 
or  j  f  iy  ,  as  can  be  confirmed  by  means  of  eqs.(2l4)  and  (215)»  /156 

In  accordance  with  eq.(232),  we  obtain  for  the  remaining  summands  in 


eq.(234),  i.e.,  for 
by  means  of  eq.(2l3)> 


(236) 


=  ‘Hr”'*  t 


+  - '  *!/  T  “ 

2 - -  -4.  ) 


■T.  -W, 


*j 

-A  Y 


’A 


(237) 


Since,  according  to  eq.(233)»  we  have 


I  y  >  -  2  >  >“ 


(238) 


the  sum  does  not  occur  in  eq»{23u)  for  ^be  case  that 


f“  0,1,.. ^  » 


(239) 


so  that  here  '^'^(t)  according  to  eq.(235)  [resp.  eq.(206)]  is  constant^.  Conse¬ 


quently,  we  now  need  consider  only 


^  ■  j  —  t  0  “  1 


Next,  we  calculate,  in  accordance  with  eqs.(23it),  (235),  and  (237), 


(240) 


(t).  '^JL  i  «)  + 

♦  'H-  -•*  »  . 


-A  -a)  ‘Of 

•<=  ?  .  >».=•  \ 


'  /«•  -A  *  -A  /  'fv+OtA  '  i 

Cr-^^x.  =,  y  r«. 


(241) 


Then,  we  select  the  arbitrary  '’d^  such  that  the  parentheses  vanish,  i.e.. 


J.  -Jt  “  J-  /4.  »»  =  *y  *<.1  , 


(242) 


where  only  eq.(27)  for  v  =  v  must  be  taken  into  consideration.  Then,  the  only 
remaining  terms  of  eq.(24l)  are 


<M=  *1-  t  ^  ^ 


/»  -A 

%1,X 


(2/,3) 


Here,  the  sxim  occurs  only  if  [see  eq.(23l)J 


i.e.,  because  of  eq.(226),  if 


5-t,x  *  ^c; 


Consequently,  according  to  eq.(24l),  the  quantities  H'o  ,  V'l  ,  ...,  'V^~i  are  con¬ 


stant  while 


(244) 


is  not  constant  [see  eq.(l98)].  In  view  of  theorem  3  (or  theorem  20)  this 


At  j  “  iy  >  these  terms  do  not  occur. 
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proves  theoron  2U  under  the  additional  stipulation  that  x^(t)  in  eq.(223)  con¬ 
tains  no  t-powers  higher  than  t*’,  with  K  from  eq.(227)»  For  this,  it  must  first 
be  established  that  no  v  of  the  series  v  +  i,  X  can  be  a  resonance  index 

and  thus  need  not  furnish  any  t-power.  For  such  a  v,  according  to  eqs.(7/4.) 
and  (150),  we  would  have 


which  constitutes  a  contradiction  to  eq.(220).  For  the  resonance  indices  of  the 
series  X  +  1,  3,  according  to  the  definition  in  Section  U,  the  highest  in¬ 

dex  has  the  greatest  mv  =  my  which  means  that,  according  to  eq.(220),  we  must 
have  my  <  m=.  Here,  it  shoilLd  be  noted  that  the  case  m  =  m^  =  m^  in  eq.(220) 
had  been  taken  care  of  already  when  treating  eq.(2l9)»  This  will  immediately 
produce  the  bound  (227)  since  the  indices  v  >  p  furnish  orJy  periodic  components 
to  a  normal  solution  x(t)  of  eq.(l). 

As  an  application  to  theorem  2A»  we  can  use  the  example  on  p.l09  with  /1 58 
j  =  Y,  in  which  case  we  can  assume  =  0,  2,  A  as  resonance  indices.  We  here 
have 

7,  5* 

w  *  *  ? 

so  that,  according  to  eq.(225),  it  can  be  computed  that 

w.Wj-5+(8-5)-10 

and,  according  to  eq.(226) 


1  .  8  -  5  -  5. 

In  Fig. 12,  the  dotted  line  indicates  the  power  order  of  the  derivative 
of  the  solution  x(t)  according  to  theorem  2A,  while  the  solid  line  shows  the 
minimal  power  order  according  to  Fig. 9a;  the  dot-dash  line  gives  ^x(t).  This 
latter  line  generally  also  is  valid  for  the  solution  (223)  with  o'  from  eq.(225), 
i.e.,  for  ou  =  10,  and  can  drop  below  the  value  ms  =  7  only  at  a  special  selec¬ 
tion  of  the  constants  for  the  selection  (2A2),  this  can  drop  to  the  minimal 


order  of  ^  =  5 


So  as  to  obtain  an  example  also  for  the  case  (228),  our  above  statonoits 
are  modified  such  that  we  lower  nt  =  6  and  add  v  =  5  as  a  resonance  index  to 


Fig.l3 


V  =  0,  2,  and  k.  We  now  have  again 


“  -  -  10;  1  -  5, 


y  -  2,  m  «  7,  n  »  vhereas 

y  J  y  |»  7  ? 


Consequently,  we  have  in  eq.(227). 


U  «  Max  (5;6)  «  6 


and,  according  to  eq.(228), 


1^  =  5  -  (6  -  5)  =  U 


The  new  conditions  are  illustrated  in  Fig.13. 

Section  11.  Method  for  the  Formation  of  Examples  /I60 


Theorem  25:  The  functions 

for  (2/4,5) 

where  the  ®p_(t)  are  n- times  continuously  differentiable  and  periodic  with  the 
period  P,  form  a  fundamental  system  of  a  differential  equation  of  the  n^^  order 
having  coefficients  periodic  with  P,  provided  that  the  Wronski  determinant  of 
the  functions  yi ,  y^  differs  from  zero  for  each  value  of  t  within  the  in¬ 
terval  0  ^  t  ^  P.  Here,  the  matrix  515  ?)~^(0)?)(P)  contains  only  elementary 

conponents  of  the  order  1. 

Proof;  Obvioxisly,  the  functions  yg^(t)  satisfy  the  differential  equation 
of  the  n^^  order: 


*  e*  tn. 


V  /’■“  *  (*'•»  ♦  "it.) 

y  >e’  ft  f  f  ^  V  '*') 

la'”' «%“  % -t “rt  ^  f,Hv  C'ti  t-t  Cj 


(2/46) 


If  the  ejqjonertial  functions  are  canceled  from  the  columns,  a  differential  equa¬ 
tion  with  coefficients  periodic  with  P  is  obtained,  in  which  case  the  factor  of 
the  highest  derivative  which,  except  for  the  sign,  coincides  with  the  Wronski 
determinant  of  the  functions  yi ,  ...,  ,  according  to  definition  is  a  continu¬ 

ous  fimction  differing  from  zero.  In  addition,  the  following  is  valid  for  the 
fundamental  system: 


t  »• 


(2/4?) 
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-ft  ♦ 


/  •  •  •  #  , 


i  ^,9 

Tit~ft  '-  vV  ♦»tlt/«,  tt-tt  ^ 


^Cttr)*1^).|  * 


From  this,  it  follows  that  [see  (Bibl.l),  eq.(lO)] 


(247) 


/161 


(248) 


(249) 


Consequently,  the  matrix  ^  consists  of  elementary  components  of  the  order  1. 
Theorem  26;  The  fvinctions 


(M  »  e  ^  (t> 

•*%/**) 

=  /*  [  I?  Y,"*+  ^  t 


(250) 


where  the  cci(t),  ...,  <i^(t),  are  m- times  continuously  differentiable  fiinctions 
with  the  period  P,  form  a  fundamental  system  of  a  differential  equation  of  the 

X  'U 

m  order  having  coefficients  periodic  with  P,  provided  that  the  Wronski  detei^ 
minant  of  the  functions  yi(t),  y, (t)  differs  from  zero  for  each  value  of  t 

within  the  interval  0  ^  t  ^  P.  The  corresponding  matrix  $  consists  of  exactly 
one  elementary  component  of  the  order  m. 

Proof;  The  differential  equation  of  the  m^h  order,  which  is  ainalogous  to 
eq.(246),  after  reducing  by  e"°*  and  introducing  the  operator 

T^  -  d 

D  "  "Hr 


(251) 


as  well  as  the  notation 


•f  to  /or] 


(252) 


-  siunmands  in  which  the  index  cp  ^  0  must  be  replaced  by  0  -  can  be  written  in 
the  following  form; 

(253) 

. A  tLll 

The  first  row  is  obtained  directly,  while  the  other  rows  are  obtained  by  com¬ 
plete  induction  with  reference  to  the  row  index.  If  the  power  factors  are 
eliminated  by  forming  column  combinations,  a  differential  equation  having  co¬ 
efficients  periodic  with  P  will  be  obtained.  Again,  the  factor  of  the  highest 
derivative  y^"^(t),  except  for  the  sign,  is  equal  to  the  Wronski  determinant  of 
the  functions  yi(t),  ...,  y.Ct),  i.e.,  is  continuous  and  different  from  zero. 

In  addition,  analogous  to  eq.(2/i7),  we  have 


igto. 


3,  ‘h  ' 

.  3,‘  • 


»  3* 
•> 

•j  J’fM 
* 

■j  rf'.w 


(25A) 


-i,'  '-J.  *4t'  Jj  ' 


Z162 


'‘'  'f~.'i!i  -•?«*" 


m-| 


’^*  {»-*>■  ■*•*  ^ 


(255) 


•’}' 


p 


we 


i  T 


i-'  • 


i 


From  this,  again  in  accordance  with  another  paper  [(Bibl.l),  eq.(lO)],  it 
follows  that 


^=e 


4  P 


1  p  t: 

1  ^ 


1  r  c: 

t'. 


-*»; 

1 

Cm-»»f  ’ 


1  T  .. 


1 


«V 


«rp 


e 


(256) 


Consequently,  the  matrix  $  consists  of  a  single  elementary  conponent  of  the  /l6/i 
order  m.  Theorem  25  for  M<  =  1  represents  a  special  case  of  theoran  26  for  m  =  1. 
From  this,  we  directly  obtain  the  following  theorem: 

Theorem  27:  If,  in  forming  the  differential  equation  {2I4.6)  or  eq.(253) 
from  several  function  systems  of  the  type  of  eq.(250),  i.e.,  for  example, 

a  differential  equation  of  the  order 

n  ■  IQ^  +  TO^  +  •  •  • 

is  obtained,  with  continuous  coefficients  P,  and  one  matrix 
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(257) 


where  )Pv  is  of  the  order  mv  • 


Section  12. 


In  accordance  with  theorem  26,  with  m  =  2  and  cr  =  0,  we  start 


from  the  fundamental  system 


7/.  *  X  ■  sin  t ,  y  •  t  •  sin  t  -  5coat  +1. 


(258) 


These  functions  are  solutions  of  the  differential  equation  [see  eq.(253)  or  the 


formula  analogous  to  eq.(246)] 


t  ..2  Sint  cost  .  Sint  a. 

L  ay - ^ — 2 - y  * 

3  ♦  sin  t  -  cos  t 


(259) 


4  ♦  COS  t _ 

3  ♦  sin^t  -  cos  t 


y  •  o 


in  which  the  denominator  is  positive.  The  corresponding  solution  matrix  ?i(t) 


reads 


with 


r  f" 

^  ^  ^  ''“5  cos  tit 

^^(t)  -  |(t)  •  - -  - p_ 

cos  t  4  sin  t  0(1 


■[:;) 


(260) 


(261) 


The  transposed  reciprocal  matrix  3(t)  -  (^~^(t))^  will  then  be 


-*-r 
•  t 


0(0  -  2J-(0  •  e  ‘  ‘ - ^ - 

^  '  3+8in  t-cos  t 

4  sin  t  -  cos  t  1  0 

•  • 

—1+3  cost  Sint  —  t  1 


(262) 


—'1+3  cos 


The  elements  of  the  last  row 


^  -t.  Sint  -  1.3  cost  ^ 

Z/-X  ■  — - - -  and  *, 

'  '  3.8in  t  -  cos  t 


3+sin  t-cost 


(263) 
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form  a  fundamental  system  for  the  homogeneous  differential  equation  adjoint  /166 


to  eq#(259): 

a  |°st  *  Sint  *J'^  yoa  t  ( 26^ ) 

13  +  sin  t  -  coat  '  5*8in‘^t-co8t 

Next,  we  consider  the  solutions  x(t)  of  the  inhomogeneous  differential  equation 


2 

T  BS1  2  Sint  cost  ♦  Sint  4  ♦  cos  t 

Ltxja  X - -  X*  ♦  - 2 - 

3  *  sin  t  -  cost  3'«'Sin  t-cost 


with 


(265) 

(266) 


(267) 


f(t)  ■  3  atzLt  -  coeta 

Since  we  have  a  solution  Zfi]  ,  periodic  with  P,  we  car.  form  the  integral 

tr  tw 

I  f(t)dt  •  J  sin t dt  •  o 

cir.d  find  that,  for  the  inhomogeneous  differential  equation  (265)»  the  exception¬ 
al  case  is  present.  Then,  the  general  solution  x(t)  is  obtained  in  accordance 
with  the  variational  method  for  the  constants  [see  eq.(18)],  as  follows; 

x(t)  ■  -  t  Sint  ^  2  5  costt  -  cost)  + 

♦  t  •  ""cg  Sint  +  "’cp(1-3  cost)  ♦  sint. 


(268) 


with  the  derivative 


x'(t)  -  -  t  cost  -  sintt  +  t  •  "'c  cost  ♦ 

o 


1  i 

4.  Cp  •  4  sin  t  ♦  c-  cos  t  . 


If  ^  q,  =  1  is  selected,  then  x(t)  is  periodic  with  2P  at  arbitrary  ^  ci .  How¬ 
ever,  if  ^  Co  ^1,  the  power  order  of  x(t)  will  be  m  =  1,  Since  the  factor  of  t 
in  x(t)  is  a  nonconstant  function,  periodic  with  2rr,  also  the  derivative  x  (t) 
has  the  same  power  order  m*  =  m  =  1  (see  theorem  3,  resp.  theorem  20), 

We  then  turn  to  the  differential  equation  /167 


U*i) 
114m  ^ 


UfO  i 


I/) 
JT  * 


(269) 


at  j  ^  1  which  is  transformed,  with  ~  into  the  reduced  differential 

equation  (265).  First,  we  construct  a  fundamental  solution  matrix  ^{t),  accord¬ 
ing  to  the  method  described  in  Section  5,  for  the  differential  equation  homo¬ 
geneous  to  eq.(269).  By  selecting  successively  j  =  1,  2,  3,  we  obtain  the  solu¬ 
tion  matrices  [see  eqs.(108),  (109),  (II4)] 


1 

(t 

AS 

# 

1  t  • 
1 


with  a  , 


(270) 


1  0 

1 

- 

O’ 

# 

1  t 

1 

a  tf 

1 

a  t 

1  t 

1 

-  -  T 

i  th  *  e  and  Ji  ss  - t 

f  a 


1  a 

a  0»at  *  ® 

1  *  £ 

a  ii 

1! 

1 

1  * 

II 

a  -  <xt 

1 

a  t 

O’ 

AS, 

f 

O’ 

«  * 

1 

•  • 
•  li  a  a 


-t  «.  a-*  '  ; 

aith  and  s  _ £  i-L 

J- 


(271) 


(272) 


Here,  according  to  eqa(260),  we  have 

* 

Aa 

i  - 


C^t  ♦ 


nm. 

(273) 


Next,  we  determine  the  matrix  ®  for  the  transformation  [see  eq.(123)] 


which  brings  the  matrix  fi  to  the  Jordan  normal  form  ft  . 


This  then  yields  [see 


In  eq.(277),  it  must  also  be  noted  that 


j  =  1,  ii  =1,  i.e.,  j  -  ii,  Y  =  1  [see  eq. (195)3 


We  read  from  this; 

This  result  coincides  with  theorems  19  and  23 •  Analogously,  the  following  is 


valid  with  respect  to  eq.(278); 

j  =  2,  j  >  ii ,  Y  =  1 

Equation  (278)  indicates  that 


im 


tya  rtVft  w  »  *  *i«#/  m*i 

This  resxalt  also  coincides  with  theorems  19  and  23*  With  respect  to  eq.(279), 
we  have 


i«i .  i>s*^ » * 

“-y  «-W 


VtMtJ  V«*i 


which  again  coincides  with  theorems  19  and  23 • 

The  reciprocal  transposed  matrices,  conjugate  with  eqs.(277)  -  (279), 


- 


will  then  read 


(280) 


(281) 


ill 


• 

• 

# 

1 

• 

• 

t.cMtcz't'-f 

UJ-^  * 

-vr 

~t  1 


— 1 - 


1 

1 

-t  1 
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From  this,  one  can  read  from  the  last  rou 


(283) 


i.e.,  independent  of  j  in  accordance  with  the  proof  of  theorem  The  corre¬ 

sponding  periodic  solution  zctm  *  in  all  three  cases,  will  be  (see  theorem  13 ) 


A  ^  aln  t _ _ 

*[PD’  *  Dl  3  +  ain-'t  -  cos  t 


(284) 


Naturally,  the  function  Znou  can  also  be  calculated  in  accordance  with  the 


formula  (I63): 


♦  %  +  % 


(285) 


In  accordance  with  the  last  row  of  the  matrix  Y  [see  eq.(262)J,  we  have  the 


functions 


—  1+5  cos  t 


5  +  ain^t  -  cost 


and  Tp 


3+8in  t-C08t 


and,  in  accordance  with  the  power  row  of  the  matrix  $  [see  eqs.(270)-(272)l , 


the  functions 


cost  and  jYpj-  -2  .  sin  t, 


As  readily  verified,  this  will  also  yield  eq.(283). 

In  addition,  it  is  possible  to  read  from  the  representation  of  eqs.(277)  /172 
to  (279)  resp.  eqs«(280)  -  (282),  the  orders  mv  of  the  elementary  components, 
resulting  in  the  scheme 
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(286) 

which  coincides  accurately  with  the  scheme  of  the  minimal  orders  for  the  reso¬ 
nance  subcase  [see  eq.(185)]. 

For  our  special  case  (266),  the  quantity  u  =  l  in  accordance  with  eq.(267) 
is  the  exceptional  index.  Conversely,  a  computation  of  the  integrcil  [  see 


eq.(283)] 

/  hudt  *  VT  ^  0  J 


(287) 


shows  that  v  =  0  is  a  resonance  index. 


is  obtained  for  the  overall  minimal  power  orders  of  the  solutions  x(t)  for  j  = 

=  0,  1,  2,  3>  etc.: 

”  =  >^0  '  (289) 


The  general  solution  x(t)  of  eq.(269)  can  be  determined  in  two  different  ways: 
either  by  a  successive  integration,  starting  from  x(t),  to  eq.(288)  or  by  /I73 
constructing  a  particular  solution  in  accordance  with  the  method  of  variation  of 


the  constants  on  eq.(269)  and  adding  the  general  solution  of  the  conjugate  homo¬ 
geneous  differential  equation  which  must  be  taken  from  the  matrix  ^(t)  [see 
eqs.(277)  -  (279)3. 
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By  successive  integration  of  eq.(268)  (i  =  1,  2,  3),  we  obtain  the  solutions 


3  »  ij  x(t)  >  t  >  cost)  -  (2  Sint  —5  sin2t)  ♦ 
t  (1  -  cost)  -f  (-2^0^  Sint  -  coat  ♦  s^)/ 

J  “  2J  x(t)  "  ^  ^  ^  Sint  +  (3  cost  —  j  cos2t) 

2 

ly  +  t(a^  -  Sint)  ♦  cost  -  ^c^  sint'««2}^ 

j  -  3I  x(t)  .  ^  -  t  cost  ♦  (4  Sint  -  ^  sin2t)  ♦ 

t^  1  2 

TT  ^  *  *2^  *  (^c^cost  ♦  a,)* 


(290) 


(292) 


The  variational  method  for  the  constants,  using  the  argument 

x(t)  ■  ®x(t)  ♦  \(t)  (291) 

and  [see  eq.(20)] 

i 

x(t)  '  jjtjct)  •  /  a^t)  f(t)<lt  far  y-  o,n  (292) 

would  yield 

s  y  «Ut  f  1  f  t 

(293) 

xtn  s  t  (  (x’t)-  -  j  y 

f  t  V,  t 

which  coincides  with  the  first  equation  of  the  system  (290)  when  taking  into  /17/a 
consideration  that  ^Cq,  ^ci,  and°ci,  or  ai  are  arbitrary  constants.  Analog- 


(293) 


ously,  we  obtain  for 


°x(t)  -  Sint  t  -  4  Sint  +  J  8in2t  ♦  ®c^), 
1  S  1=2  . 

x(t)  •  2^5T  *  ^  sint)  +  (3  cost  —  ^  cos2t  + 

2  1  1;^  1 

4  sinS  -  5  8int.ain2t)  +  + 

t(^c<j  -  **Cjj  Sint)  cost  - 


(29Zt) 
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and  for 


i-3 


1  ^ 

Sint  ♦  02)* 

S  t2  1 

x(t)  •  2  '  7T  *  ^  '*'  ~  B 

^  t£  .  ♦  t  ("'c^  —  "*c^slnt)  ♦ 

“’c^  coat  -  -  ®c^)  slnt  ♦  "’cg 


®x(t)  -  -  co«t  t  -  4  Bint  ♦  i  aln2t  ■♦•  °c^), 

s  t5 

1^(t)  -  jC^  ♦  *  cost)  ♦  4  Bint  - 

^  ain2t  *  ^  aiii2t  coat  ♦ 

.  *3  1 

'c^(^  ♦  t  coat)  ♦  ®1^5T  *  ■*■ 


'e,t  ♦  'c, 


(295) 


C  i 

x(t)  ■  2  *  ^  ^  cost  ♦  (4  Bint  -  8in2t)  ♦ 

3  2 

|y  '’Cq  ♦  It  ^C-,  ♦  * 


Cc^  -  ®c<,)  coat  ♦  Cj. 

Again,  we  obtain  coincidence  with  the  two  last  relations  of  the  syston  (290),  /175 
Solutions  with  the  indicated  minimal  order  (289)  are  readily  constructed 
[see  eqs.(293)  -  (295)]: 

For  j  =  0,  we  only  must  put  (as  already  known)  ^  =  1  in  eq.(268). 

For  j  =  1,  the  selection  of  constants  in  eq,(290)  is  arbitrary. 

For  j  =  2,  we  put  ^  cq  =  while  the  other  constants  are  arbitrary. 

For  j  =3,  we  put  = - 1-  and  ^ci  =0,  while  the  other  constants 

are  arbitrary. 

If  we  start  from  a  fixed  solution  x(t),  power  orders  are  obtained  for  the 
sequence  of  the  functions  x(t),  x'  (t),  x”(t),  x''*(t),  etc,,  whose  patterns  can 
be  determined  in  accordance  with  eq,(290).  Let  us  consider  the  case  j  =  ?  and 
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select  a  definitely  determined  x(t),  resvilting  in  the  power  orders  of  the  deriv¬ 
atives  x^^^Ct)  for  four  different  cases  as  shown  in  Figs.l^  -  17  [see  eq.(295)J» 


In  all  these  cases,  the  following  is  valid: 

j  =  3,  ii  =  1,  Y  =  1  [see  eq. (195)3 

(ib~nb=l  [v  =  0  resonance  index,  see  eq.(  287)3 

J 

Case  It  ^  Co  f  0. 

Here,  we  have  u)i  =  3,  -ty  ^j-iy  ~2  [see  eqs.(20/4.)  and  (296)3. 
taining  pattern  corresponds  to  Fig. 7. 


Power  orders 


i 


Fig.l4 


Case  II:  '*■  ^  ■  0,  but  ''' ci  f  0, 

2 

Here,  we  have  uh  ^  2,  -ty  =2. 

The  following  pattern  corresponds  to  Fig. 8: 


H  «■*. 


Power  orders 


\vy _ 


#  1  <  W 

Fig.l5 


Case  III:  ^  Co  ^  =  0,  ^  ci  =0,  and  ^  Cg  arbitrary. 

Here,  we  have  (Uj  =  1,  -ty  =  =  Mn  (j  -  ii ,  )  =  Min.  (2,  l)  = 

eqs.(20i!,.)  and  (296)3. 


(296) 

The  per- 

/m 


1  [see 
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In  accordance  with  Fig. 6,  we  obtain  here 


Power  orders 


•  I  2  i»i  • 

Fig.  16 


Case  IV:  Cq  -1  =  0,  and  the  remaining  c  arbitrary. 


Here,  we  have 


“  ?*  ly"  J  “  ^4'“ 

y-  o  (see  eq.(220)),  i  -  i-  -  O  (*«  «9(''52)), 

y 

««y»n  =1  (see  eq.(296)), 

5  5°  ° 

l«j-i  »5-o»5  (»««  eq-C226)). 


(297) 


The  constant  ^  co  is  so  selected  that  [see  eq. (291)3  the  coefficient  of  t  van¬ 
ishes  in  the  siim  x"*(t)  =  °x"*(t)  +  ^x”’(t)  (see  theorem  24). 


«-=S-H 


Power  orders 


'\  V® 


..  V* 

*•.  _<*) 

— 

t  I  t  >•»  »-i»f 


Fig.17 


See  also  Fig. 10. 


In  this  case,  the  condition  (242)  must  be  satisfied.  We  will  check  this 


conr'-’’ +,ion;  the  only  remaining  item  is 


(298) 
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where  is  the  coefficient  of  the  highest  power  in  [see  eqs.(l22)  resp. 
(199)  and  (26)],  i.e.,  -J-  c<j  according  to  eq.(295),  whereas  is  the  con¬ 
stant  factor  of  the  highest  power  in  °x(t),  i.e.,  equal  to  according  to 
eq»(295)»  In  accordance  with  eq.(203)  resp.  (26),  the  coefficient  of  °x(t)  can 
be  computed  [see  eq.(279)]t 

\(t)  .  cost), 

o  7 

so  that,  necessarily,  cb  =  .  Consequently,  the  condition  (298)  reads  as 


follows : 


— ^  ^  ~  which  yields  ^  qj  =  1 


/m 


as  had  been  assumed  above,  for  the  case  IV. 

The  four  traces  from  Figs. 14  -  17  are  conpiled  in  Fig. 18,  together  with 
the  minimal  solution  of  Fig. 18  (the  latter  is  shown  as  a  solid  line). 


It  is  readily  demonstrated  that  all  possible  cases  are  covered  by  the 
cases  I  -  IV.  No  solution  x(t)  exists  which,  simultaneously  with  all  three 
derivatives  x*  (t),  x’'(t),  and  x”*(t),  would  have  the  corresponding  minimal 
order,  as  represented  in  Fig. 18  by  the  solid  line. 

As  a  second  exanple,  we  will  consider  the  differential  equation 


(299) 


.  2  cost 

+  3-2  sint 


-  r{t) 


(i  ^  o) 
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with  two  different  right-hand  sides :  /179 

(a)  f(t)  «  aint  (2  slot  -  5)t 

(b)  f(t)  •  (sint  -  1)  (2  Bint  -  3). 


The  general  solution  of  the  inhomogeneous  reduced  differential  equation  (j  =  O) 
has  the  following  form:  ^  ^ 


with 


y(2)<<=> 


(300) 


^  ^  A 

“  T^-lj(t)  •  ^  coat, 

yp5(t)  -  t  •  ♦^^jCt)  - 

'I  'I 

=  t  •  ^  cost  (y  aint  ♦  1), 
7(2)^^^  -  “  8i“t. 


(301) 


According  to  theorems  25  -  27»  as  already  expressed  in  the  notation, 

m  ■  2.  -  1,  -  1,  1-5  »*  nonexistent, 

'  ^  ^  (302) 

<  “l  -*‘2  ■  ®  » 

m  ^  cost,  ■J  sint  +  1,  “  sint. 


Fvirther,  for  the  case  (a)  we  calculate 

X  (t)  •  t(^  sint  -  cost)  ♦ 

sint  +  3  cost  +  ^  sint  cost  -  3) 

and,  for  the  case  (b),  2 

X  (t)  »  jj—  cost  -f  t  (sint  -  cost  +  3)  ♦ 

•tC^  cost  -  j  sint  ♦  j  sint  cost  -  3), 

By  means  of  the  matrix  g(t)  =  (^^“*(t)y  ,  the  following  is  obtainea: 

I*(1)  ^  *  2  si^t^-;  ^ 

.  ^f3  cost  -  sint  cost) 

♦  2  sint  -1 

A  A  —  S 

2  sint  -  5  ’ 


(303) 

/180 

(30a) 


(305) 


159 


For  the  case  (a),  this  will  yield 


1 


2 

♦  3  slot  *■  cos  t 
2  Sint  -  3 


tr  /*■ 

J  f(t)dt  ■  -  /  5  Bint  dt  -  o  , 


and 
<r  •; 


y  •>  1  exceptional  index 

tw 


J  ^pj(t)f(t)dt  »  J  Bint  ('l-f3  Bint  +  cob  t)dt 

V«  2  resonance  index 


and,  for  the  case  (b), 


tv 

**  -  -  3  J  (Bint-  1)dt 

0 


/ 


a6V^o,  y«  1  resonance  index 

fr 

f»  -  •  f  ? 

f  ^pj(t)f(t)dt  «  J  (sint  -  1)(1^5  Sint  +  cos  t)dt 

?  St 


4  and 


«o,  >*■  2  exceptional  index 

The  solution  Zcton  (t)  is  calculated  according  to  the  formula  (163 ): 


'foj  jTcu 


The  still  to  be  determined  functions  jTcii  >  defined 
tions  with  zero  mean  value  from  the  system  of  differential  equations 
in  our  case,  has  the  following  appearance: 


As  solutions,  we  obtain 


II 

A 

’»’o, 

'  t' 

Ij  - 

A 

* 

'  ^.7 

%  - 

X 

S  -••At’ 

Y 

'(f)  *  J 

Sint , 

0,  ^^2) 

■  -  cost 

(306) 


(307) 


(308) 

as  fvinc-  /181 
(97)  which, 


(309) 
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On  substituting  this,  together  with  eq.(305),  into  eq.(308),  we  obtain 


•  -  .  7  coat  ^  3  aint  coat 

ffafl  2  Sint  -  5 


(310) 


In  addition,  according  to  eqs.(l62)  and  (305),  we  have 


as  well  as 


2  Sint  -■ 


®C23  ■  'f[22 


1^3  Sint  4  cos  t 
2  Sint  -  5 


(311) 


(312) 


In  both  cases  (a)  and  (b),  the  exceptional  case  exists  for  v  -  0. 

Thus,  the  scheme  of  the  minimal  orders,  according  to  eq.(185),  will  be 


j  .  o  12  3  »  3 

o  1111  1 

2  2545J^1 

1  1111  1 


From  this,  in  the  case  (a),  we  obtain  for  j  -  0,  1,  2,  3»  •••  the  sequence  /182 


of  minimal  orders  [see  eq.(184)]: 


In  the  case  (b),  we  obtain  accordingly 


Then,  we  determine  the  general  solution  x(t)  for  j  =  4*  yielding: 
In  the  case  (a): 


(313) 


(314) 


x(t)  ■  -  Sint  -  cost)  +  Sint  ♦ 

+  |.C08  t  +  jj^sin2t)  ♦  5T  ®1  fr  *2  It 

A  ^  \  1  '1 

+t{  cost  ♦  CjjSint  +  ®/|*  5  cost  ♦ 


Sint 


In  the  case  (b): 
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5  4  2 

x(t)  cost  ♦  t  (-  5  Sint  -  cost)  ♦ 

^  3 

+  Sint  -  ^  cost  ♦  ^  sin2t)  ♦  5T  *  *1  ?T 
2 

t*2  It  *  ^  (“'c^  •  j  cost  ♦  Sj)  ♦  (-'’c^  Sint  ♦ 

^C2‘  j  cost  +  Sint  ♦  s^). 

In  the  case  (a),  four  subcases  must  be  differentiated: 

^Co  - 

^*0^  _  5  -  o,  but  O. 

-  5  •  s^  ■  o»  6ut  *2^®’ 

-  3  ■  s^  “  *2  ■  o» 

from  which  the  following  curves  for  the  power  orders,  analogous  to  Fig. 18,  are 


obtained: 


Power  orders 


Z182 


•>» »  4  p. 


‘  \1S\  \ 
\  • 

,  ^ 


•  f  <  s  i*f  r  wF-imt 


Fig.l9 


In  the  case  (b)  only  one  possible  case  exists  with  arbitrary  constants, 
yielding  the  pattern  shown  in  Fig. 20; 


Power  orders 


9  *  z  3  r  —»-/ 
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It  should  be  mentioned,  in  addition,  that  it  woxild  not  have  been  necessary 
in  the  above  examples  to  check  the  index  v  =  0  as  to  resonance  since  the  minimal 
order  nfe  can  never  be  greater  than  or  IH3 . 

Example  3;  As  the  third  exanple,  let  us  consider  the  differential  eqviation 


»  *  ainSt  ♦  2  coa2t  ^(d+2) 
sin2t  >  4  *  ♦ 


I  eoaft  -  sln2t  ^  gln2t  -  6  (j)  -f(t){j  io) 

T  sln2t  *  4  *  iln2t  +  A  *  ^ 


(315) 


with  the  right-hand  sides  /18L 

(c)  r(t)  «  coat  (4  2  Sint  cost), 

(d)  t(tj  «  (4  -f  2  Sint  cost). 


The  general  solution  of  the  inhomogeneous  reduced  differential  equation  has  the 


following  form: 

^(t)  -  x(t)  ♦  ''cg  ^(2)^^^ 

with 

<  ^  ♦tp3(t)-tco8t  ♦  Sint 

72<.t)  -  tgCt)  •  s^  -  s*, 

where,  in  the  case  (c). 


(316) 


(317) 


x*(t)  ”  5}^  cost  -  t  cost  ♦  J  Bint)  + 

•t(^  Sint  -  ^  cost  -  I  sin^t  - 
”5^  cos^t  -  ^  Sint  cos^t  -  j  sin^t  cost) 

and,  in  the  case  (d). 


(318) 


(t)  -  I  e^- 


t  cost  {—  Sint  +  w  sin2t  - 


Thus,  we  have  mi 
is  present;  for  v  =  1, 


C082t  -  Sin^t  -  (319) 

=  2,  ffte  =  1  in  which  case,  for  v  =  2,  the  principal  case 
no  ii  exists.  For  the  matrix  3(t)  =  ,  we  obtain 
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'*  -5  ^  .  sict  ♦  cost  ♦ 

* ( 1 ) *  ”  ^[13 *  "1 )  “  4+2  Sint  cost 

.  slnt  -  2  cost 

♦ — -  t 

4-»2  Sint  cost 

A  ^ 

z  *  *  sirxt  ^  cos^ 

4  ^  g  Sint  "coat  • 


z-i  •  ^  ® 


1  *003  t  .  e 

U  *  2  Sint  cosi 


-t 


(320) 


In  the  case  (c),  we  calciolate 


t* 


J  fCOdt  -  -  i  (sint  ■*  cost)  cost  dt  • 

*  -  TT  ^  o,  V-  *1  resonance  index 


/m 


(321) 


and,  in  the  case  (d). 


fw 


t-r 


f(t)dt  -  —J  (sint  +  cost)dt  »  o, 
y»  1  exceptional  index 


(322) 


For  the  orders  of  the  elementary  con?3onents  with  variable  j,  we  obtain  the  scheme 


■ 
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“2 
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1 

In  the  case  (c),  we  obtain  v  =  1  as  the  resonance  index.  For  establishing  that 
V  =  0  is  the  exceptional  index,  it  is  only  necessary  to  determine  Z([o])  which, 
according  to  eq.(l63),  yields 

^  ^  —2*2  sint  cost 

4  +  2  sint  cost 

which  means  that  v  =  0  is  the  exceptional  index.  According  to  eq.(18Z4,),  the 
following  values  are  ottained  for  the  minimal  orders: 

«  -  «  2,2,2,  etc.  (323) 

Thus,  in  analogy  to  Fig. 18,  the  pattern  shown  ^n  Fig. 21  will  be  obtained. 
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For  each  curve,  a  normal  solution  x(t)  can  be  defined.  /186 

Case  (d):  With  the  period  P  =  2n,  the  quantity  v  =  0  is  the  resonance 
index,  v  =  1  the  exceptional  index,  and  v  =  2  the  principal  index.  Since,  for 
the  reduced  differential  equation,  the  resonance  case  does  not  exist,  the  follow¬ 
ing  can  be  calculated  according  to  theorem  18: 

t*  iJt 

i  f(t)dt  -  J  x(t)dt  -  -4Tr^  o 

with  the  solution  x(t),  periodic  with  2n  [see  eq.(3l6)]  at  *Co  = - r-  and 

s 

^ Cl  =1.  For  the  minimal  orders,  the  following  values  are  obtained: 

m  -  »  0,1, 2, 5,  etc. 

For  the  power  oi-ders  with  two  case  differentiations,  the  pattern  of  all  possible 
cases  is  obtained.  (In  the  dotted  curve,  (t'l  =  1  has  been  assumed  while,  in  the 
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dashed  curve,  the  minimal  solution  uji  =  0  has  been  assumed. ) 


Example  Ui  In  the  differential  eqviation 

Sint  +  »  *  cost  ^(j+2)  ^ 


2  *  cost 


Sint  A  _ 2 - 

»  2  +  cost  2  ♦  cost 


2  ♦  cost 

-  2(2+cost) 


•with  j  ^  0  which,  for  j  =  0,  has  the  general  solution 


x(t)  -  ^  cost  +  J  t  Sint  +  2  cost  -  2  + 


+  ”'c  (t  Sint  1)  +^c^sint  t  cost  +  C2*cost 


we  obtain 


‘tjU  •  Sint,  ^21”  cost,  ®’ 

^  A 

«  2  mit  «  1  ond  nonexistent, 


im 

(324) 


(325) 


(326) 


Here,  our  scheme  for  the  minimal  orders  reads 
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According  to  eqs.(l64)  and  (166),  we  calculate 


^£13"  IciJ" 


-  sint 
t2+cost ) 


md 


2  cost  ■»  1 
2(2  ♦  cost! 


(327) 


This  will  yield 


tr 


and 

tr  ^  . 


r(t)dt  =  2 

)  f(t)dt  =■ 


I 


$ 

/» 

; 


Sint  (it  .  o 
(2  cost  +  1)dt 


f 


2  r#  o. 


(328) 


Consequently,  the  exceptional  case  is  present  for  v  =  1,  while  the  resonance 

case  occurs  for  u  =  2.  Thus,  the  overall  minimal  order  [see  eq.(18/4)]  will 

be  m  =  Max(mv )  =  2  for  all  j.  Here,  it  is  of  no  inportance  whether  v  =  0  is  /188 

(  V  r  e  ».) 
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the  resonance  index  or  not  (see  case  IV  on  p.l08).  All  solutions  x(t),  no 
matter  how  large  j  might  be,  show  the  same  behavior  with  respect  to  the  sequence 
of  the  power  orders  of  x,  x* ,  x",  etc.  This  will  yield  the  pattern  shown  in 
Fig. 23. 

1  Power  orders 


tit 


Fig.23 


As  a  final  exanple,  let  us  consider  the  general  differential  equation  with 

constant  coefficients  ,  ACn-D  Arn-i-^^ 

t  [xj  s  x'‘  X*-”  J  x  -  ilt) 

f  o),  (329) 

where  f(t)  is  to  have  the  period  P.  It  can  be  readily  proved  that  all  solu¬ 
tions,  periodic  with  P,  of  the  reduced  homogeneous  differential  equation  are 
functions  with  a  mean  value  of  zero.  If,  namely,  a  function  y(t),  periodic 
with  P,  is  inti^sduced  into  the  homogeneous  reduced  differential  equation,  all 
terms  which  are  multiplied  by  the  derivatives  y* ,  y”,  ...»  have  a  mean 

value  of  zero.  Consequently,  also  the  last  term,  aa-jy(t),  must  have  a  zero 
mean  value. 

It  follows  from  the  above  statement  that  the  functions  /18 

-  7(y)(t)  (  V-  1,...,  f)  [see,  for  example,  eq.(250)] 

have  a  mean  value  of  zero.  We  must  now  prove  that  all  functions  $p,(t)(M.  =  (v)  + 
+  1,  ...,  [v];  V  =  1,  ...,  p)  must  have  a  mean  value  of  zero. 


For  exan^sle,  the  differential  equation,  homogeneous  xio  eq.(329),  is  solved 


by  means  of  the  argument  y(t)  -  This  is  assumed  to  yield  the  character¬ 
istic  equation  *  /♦ 

„  (330) 

is-  * 

with  pairwise  different  values  cti ,  ,,,,  ag  where,  in  addition  to  a  complex  o',,, 

also  the  conjugate- complex  root  occurs.  Then,  the  general  solution  of  the 

differential  eqxiation  homogeneous  to  eq.(329)  will  read 


t  at*  t 

y(t)  -  ”  '^it> » 

where  Pv(t)  represents  an  arbitrary  polynomial  of  the  degree  mv  -  1: 


(331) 


fA— 1)! 


(332) 


The  differential  equation,  homogeneous  to  eq.(329),  has  solutions  with  the 

M. 

period  P  exactly  when  the  solutions  ...,  org  of  the  characteristic  equation 
/  X  2ni 

(330)  contain  whole  multiples  of  — ^ — «  Let  us  assume  that  this  is  true  for 
the  values  ofi ,  nti'  while  the  remaining  a  are  no  whole  multiples  of  — - — . 

p 

We  repeat  here  that  ckv  f  0  since,  otherwise,  the  differential  equation  homogene¬ 
ous  to  eq.(329)  would  have  to  have  a  constant  solution  which  contradicts  the 
stipulation  of  an-j  f  0.  For  the  eigenvalues  =1,  ...,  p,  ...,  s),  the 


following  bound  [see  (Bibl.l)  eq.(20)] 

- 1  < 

at  first,  is  rot  valid.  Generally,  we  put 


Y.  i” . . . 9- 


(333) 

/190 

(33L) 


According  to  eqs.(33l),  (332),  and  (33A)»  the  first  row  of  the  fundamental 
solution  matrix  ?)(t)  has  the  form 


>  ‘{a: 


(335) 


with 
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336) 


[see  eq.(95)]. 


From  this,  we  obtain  directly 


■.b>2,(y- . ,?). 


(337) 


That  ('^  ^  •••>  ^  mean  value  of  zero  was  determined  above; 

naturally,  the  remaining  =  0  also  have  a  mean  value  of  zero. 

The  above  statements  yield  the  following  theorem: 

Theorem  28:  In  the  case  of  a  differential  equation  with  constant  coeffi¬ 
cients,  no  iv  exists. 

Consequently,  the  matrix  fi'  automatically  is  in  the  Jordan  normal  form 
[see  eqs,(llZj,),  (123)];  Tor  the  transitory  transformation  matrix,  5  =  ®  is 
valid# 

Consequently,  the  scheme  for  the  minimal  orders  has  the  following  form: 

/191 


(338) 


In  general,  at  low  values  of  j  and  if  the  resonance  case  is  present  for  the 
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reduced  differential  equation,  the  resonance  subcase  cannot  be  of  any  signifi¬ 
cance  for  the  index  v  =  0.  However,  as  soon  as  j  is  sufficiently  large,  up  will 
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predominate  over  all  other  orders  mi,  mi  gmd  it  becomes  of  inportance  to 

determine  whether  v  =  0  is  a  resonance  index  or  not.  For  this  we  need  the  solu¬ 
tion  Zitoij  (t),  periodic  with  P,  of  the  adjoint  homogeneous  differaitial  equa¬ 
tion,  which  already  is  defined  for  j  =  1, 

^  [®IE<30  ]  *  ~  it  ^  [*irou]  ■  °  •  (339 

First,  we  will  demonstrate  that  Ztcoij(t)  must  be  constant.  For  this,  we  start 
from  eq.(329)  with  f(t)  of  zero  mean  value.  By  j  integrations  and  making  use 


of  the  auxiliary  theorem  in  Section  3,  ve  obtain  an  equation 


L  .  a,  ..... 


P  (t) 


(340 


with  an  F(t)  of  zero  mean  value.  Since  x(t)  has  the  power  order  of 

m  -  -^*91  (A,  ) 

the  same  power  order  will  be  valid  also  for  x(t).  Since  we  have  nt  =  j  in  ac¬ 
cordance  with  eq.(338),  it  follows  that  v  =  0  is  the  exceptional  index  for  j  > 
>  in.  Consequently,  it  then  is  necessary  that 


4T  so 


The  assumption 


Z192 
(342 


with  zitoi)  of  zero  mean  value  and  constant  k  would  lead  to  a  contradiction 


with  eq.(34l)  if 


F(t)  =  ziion  (t) 


is  selected.  Therefore,  it  is  obvious  that  in  eq.(342),  we  have  ztcojj  (t)  =  0. 
The  constant  k  is  determined  in  accordance  with  eq.(l74)  as 

*r>v.Ti(t)  -  k  -  _ 1 _  . 


From  eq.(179)  xe  then  find  that  v  =  0  is  the  resonance  index  for  the  differ- 
tial  eqviation  (329)»  provided  that  f(t)  has  a  mean  Vcilue  differing  from  zero 
and  will  be  the  exceptional  index  if  the  mean  value  of  f(t)  is  zero  (see  also 
footnote  on  p.35)» 
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